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The Fermi Gamma-ray Space Telescope, now in its fourth year after 
launch, continues to make important discoveries and establish new results 
in all directions: pulsar astronomy, cosmic-ray physics, AGN and black-hole 
astrophysics, galactic astronomy, gamma-ray bursts (GRBs), limits on dark 
matter and Lorentz invariance violations, 7-ray astronomy of the Sun, moon, 
and Earth, etc. In this chapter, I survey results at medium energy 7 rays, 
from some tens of MeV (at energies above nuclear de-excitation 7-ray lines) 
to sa 100 GeV where the ground-based 7-ray Cherenkov detector arrays be- 
come more sensitive. As shorthand, Fermi and medium-energy 7-ray astron- 
omy is referred to here as "GeV astronomy," and ground-based Cherenkov 
7-ray astronomy ^ 100 GeV as "VHE astronomy." 

The Fermi results already provided considerably more material than could 
be presented in the nine lectures that I gave on this subject at the Saas-Fee 
school on "Astrophysics at Very High Energies," held 15-20 March 2010 in 
Les Diablerets, Switzerland. Happily, though, Professor Lars Bergstrom gave 
a brilliant series of lectures that covered dark matter, so the absence here of 
extended discussion on dark matter and new physics in GeV astronomy reflects 
Prof. Bergstrom's better capabilities to address this subject. My lectures and 
this book chapter are therefore restricted to astrophysical and astroparticle 
sources of GeV radiation rather than to 7 rays with origins in exotic particle 
physics and dark matter. 

Even while the school was in progress, news appeared of a new type of 7- 
ray emitter of GeV photons that was identified with the symbiotic binary Nova 
V407 Cygni. This extraordinary system reveals an explosive shock evolving on 
a timescale of days to weeks, rather than the hundreds of years for supernova 
remnants (SNRs). I showed the 2010 March 18 ATEL [l] announcing the Fermi 
result in my last lecture, capturing the real-time recognition of a new type 
of galactic 7-ray source triggered by a thermonuclear explosion on a white 
dwarf fed by its binary red giant's wind. The discovery has been published by 
the Fermi Collaboration [5], and is already triggering a new line of research 
strongly tied to MeV line astronomy and white-dwarf physics (cf. citations 
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to [2])- So these lectures and this write-up capture an early glimpse of the 
state of knowledge of astrophysical sources of it 100 MeV and GeV radiation 
obtained with the Large Area Telescope (LAT) on Fermi at « 2 - 3 years into 
the mission, weighted by the extragalactic interests of the author. 

Alongside the LAT on Fermi is the Gamma-ray Burst Monitor (GBM), 
sensitive to GRBs and bright transients in the 10 keV - 30 MeV range. 
This review can only briefly mention important GRB results made with 
the GBM — by itself and with the LAT — and related GRB science employ- 
ing Swift, INTEGRAL, and other detectors. Indeed, multiwavelength sci- 
ence is value-added science, and the possibilities to uncover the underlying 
physics of the powerful compact systems that are at the heart of high-energy 
astronomy are multiplied by radio/microwave/sub-mm/IR/optical/UV/X- 
ray/MeV/TeV/neutrino/gravitational-wave data correlated with the GeV 
window, now observed with unprecedented clarity due to the LAT on Fermi. 

The GeV field is in full discovery mode, not only due to Fermi but 
also thanks to AGILE, an EGRET-like sentinel of bright 7-ray transients, 
and to ground-based VHE observatories, particularly HESS, VERITAS, and 
MAGIC Q The i^F^ spectral energy distributions (SEDs) based on simulta- 
neous and overlapping data sets are providing valuable information about 
Galactic sources, blazars and radio galaxies, and starburst and normal galax- 
ies. GeV astronomy with Fermi is still in the midst of an active phase as 
Fermi accumulates data and increasing time makes faint sources visible and 
detection of rare cosmic transients more likely. 

This chapter will be divided into sections that follow the course of lectures 
delivered at the Saas-Fee course. Though now somewhat out-of-date, these 
lectures can be found on my website in ppt format]^ The topics of the lectures 
and the corresponding sections of this chapter are: 

1. GeV instrumentation and the GeV sky with the Fermi Gamma-ray Space 
Telescope 

2. First Fermi Catalog of Gamma Ray Sources and the Fermi Pulsar Catalog 

3. First Fermi AGN Catalog 

4. Relativistic jet physics and blazars 

5. 7 rays from cosmic rays in the Galaxy 

6. 7 rays from star-forming galaxies and clusters of galaxies, and the diffuse 
extragalactic 7-ray background 

7. Microquasars, radio galaxies, and the extragalactic background light 

8. Fermi observations of GRBs 

9. Fermi acceleration, ultra-high energy cosmic rays, and Fermi at 2 



^ AGILE: Astro-rivelatore Gamma a Immagini LEggero (Gamma-ray Light Imag- 
ing Detector); HESS: High Energy Stereoscopic System, in Namibia. VERITAS: 
Very Energetic Radiation Imaging Telescope Array System, in Arizona. MAGIC: 
Major Atmospheric Gamma-ray Imaging Cherenkov Telescope in La Palma, Ca- 
nary Islands 

heseweb. nrl.navy. mil / gamma / ~dermer /default . htm 
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Besides a discussion of the results of the Fermi Gamma-ray Telescope, 
I also include here some high-energy astrophysical theory essential for anal- 
ysis of 7-ray data that builds on the research presented in my book with 
Govind Menon, Ref. "High Energy Radiation from Black Holes: Gamma 
Rays, Cosmic Rays, and Neutrinos," published by Princeton University Press 
in 2009. The book itself is focused on theory rather than observation, and 
develops the hypothesis that the most energetic and powerful radiations in 
nature are made by particles energized through Fermi acceleration processes 
in black-hole jets powered by rotation. 

It is not possible, even at this early stage of the Fermi mission, to ade- 
quately summarize all the results from Fermi. But together with the accompa- 
nying lectures and book, this article provides a broad overview of some recent 
astrophysical advances in the Fermi era of 7-ray astronomy. 

1 GeV Instrumentation and the GeV Sky with the 
Fermi Gamma-ray Space Telescope 

1.1 Historical introduction 

The year 2010 represents a highpoint in high-energy astronomy. Astronomical 
observatories at multiwavelength electromagnetic, neutrino, cosmic-ray, and 
gravitational wavebands are operating and collecting data. X-ray astronomy 
has sensitive pointing X-ray telescopes, Chandra and XMM Newton, deployed 
in space to observe catalogued sources. Broadband X-ray and soft 7-ray ob- 
servatories like Swift, INTEGRAL and Suzaku are available to measure the 
X/7 spectra of compact objects. The Rossi X-ray Timing Explorer (RXTE) 
continues to operate and provide a monitor of the brightest X-ray sources 
in the sky. Already mentioned were AGILE and the ground-based 7-ray air 
Cherenkov telescopes. The Cherenkov Telescope Array (CTA) consortium [4] 
is planning to start building as early as 2014 with array completion as early 
as 2018. Lacking at the moment, however, is an operating all-sky water-based 
Cherenkov telescope successor to MILAGRO in the 1 - 100 TeV range. This 
gap will soon be filled by the High Altitude Water Cherenkov (HAWC) ex- 
periment [5] on the Sierra Negra mountain near Puebla, Mexico, which uses 
300 tanks rather than a single pond as utilized by MILAGRO. 

The Pierre Auger Observatory, located in the Mendoza province of Ar- 
gentina and covering an area the size of Rhode Island, had its third data 
release in 2010 ^ , giving the spectrum and composition of ultra-high energy 
cosmic rays (UHECRs) with energies E ^ 10^^ eV. The IceCube Neutrino 
Observatory, most sensitive to astrophysical neutrinos with energies of 100 
TeV - 100 PeV (« lO^^* - lO^^ eV), has deployed all 86 of its strings in the 
latest austral summer (December, 2010), and has developed the DeepCore 
subarray that is sensitive to lower energy, « 10 GeV - 100 GeV, neutrinos. 
LIGO, the Laser Interferometry Gravitational- wave Observatory, hoping to 
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detect gravitational radiation from coalescing compact objects, is operating 
at design sensitivity. Development of an order-of-magnitude more sensitive 
Advanced LIGO has been approved, with completion expected for 2017 or 
thereafter. The NASA Laser Interferometry Space Antenna (LISA) is sup- 
ported in the recent Astro-2010 study, though ESA is developing a separate 
space-based gravitational wave facility. Likewise, Constellation-X has evolved 
into IXO/ATHENA0 The large-area X-ray timing mission RXTE will soon 
be ended, with the ESA Large Observatory for X-ray Timing (LOFT) mission 
taking its place. Here in the US, a focusing hard X-ray telescope, NuSTAR 
(Nuclear Spectroscopic Telescope Array), in the 5-80 keV range, and the 
Gravity and Extreme Magnetism Small Explorer, GEMS, a NASA mission to 
study X-ray polarization of astrophysical sources in the 2-10 keV range, will 
soon be launched. 

Brief history of GeV astronomy 

The progress of GeV astronomy in the range from w 10 MeV to w 100 GeV 
followed a period of remarkable advances starting over 40 years ago that culmi- 
nated with the launch of Fermi. Prior to the Fermi-LAT, the most important 
detectors and some of their achievements in the development of medium- 
energy 7-ray astronomy are the following: 

• 1967-1968, OSO-3, the Third Orbiting Solar Observatory, carried a Cherenkov 
counter experiment sensitive to > 50 MeV 7 rays, of which 621 were de- 
tected [7]. It discovered one source, the extended 7-ray emission of the 
Milky Way. 

• 1972-1973, SAS-2, the Small Astronomy Satellite-2 [S], a spark chamber 
experiment sensitive to 7 rays with energies between « 30 MeV and 200 
MeV (and an integral flux above 200 MeV). It detected « 8000 celestial 
7-rays, making the first 7-ray identifications of the Crab and Vela pulsars, 
Geminga (7I95 -I- 5, then unidentified), and the Cygnus region, and an 
association with Cygnus X-3 was suggested. A north-south asymmetry in 
the Galactic 7-ray plane emission was noted and attributed to the massive 
stars in the Gould belt. An isotropic 7-ray background radiation was first 
reported from analysis of SAS-2 data HH] ■ 

• 1975-1982, COS-B, the Cosmic ray Satellite (option B), a magnetic-core, 
wire-matrix spark chamber sensitive to 7 rays with energies from ^ 30 
MeV - « 5 GeV), with an effective area of « 50 cm^ at 400 MeV pi] . 
Its orbit resulted in a large and variable background of charged particles. 
During its lifetime, it detected ~ 200,000 7-rays, with the COS-B Car- 
avane Collaboration announcing the discovery of 25 sources, most along 
the Galactic plane. These included 2CG 135-1-01, now identified with LSI 



^ ESA's Cosmic Vision International X-ray Observatory/Advanced Telescope for 
High ENergy Astrophysics, sensitive to photons with energies 0.1 - 40 keV 
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+61° 303, and the first extragalactic source of ^ 100 MeV 7 rays, 3C 273 

• 1991-2000, EGRET, the Energetic Gamma Ray Experiment Telescope on 
the Compton Gamma Ray Observatory, was a spark chamber experiment 
with large, ^ 1200 cm^ effective area between 200 MeV and 1 GeV, accom- 
panied by excellent background rejection [13|. During its 9 year mission 
lifetime, it collected > 2 x 10® 7 rays and discovered that blazars and GRBs 
are luminous 7-ray sources. Because of the importance of this experiment 
in the development of GeV astronomy, we describe EGRET in more detail 
below. 

• 2007-, AGILE, with a Gamma Ray Imaging Detector (GRID) and small 
calorimeter giving 7-ray sensitivity from 30 MeV to 50 GeV, roughly com- 
parable to EGRET. Along with the GRID is the accompanying 18 ~ 60 
keV super- AGILE hard X-ray survey instrument [2]. The first AGILE 
GRID catalog consists of 47 > 100 MeV 7-ray sources with significance 
> 4(T from data taken between July 2007 and June 2008 [151 . The brightest 
sources in the catalog, e.g.. Vela, Crab and, during periods of outburst, 3C 
454.3, exceed integral photon fluxes of 2 x 10~® pli(> 100 MeV)/cm^-s, 
while the weakest are fainter than 50 x 10^® pb(> 100 MeV)/cm^-s. 

The EGRET experiment on CGRO 

The Compton Gamma Ray Observatory, or GGRO, was a pioneering 7-ray 
space observatory (5 April 1991 - 4 June 2000) consisting of 4 instruments, 
OSSE, the Oriented Space Scintillator Experiment (sensitive from « 50 keV 

- 10 MeV); BATSE, the Burst And Transient Source Experiment (w 20 keV 

- few MeV); COMPTEL, the Compton Telescope (w 800 keV - 30 MeV); 
and EGRET, the Energetic Gamma Ray Experiment Telescope (w 30 MeV 

- 10 GeV). egret's detection method utilized a gas-filled spark chamber 
that tracked the 7 ray after it converted to an electron-positron pair by pair- 
production interactions with nuclei in the thin Ta foils in the gas-filled spark 
chamber. Directional information was obtained by time-of-flight coincidence 
with a scintillator array in the lower spark chamber assembly. Charged par- 
ticles were vetoed by an anticoincidence shield and total 7-ray energy was 
measured with EGRET's total absorption scintillator counter, the TASC, con- 
sisting of crystal scintillators and photomultipliers [16] . 

The Third EGRET (3EG) catalog [17^, made from data accumulated be- 
tween 1991 April 22 and 1995 October 3, consists of 271 sources of > 100 
MeV emission, including a Solar flare in 1991, the Large Magellanic Cloud, 
five pulsars, one probable radio galaxy, namely Centaurus A, and 66 high- 
confidence detections of blazars. The catalog also lists 27 lower confidence 
potential blazar detections, and contains 170 unidentified EGRET sources 
lacking associated sources. Five GRBs were detected with the spark chamber 
on EGRET. Figure 1 shows a skymap of the 3EG sources. 
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The EGRET Field of View (FoV), defined roughly by the solid angle within 
which the effective area is greater than 1/2 of the on-axis effective area, was 
« 0.5 sr, or w 1/24*'* of the full sky. EGRET operated in a pointing mode, 
and targeted one region of the sky for two weeks, representing ~ 10^ s after 
time for Earth occultation and time spent in the South Atlantic Anomaly are 
subtracted. The Point Spread Function (PSF) at 100 MeV was « 5.7°, with 
the PSF improving roughly as « E^-^/^ [Tgl [T3]. The first 18 months of the 
mission were devoted to a full-sky survey. 

Important results from analysis of the full 6 years of data (degradation 
of spark chamber gas led to smaller apertures after 4 years into the mission) 
include the following: 

1. Diffuse extragalactic background with intensity w 1.5 keV/cm^-s-sr, cor- 
responding to an energy density ~ 10^^^ erg/cm"^; 

2. The galactic flux is as much as 50 x brighter than the extragalactic flux, 
and much softer, but not as soft as expected if the 7 rays are formed 
by secondary nuclear production by Galactic cosmic rays with the same 
spectrum as observed locally (the so-called "EGRET excess"); 

3. Typical fluxes of EGRET sources are between ~ (10"^ and 10"^) ph(> 
100 MeV)/cm^-s, with a typical 2-week on-axis limiting flux at w (15 - 
25) xlO-s ph(> 100 MeV)/cm2-s; 

4. Galactic sources, including young radio pulsars; 

5. Sources with significant flux variability. 

Temporal variability is an essential characteristic of GRBs and blazars, Solar 
flares, and now V407 Cyg. The large FoVs of the LAT and Swift, and larger 
still with BATSE and GBM, are crucial for study of 7-ray transients. 

Point source sensitivity of EGRET 

The best description of results leads to the question of units in GeV astronomy. 
Because medium-energy and high-energy 7-ray astronomy is challenged by 
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limited signal counts, an integral photon flux is the most natural unit. For 
the EGRET experiment, units of 10"^ ph(> 100 MeV)/cm^-s are suitable, as 
this value roughly separates signal-dominated and noise-dominated detection 
in EGRET, as we now show. Units 10~® ph(> 100 MeV)/cm-^-s are more 
suitable for Fermi sources. 

As measured with EGRET, the diffuse, or at least unresolved intensity 
of photons at 100 MeV, is measured to be el^ « 1.5 keV/cm^-s-sr at = 
emeC^ « 100 MeV. Writing 

eh = 1.5 i — r^-' , (1) 

eioo cm^-s-sr 

with photon number index and eioo = 100 MeV/mgC^, implies 

dN{> e) _ r°° e'l,, _ 1.5 x 10"^ e _^ ph( > 100 MeV) 
dAdtdQ meC^e'^ ~ I eioo cm^-s-sr 

(2) 

EGRET measured a hard spectrum, with F^ « 2.1 [19J, whereas LAT mea- 
sures a softer spectrum, with F^ w 2.4 120 , and with a 100 MeV intensity only 
60% as large as the EGRET intensity^ For the EGRET intensity, therefore, 
dN{> eioo) /dAdtdfi « 1-5 x 10"^ ph(> 100 MeV)/cm2-s-sr. 

The EGRET PSF at 100 MeV, as previously noted, is « 5.7°. In compar- 
ison, the LAT PSF (68% containment radius) at 100 MeV is « 3.5°, with the 
PSF dropping to ~ 0.6° for 1 GeV photons, and w 0.15° for 10 GeV photons 
(for conversion in the thin layers) [21]. The EGRET PSF at 100 MeV repre- 
sents about 7r(5.7°)^ sa 100 square degrees, or « 1/400*'* of the full sky. Thus 
the flux from each patch corresponding to the EGRET PSF is « 5 x 10^*" 
ph (> 100 MeV)/cm^-s. A 7-ray source is signal dominated for EGRET when 
its flux is > lO"*' ph (> 100 MeV)/cm2-s. By contrast, it is noise dominated 
when its flux is < 10"^ ph (> 100 MeV)/cm2-s. 

The time At needed to accumulate 100 photons with EGRET, consisting 
of 50 signal S and 50 background B photons, and to give a detection at the 
« Sj^lB = ha level, is found through an expression for the integral photon 
flux F(> e) of the source. For a source at the level of 

p_ ^ -^(> '^lOo) /ON 

^ 10-8 ph(> 100 MeV)/cm2-s ' ^ ' 

the number of detected photons is w lO^^F-g x Z\i x 1000 cm^, so that « 50 
ph can be detected from a bright source at the level of F_g ~ 10^ during an 
EGRET observation period of At « one day (during half this time, the Earth 



A number index F^ = 2.41 ± 0.05 and intensity normalization /(> 100 MeV) 
= (1.03±0.17) X 10~^/cm^-s-sr for the intensity spectrum of the isotropic diffuse 
7-ray background are measured from the first year Fermi data ,20,. The lower 
diffuse extragalactic flux measured with Fermi compared to EGRET is partially 
but not entirely due to resolving more point sources out. 
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is occulted). A nominal 2- week observation period with At « 10^ s on source 
gives a 5a limiting sensitivity for sources with i^_g ?s 2 x 10^^ ph(> 100 
MeV) /cm^-s. 

These sorts of arguments can be used to estimate the time needed to make 
a detection and resolve temporal variability with EGRET, LAT, and counter 
detectors with broad FoVs, including neutrino telescopes. The Fermi LAT be- 
comes noise dominated at much lower flux levels than EGRET, with sources 
regularly detected at the F_8 < 1. Note that the background is smaller at 
higher energies as a result of the smaller PSF, but the smaller flux usually 
reduces significance for a given observation time except for hard-spectrum 
sources. A better approach for characterizing detection significance is likeli- 
hood analysis [22] (described below), but the detection significance can be 
simply estimated as outlined above ^5] . 

The above estimates apply to high galactic latitude, \b\ > 10° sources 
where the 7-ray sky is dominated by extragalactic sources and unresolved 
isotropic 7-ray background. At lower galactic latitudes, the diffuse 7-ray emis- 
sion from cosmic-ray interactions with gas and dust makes source detection 
more difficult. Subtraction of the diffuse emission and nonuniform and vari- 
able background requires a Galactic model for cosmic-ray/gas interactions. 
The sensitivity of source detection to background model is seen in the alter- 
nate 3EGR analysis of the EGRET data by Casandjian & Grenier [21] . They 
do not confirm 107 3EG sources, most in the vicinity of the Gould Belt, and 
find 30 new sources in the full 9-year data set. 

Long after EGRET's effective lifetime had expired (though BATSE, OSSE, 
and COMPTEL were still collecting valuable data), CGRO was deorbited into 
the Pacific Ocean in 2000. The launch of INTEGRAL in October 2002, Swift 
on November 20, 2004, and AGILE on April 23, 2007, helped fiU the 7-ray 
gap leading to Fermi. 

1.2 Fermi Gamma- ray Space Telescope 

The Fermi Gamma-ray Space Telescope consists of two major detector sys- 
tems, namely 

1. The Large Area Telescope (LAT), a pair-conversion telescope sensitive 
between 20 MeV and ^ 300 GeV, with the higher energy limit a result 
of the vanishing small detection rate given Fermi's ~ 1 m^ aperture. The 
Fermi-LAT has opened the previously unexplored sa 10 GeV - 100 GeV 
window, as self-vetoing due to particle shower backsplash in the anti- 
coincidence detector may have reduced EGRET's effective area above k, 5 
GeV [25|. Fermi nominally operates in a survey mode, and given its large 
FoV, scans the entire sky every 3 hours. 

2. The Gamma-ray Burst Monitor (GBM), sensitive in the 8 keV - 40 MeV 
range, consisting of 12 NaT detectors sensitive between 8 keV and 1 MeV, 
and 2 BGO detectors sensitive in the 0.15 MeV - 40 MeV range. The 
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scintillator detectors surround the base of the LAT, and view the entire 
unocculted sky which, at the nominal 565 km altitude of Fermi, represents 
« of the full sky. 

Some of the science questions that Fermi was designed to answer are: 

• How do supermassive black holes in Active Galactic Nuclei (AGN) create 
powerful jets of material moving at nearly light speed? What are the jets 
made of? 

• What are the mechanisms that produce GRB explosions? What is the 
energy budget? 

• What is the origin of the cosmic rays that pervade the Galaxy? 

• How does the Sun generate high-energy 7 rays in flares? 

• How has the amount of starlight in the universe changed over cosmic time? 

• What are the unidentified 7-ray sources found by EGRET? 

GLAST becomes Fermi 

After a number of delays, a Delta H 7920-H(eavy) rocket carrying the Gamma- 
ray Large Area Space Telescope — GLAST — payload was launched from Cape 
Canaveral Air Station on 2008 June 11 at 12:05 pm EDT into a 565 km 
altitude circular orbit with 25.6° inclination and a 96 minute period. GLAST 
completed a 60 day checkout period by early August, and released its first- 
light image on 26 August 2008, when it was renamed the Fermi Gamma-ray 
Space Telescope, after the Italian- American physicist Enrico Fermi (1901 - 
1954). 

The first-light image, based on only 4 days of observation, already reveals 
dozens of sources to the eye, with the blazar 3C 454.3 almost as bright as 
the Vela pulsar, which is the brightest persistent 7-ray source in the sky 
(F_8 = 1060). Limiting fluxes of EGRET sources are at the level of F^s w 20 
- 30 for a two-week observation, with a corresponding all-sky flux limit of 
F_g « 15 - 30 for a year of observing (given the FoV of EGRET). For a 
source with a flat vF^ SED, or a number index F-^ = —2, Fermi reaches 
F_8 ^ 1 in one year over the entire high Galactic latitude \h\ > 10° sky. Due 
to its energy-dependent effective area and PSF, limiting fluxes are strongly 
dependent on both source and background spectrum. Hard-spectrum sources 
with number index « — 1.5 are detected with comparable significance as soft 
spectrum sources, but at integral photon fluxes as low as F_g w 0.1. 

The LAT operates in a nominal scanning mode whereby the spacecraft 
rocks about the zenith. The rocking angle was equal to 39° in the first part 
of the Fermi mission, and then increased to 50° after 3 September 2009. The 
larger rocking angle gives a more uniform exposure, but with the loss of data 
due the increased fluorescence 7-ray emission from cosmic-ray bombardment 
of the Earth's atmospherej^The LAT observes the entire sky every two orbits. 



^ This is also called, colloquially and inaccurately, "7-ray albedo," or just "albedo". 
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or « 3 hours, by rocking north and south of the zenith on ahernate orbits, with 
each point in the sky receiving « 30 min exposure during this time. Onboard 
triggering has also been enabled for the Fermi spacecraft to autonomously 
slew. As of October 2011, more than 30 autonomous repoint requests (ARRs) 
have taken place, resulting in 5-hr pointing mode observations in response 
to bright GRBs detected with the GBM (the length of the ARR repoint 
was reduced in the third year of the mission). Several dedicated Targets of 
Opportunity (ToO) pointings have been executed, including one to 3C 454.3, 
two to the flaring Crab, one to Cyg X-3, and two Solar pointings]^ A regular 
schedule of nadir pointings to look for 7 rays from terrestrial 7-ray flashes, or 
TGFs, is underway. 

The Fermi/GLAST LAT Collaboration is an international organization 
originally comprised of institutions in France, Italy, Japan, Sweden, and the 
United States. There are about 120 full Collaboration members, and a total 
of about 200 scientific members, including affiliated scientists, postdocs, and 
students. The principal investigator of the project is Professor Peter Michelson 
of Stanford, and the project scientist during the development, commissioning, 
and early science phase was Dr. Steven Ritz, now professor at UC Santa Cruz. 
Dr. Julie McEnery of NASA's Goddard Space Flight Center is now Fermi LAT 
project scientist. A review of the first-year LAT results is given by Michelson, 
Atwood, and Ritz |26j , and the detailed instrument description can be found 
in the paper by Atwood et al. (2009) [21 . 

In recognition of the importance of the Fermi mission to high-energy as- 
tronomy, Peter Michelson, Bill Atwood, and the Fermi Gamma Ray Space 
Telescope LAT team were awarded the 2011 Rossi Prize of the High En- 
ergy Astrophysics Division of the American Physical Society, "for enabling, 
through the development of the Large Area Telescope, new insights into neu- 
tron stars, supernova remnants, cosmic rays, binary systems, active galactic 
nuclei, and gamma-ray bursts." UCSC physicist Bill Atwood was also recently 
awarded the 2012 W.K.H. Panofsky Prize in Experimental Particle Physics 
by the American Physical Society. 

LAT instrument description 

The LAT is a pair-conversion telescope composed of 16 tracker-converter tow- 
ers, a calorimeter, an anticoincidence shield, and an electronics system to 
filter signal from background. The conversion of 7-ray photons to electron- 
positron pairs takes place when a 7 ray produces a pair in one of the 16 
high-Z tungsten planes in each tower. The photon direction is measured with 
a precision silicon-strip tracker with a 228 /i pitch (separation between strips), 
totaling 8.8 x 10^ channels distributed over 18 tracker planes in each of 16 
tracker modules. The Si tracker follows the e^-e" pairs to give direction, and 
the photon energy is measured in a by a cesium-iodide calorimeter in a ho- 
doscopic assembly that images the shower, allowing for additional rejection of 

^ fermi.gsfc.nasa.gov/ssc/observations/timeline/posting/ 
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non-photon-like events. A segmented anticoincidence detector (ACD) rejects 
background of charged cosmic rays. Segmentation of the ACD helps prevent 
self-vetoing at high energies. 



Y I incoming gamma ray 
I 
I 




electron-positron pair 



Fig. 2. Schematic of the Large Area Telescope on Fermi (from Atwood, et al., 2009 

my 

Table[T]gives the performance characteristics of the Fermi LAT. Depending 
on the source type, different classes of events are adopted to maximize sensi- 
tivity. Three hierarchical analysis classes based upon analysis cuts and event 
selections are generally used for Fermi analysis. For the most stringent diffuse 
event class used to identify 7 rays, the on-axis effective area increases from 
1500 cm^ at 100 MeV to « 8000 cm^ at 1 GeV, and is roughly constant at 
higher energies. The transient event class is used for the study of GRBs and 
brief transients where the amount of background can be greater due to the 
smaller time over which data is taken. The source class, intermediate to the 
other two, is where the residual instrumental plus environmental background 
is comparable to the extragalactic diffuse emission measured with EGRET. 

LAT instrument response 

Even though two tracker towers and three calorimeter modules of the LAT 
were beam tested, the original science tools provided for instrument response 
were the Pass 6_v3 response functions based on extensive GEANT4 Monte 
Carlo simulations of the satellite that followed secondary particles and photons 
through the digitization and software filters to determine — assuming tracker 
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Table 1. Performance Characteristics of the Fermi LAT 1211 



Parameter 



Value or Range 



Energy range 

Effective area at normal incidence (60° off-axis) 
100 MeV 
1 GeV 
10 GeV 

Energy resolution {la, on-axis): 

100 MeV - 1 GeV 

1 GeV - 10 GeV 

10 GeV - 300 GeV 
Single photon angular resolution (68% containment radius): 

100 MeV 

1 GeV 

> 10 GeV 
Field of View (FoV): 
Timing Accuracy 
Event readout time (deadtime) 



20 MeV - 300 GeV 

3700 (700) cm^ 

< 8300 (2900) cm^ 

< 8400 (3100) cm^ 

15% - 9% 
15% - 9% 
8.5% - 18% 

3.5° 
0.6° 

< 0.15° 
2.4 sr 
300 ns 
26.5 /is 



azimuthal symmetry — effective area, energy uncertainty, and PSF (or con- 
tainment radius). The event reconstruction, filters, and different event classes 
are described in |21) . 

The P6_vll corrections to the P6_v3 response functions using a few months 
of on-orbit data [27] showed that the earlier response functions gave smaller 
angles for 68% and 95% containment radii above w 5 GeV than given by the 
on-orbit calibration. The PSF is best fit with a King model profile of the form 

/Ki„g(0,a,7) = (l--)(l + -^r^ (4) 

where 6 is the angle between the incident, "true" photon direction and the 
reconstructed direction, and g is a fitting parameter. The wings of the PSF 
follow a power-law rather than exponential behavior, and are well fit with the 
sum of two King functions. 

The first-year analyses and the IFGL, and second- year analyses leading 
to the 2FGL, generally employ the P6 and P7 instrument response func- 
tions, respectively. Updated instrument performance that improves on pre- 
flight Monte Carlo and muon calibrations by using inflight corrections to the 
instrument response functions are found at 

Utilization of Fermi-LAT to energies as low as 30 MeV is possible with 
the LAT Low-Energy (LLE) technique [55]. This approach was developed for 
transients, such as Solar flares and GRBs, where a direction is known and 
a short time window can be defined. During this period, less discriminating 



www.slac.stanford.edu/exp/glast/groups/canda/lat_Performance.htm 
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event selections are used, the analysis includes all photons with reconstructed 
arrival directions within 20° of the target, and the criteria for > 100 keV 
vetoing in the ACD is relaxed. 



1.3 Energy, flux, and luminosity 

Consistent with the notation of [3 , we define the dimensionless photon energy, 
in units of electron rest mass, as 

Flux density Fi, is usually reported in units of Jansky (1 Jy = 10"^'^ erg/cm^- 
s-Hz), so that the quantity vF^^ is an energy flux F (units of erg/cm^-s, or 
Jy-Hz, noting that 10^" Jy-Hz — 10^^'^ erg/cm^-s). The luminosity distance 
for a steady, isotropically emitting source is defined so that the energy flux 
F is related to the source luminosity (erg/s) according to the Euclidean 
expression 

If (/>(e) is the measured spectral photon flux (units of photons cm^'^ 
s^^e^^). then I'F^ — rrieC^ 4){e) . Henceforth we use the notation 

fe = (7) 

for the uFjj flux. From the definitions of F and /e, 

F= de^- . (8) 

Considering eq. (|6|, the luminosity radiated by a source between measured 
photon energies ei and £2, or between source frame photon energies ei(l + z) 
and £2(1 + 2), is therefore given by 

('€.2 /"In 62 

L,[ei(l+z), £2(1 + 2)] = 47r4TOecM de e (jj^e) = Andl d(ln£)/, . (9) 

J ei Jlnei 

Eq. ^ shows that if the i^F^ spectrum is flat with value , corresponding 
to a photon flux 0(£) cx £^^, then the apparent power of the source over one 
decade of energy is « (InlO)Lo ~ 2.30Lo, where Lq = And'j^f^. 

From these relations, one obtains the mean apparent isotropic 7-ray energy 
release 

_ 4ndlF{At) 
(l + z) ■ 

for a source at redshift z that releases average energy flux F during observing 
timescale At. The apparent versus the absolute energy releases and lumi- 
nosities depends on the jet structure and variability behavior. For a steady. 
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two-sided top-hat jet with uniform emission within angle 9 < 9j of the jet 
axis, the absolute luminosity 

Labs — fbLiso , (11) 

where the beaming factor /f, = 2x 27r diij/in = 1 — For 6j = 0.1 (5.7°), 

fb = 1/200, whereas if 6, = 0.01 (0.57°), fb = 1/20, 000 5 x 10"^ What is 
reported is the apparent isotropic luminosity L*^iso = ^T^d^F, with absolute 
luminosity implied by arguments for the jet opening angle and beaming factor. 



Variability information 

The Schwarzschild radius of a black hole of mass M is 

Rs = ^ 3.0 X 10^ f — j cm = 10"^ Mg pc, (12) 

defining Mg = M/{10^ Mq). Variations in the source flux by a large factor 
{ ^ 2) over a time scale At must, from causality arguments for a stationary 
source, originate from an emission region of size R ^ cAt/{l + z). Incoher- 
ent superpositions of emission from larger size scales and from regions that 
are not in causal contact would usually (though not always) wash out large- 
amplitude fluctuations. For high-quality data from bright flares, large ampli- 
tude variations in source flux on timescale At would, from this argument, 
imply a black-hole mass 

M,5^^. (13) 

1 + z 

This relation is even preserved in systems with relativistic outflows, unless 
specific mo del- dependent provisions are made on the method of wind dissi- 
pation. Variability timescales far shorter than the light-crossing time across 
an « IO^Mq black hole have been measured in the TeV BL Lac objects PKS 
2155-304, Mrk 501, and Mrk 421, discussed in more detail below. Another in- 
teresting fiducial is the ratio of the Eddington luminosity to the light-crossing 
time across a Schwarzschild black hole [29 , namely 

= 4^ ^ 2.5 X 10^3 erg . (14) 

We take this criterion as separating the extreme universe from the moderate 
universe. During intense 3C 454.3 flaring activity [3D], the ratio of the apparent 
isotropic luminosity and the source variability timescale strongly violate this 
limit, making this an extreme event. 
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log [E(eV)] 

Fig. 3. Spectral energy densities in intergalactic space of various radiations, in- 
cluding the CMB, the infrared (IR) and optical, X-ray, 7 ray, and the extragalactic 
cosmic ray energy density. Also shown is the energy density of cosmic rays measured 
near Earth; the transition energy between the galactic and extragalactic component 
remains uncertain. 

Extragalactic background light (EBL) 

In intergalactic space, the energy density of the EBL is dominated by that of 
the cosmic microwave background radiation (CMBR), with present tempera- 
ture TcMB = 2.72 K and energy density ucmbr = 0.25 eV/cm^ = 4 x 10^^^ 
erg/cm'^. The intensity of the EBL at infrared frequencies is difficult to mea- 
sure directly because of foreground radiations, such as Galactic electron syn- 
chrotron radiation and zodiacal light scattered by dust in our Solar system. 
Photons observed at 10 - 100 GeV interact primarily with EBL photons at op- 
tical and UV energies. The energy density of the dust and stellar components 
of the EBL is ~ 10% of the CMBR energy density at the present epoch. 

Figure |3] shows the energy density of photons in intergalactic space and 
cosmic rays in outer space in the Solar cavity near Earth. Fermi measurements 
PU] give a lower intensity for the diffuse extragalactic 7-ray background (EGB) 
than EGRET data [111131], probably due to EGRET miscalibration above « 5 
GeV EgplThe cosmic-ray particle spectrum is modulated at low, ~GeV/nuc 



The spectral energy density of the isotropic background radiation field at 100 
MeV as measured with EGRET is 
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energies by changing Solar activity. Cosmic-ray source origin changes from a 
Galactic to extragalactic origin at an energy scale currently under debate. The 
case for an extragalactic origin of the UHECRs can be made by comparing 
Galactic size with the Larmor radius 

of an ion with energy Lorentz factor 7, atomic charge Z and atomic mass 
A. The characteristic magnetic field B = 10~^i?^G G of the Milky Way is a 
few /xG in the disk and probably much less in the Galactic halo. 



1.4 Limits to the extreme universe 

The largest luminosity is obtained if the entire rest mass energy of the object 
with mass M is transformed into energy on a light-crossing time scale for the 
gravitational radius of the object, so that 

i^a. = - ^ = 3.6 X 10^9 s-i (17) 

[32I [33] . The most luminous blazar source yet detected with the Fermi L AT 
is 3C 454.3, with apparent isotropic 7 ray luminosity reaching ss 2 x 10^° 
erg/s. The most luminous GRBs detected with the Fermi-LAT, namely GRB 
080916C and GRB 090510A, reached w lO^^ gj.g/g 



2 Fermi Gamma-ray Source Catalogs and Fermi Pulsars 

Now we take a global view, and look at the Fermi Large Area Telescope First 
Source Catalog (IFGL) [33] as a catalog, or systematic collection of objects, 
before focusing on 7-ray emission from neutron stars and the first Fermi pulsar 
catalog |3S]. The IFGL catalog, taken with 11 months of data between 2008 
August 4 and 2009 July 4, expands the 3 month bright source list (BSL, or 
OFGL) [36], just as the First LAT AGN Catalog (ILAC) [37], the subject of 
the next lecture, expands the LB AS (LAT Bright AGN Sample) [35]. The 
2FGL and 2LAC build on the IFGL and ILAC, respectively. The IFGL and 
2 FGL catalogs are described in this section. 

Source detection significance for the IFGL and 2FGL catalogs is based on 
the likelihood test statistic (TS) . The likelihood ratio compares the likelihood 
of a null hypothesis with the likelihood of the model. According to Wijk's 



Because the eJe spectrum is essentially flat, the total energy density over the 
100 MeV - 100 GeV energy range gives a factor In 10^ ?^ 6.9; thus the EGRET 
energy density of the EGB is nearly 5 orders of magnitude below the present 
CMBR energy density. 
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theorem, —2 times the logarithm of this ratio approaches for large x^- The 
likelihood C is the probability of observing counts in pixel i,j given the 
expected number of counts A^- in pixel i,]. Poisson statistics gives 

^ I — ■ (18) 

Dropping a term independent of model, the logarithm of the likelihood is 

log £ = ^ n,j log(A„ ) - ^ Xij , (19) 

ij ij 

and TS = — 21og(£/£o) ^ X^j where Cq is the likelihood of the null hypoth- 



esis for the given data set. Eq. (19 1 gives a prescription for calculating TS 
when analyzing 7-ray data. Mattox and colleagues '22] originally applied this 
method to EGRET data. 

The Fermi catalogs are based on a limiting test statistic, which is TS > 25 
in the IFGL, 2FGL, ILAC and 2LAC, and TS > 100 in the BSL and the 
LB AS. This introduces systematic effects that should be taken into account 
when interpreting the IFGL. For example, a cut on TS produces integral 
flux thresholds strongly dependent on the source spectral hardness. Because 
of the steeply falling Galactic background at high galactic latitudes, high- 
latitude hard spectrum sources will have a large TS for much smaller integral 
fluxes than a soft-spectrum source detected at the same significance level. 



2.1 First Fermi catalog of gamma-ray sources: IFGL 

The IFGL comprises 1451 sources with TS > 25, corresponding to a signifi- 
cance ^ 4cr, and represents 21.22 Ms (or w 73.3% livetime) of data. Most of 
the deadtime loss is due to passage through the SAA South Atlantic Anomaly 
(« 13%) and to readout deadtime (9.2%). The exposure is uniform to within 
a factor « 1.25 between north and south and the detections are based on in- 
tegrated data rather than shorter bright periods or flares. An improved back- 
ground and calibration model is used in the IFGL compared to the OFGL. 

The IFGL is not a flux-limited survey. Although the exposure is rela- 
tively uniform — an excellent feature of the GLAST design in nominal scan- 
ning mode — there are large differences in integral photon fluxes of sources 
depending on source location and the strong Galactic diffuse emission. Conse- 
quently, the ILAC catalog is drawn from sources above 10° latitude where the 
Galactic diffuse intensity is lowj^Any conclusions are weakened when source 
identifications are not complete, for example due to the limits of a counter- 
part catalog used to make associations, or to uncertainty in distance or lack 
of redshift, as this incompleteness can bias results. 

^ The Galactic diffuse flux at 10° < lb| < 20°, averaged over longitude, is still a 
factor of 2 - 3 greater than the extragalactic diffuse at higher latitudes [39]. 
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The IFGL catalog gives source location regions and association criteria 
defined in terms of elliptical fits to the 95% confidence regions, and power- 
law spectral fits as well as fiux measurements in five energy bands for each 
source. In addition, monthly light curves are provided. Firm identifications 
with sources found in other astronomical catalogs are based on correlated 
variability, e.g., rotation or orbital period in the case of pulsars and X-ray 
binaries, or morphological similarities with counterparts at other wavelengths, 
as in the case of SNRs. 

For the catalogs and association criteria used, 630 of the sources in the 
IFGL are unassociated. Due to the sensitivity of the results to the model of 
interstellar diffuse 7-ray emission used to model the bright foreground, 161 
sources at low Galactic latitudes towards bright local interstellar clouds are 
flagged as having properties that are strongly dependent on the Milky Way 
gas model. 

Table 2. LAT IFGL and 2FGL Source Classes [34ll40] 



Description Designator Assoc. (ID) Assoc. (ID) 







IFGL 


2FGL 


Pulsar, X-ray or radio, identified by pulsations 


psr (PSR) 


7(56) 


(83) 


Pulsar, radio quiet (LAT PSR, subset of above) PSR 


24 - 


25 - 


Pulsar wind nebula 


pwn (PWN) 


2(3) 


(3) 


Supernova remnant 


snr" (SNR) 


41 (3) 


62" (6) 


Globular cluster 


glc (GLC) 


8(0) 


11 (0) 


High Mass X-ray binary 


hxb (HXB) 


0(2) 


0(4) 


Micro-quasar object: X-ray binary (black hole 


mqo (MQO) 


0(1) 


0(1) 


or neutron star) with radio jet 








Nova 


nov(NOV) 




0(1) 


BL Lac type of blazar 


bzb (BZB) 


295 (0) 


428 (7) 


FSRQ type of blazar 


bzq (BZQ) 


274 (4) 


353 (17) 


Non-blazar active galaxy 


agn (AGN) 


28 (0) 


10 (1) 


Radio galaxy 






10 (2) 


Active galaxy of uncertain type 


agu (AGU) 


92 (0) 


257 (0) 


Normal galaxy 


gal (GAL) 


6(0) 


4(2) 


Starburst galaxy 


sbg (SBG) 


2(0) 


4(0) 


Seyfert galaxy 


sey (SEY) 




5 (1) 


Unassociated 




630 


576 + 1"^ 


Total 




1478'' 


1873" 



"Indicates a potential association with a SNR or PWN. 
779 -I- 630 -I- (69) = 1409 -I- (69) = 1478. Greater than 1451 because of multiple 
class assignments. "Some of the 62 sources may also be associated with PWNe. 
"^576 unassociated plus one with uncertain class. 
"1169 + 577 + (127) = 1746 + (126) = 1873. 



The principal source classes found in the IFGL and listed in Table [2] are 
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• Sun, moon, Earth 

• 3 high mass X-ray Binaries 

• Rotation-powered and milhsecond pulsars 

• Supernova remnants 

• Globular clusters 

• > 600 blazars 

• 28 Non-blazar AGNs and Radio Galaxies 

• Dozens of AGNs of uncertain type 

• 3 Starburst Galaxies 

In addition, hundreds of GRBs and more than a dozen LAT GRBs had been 
detected by the time of publication of the IFGL. 

Association depends on the underlying catalog and the degree of confidence 
assigned in the probability of association. For blazars, four catalogs with sub- 
stantial overlap were used, giving 689 IFGL associations with sources found 
in at least one of these catalogs. A fuller discussion of the AGN catalog is 
deferred to the next lecture. Though no extended radio lobes of radio galaxies 
were reported in the IFGL, the lobes of Centaurus A have since been imaged 

m- 

Regarding other source classes in the IFGL, 41 Fermi sources are asso- 
ciated with SNRs or non-pulsed 7-ray emission, and three are sufficiently 
physically extended that their morphology counts as an identification (W44 
W51C, and IC 443; see Lecture 5). Fermi reported three high-mass X-ray 
binary systems in the IFGL, namely LS 5039, LSH-61°303, and Cyg X-3 
(Lecture 7), and a fourth was discovered later. In the IFGL, it was unclear 
whether 7 rays were made by massive O stars, but since then 77 Carinae has 
been been established as a 7-ray source [42J. Ten BSL candidates out of 205 
in the OFGL, do not show up in the IFGL, illustrating the highly variable na- 
ture of some 7-ray sources (though note that all these sources are found in the 
Galactic ridge, |^| < 60°, where background modeling is especially important). 
On top of this, 630 sources were unassociated in the IFGL. The unidentified 
Fermi sources are clustered towards the plane of the Galaxy where diffuse 
background Galactic emission and higher source density makes associations 
more tentative, though there are still many high-latitude unidentified Fermi 
sources. 

2.2 Second Fermi catalog of gamma-ray sources: 2FGL 

The 2FGL catalog [40J was released at the time of writing. This source catalog 
was derived from data taken during the first 24 months of the science phase 
of the mission, which began on 2008 August 4. Source detection is based on 
the average flux over the 24-month period. The 2FGL includes source location 
regions and fits to model spectral forms. Also included are flux measurements 
in 5 energy bands and light curves on monthly intervals for each source. Twelve 
sources in the 2FGL are found to be spatially extended. 
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The 2FGL contains 1873 sources detected and characterized in the 100 
MeV to 100 GeV range, of which 127 are firmly identified and 1170 are re- 
hably associated with counterparts of known or Ukely 7-ray source classes. 
Ahhough the diffuse Galactic and isotropic models used in the 2FGL analysis 
are improved compared to the IFGL catalog, caution flags for 162 sources 
indicate possible confusion, given the uncertainty in the underlying diffuse 
model. Table [2] lists the number of sources of various types in the 2FGL. 

Some important improvements compared to the IFGL catalog are: 

1. The 2FGL catalog is based on data from 24 months of observations. 

2. The data and Instrument Response Functions use Pass 7 event selections, 
rather than the Pass 6 event selections used in the IFGL. 

3. The 2FGL employs a new, higher-resolution model of the diffuse Galactic 
and isotropic emissions. 

4. Spatially extended sources and sources with spectra other than power laws 
are incorporated into the analysis. 

5. The source association process has been refined and expanded. 

2.3 Fermi pulsars 

In the IFGL, 56 pulsars are identified by their 7-ray pulsations, and another 7 
associations are based on expectations from pulsar catalogs, e.g., E/d? rank- 
ing. Six IFGL sources are associated with pulsar wind nebulae (PWNe) that 
lack known pulsars, and 8 IFGL sources are associated with globular clusters 
illuminated, most likely, by the superposition of millisecond pulsar (MSP) 7- 
ray emissions. The number of Geminga-like pulsars that lack detectable radio 
emission, has grown by a large factor. There are 83 7-ray pulsars, now all 
identified, in the 2FGL. 

A catalog of 46 7-ray pulsars is presented in the First LAT Pulsar Catalog, 
Ref. [35j. This catalog includes 16 "blind-search" pulsars discovered by search- 
ing for pulsed 7-ray emission at the position of bright LAT sources. Pulsed 
7-ray emission from 24 known pulsars were discovered using ephemerides de- 
rived from monitoring radio pulsars, of which 8 are MSPs. The remaining 6 
7-ray pulsars were known previously. 

EGRET pulsars 

The Crab and Vela pulsars were known prior to CGRO, and Geminga was 
known as a bright 7-ray point source. The pulsar nature of Geminga was 
only established by detection of X-ray pulsations in ROSAT data early in the 
EGRET era gS]. The EGRET pulsars, with EGRET fluxes above 100 MeV 
and 2FGL fluxes between 1 GeV and 100 GeV, ordered by brightness, are 
hsted in Tabled 
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Table 3. 7-ray Fluxes of EGRET Pulsars [T7l[35ll40] 



Pulsar 


Period (ms) Age (kyr) 


EGKEl 


Pulsar Catalog 


; 2FGL 




P 


P/2P 


F-8 


F-8 


F(l - 100 GeV)" 


0833-45, Vela 


89.3 


11.3 


834.3 ± 11.2 


1061 ±7.0 


135.8 ±0.4 


J0633+1746, Geminga 


237 


340 


352.9 ±5.7 


305.3 ±3.5 


72.9 ±0.3 


0531+21,'' Crab 


33 


1.25 


226.2 ± 11.2 


209 ±4 


18.3 ±0.15 


1706-44 


102 


17.6 


111.2 ±6.2 


149.8 ±4.1 


19.1 ± 1.7 


1055-52 


197 


540 


33.3 ±3.82 


30.45 ± 1.7 


5.0 ±0.09 


1951+32^= 


39.5 


110 


16 ± 2 


17.6± 1.9 


2.1 ±0.07 


1509-58,'= Circinus 


88.9 


150 




8.7 ± 1.4 


1.45 ±0.08 




"Also in 


units of 10"* 


ph cm~^ s~^ 







''Associated with SN 1054 
''Pulsars not reported in the 3EG [4|] 

PSR 1951+32, located in the confusing Cygnus arm region, is associated 
with the SNR CTB 80 discovered by Kulkarni et al. [JS]. It has a phase- 
averaged lyFi, spectrum rising c>c e°^^ and peaking at « 2 GeV, with a slower 
than exponential decline at higher energies |46j . 

The Circinus pulsar was detected with COMPTEL but not EGRET, and is 
unusual in having an inferred polar magnetic field Bp — 15.77 TG, compared 
to ~TG (lO^'^ G) fields for the others. Geminga is unusually close, at = 160 
pc [37], whereas the others are more likely to be at ~kpc distances. 

Elementary pulsar physics 

Over 1900 pulsars are known today, mostly through radio surveys. Pulsars are 
rapidly rotating, highly magnetized neutrons star. The neutron stars them- 
selves are formed by explosions made by the collapsing cores of massive stars 
(core-collapse supernovae), or through accretion- induced collapse of white 
dwarfs. Misalignment of the magnetic axis and rotation axis makes in the sim- 
plest configuration a rotating dipole-field geometry. Emission beamed along 
favorably defined field lines naturally makes a wide range of pulse structures 
for different obliqueness and inclination and gap sizes. The additional range 
of parameters associated with period and period derivative, along with poorly 
understand radiation physics, allows for purely empirical and kinematic pulse 
profile fitting. 

Neutron stars are predicted to have masses near I.4M0 and radii ^ 15 
km. The two properties of the pulsar that can be very precisely measured are 
the period P and period derivative P. Besides these observables, theory gives 
the mass of the neutron star, Mjqs^ £^nd its radius Rns |48j . The uncertain 
equation of state of nuclear matter determines whether neutron stars with 
given masses can exist, whether collapse to a black hole occurs, or whether 
other degenerate quark phases exist. 
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Electron degeneracy pressure in a white dwarf cannot support the gravi- 
tational force of a degenerate core exceeding the Chandrasekhar-mass Mq = 
IAMq. Either through nuclear burning or accretion, neutron stars can be 
formed. Neutron stars with masses M > Mq can test black-hole forma- 
tion theory, but the uncertainty in orbital inclination usually makes for large 
uncertainties in the neutron-star mass. By means of radio timing measure- 
ments of the Shapiro dela^p"] of the binary millisecond pulsar J1614-2230, a 
1.97 ± O.O4M0 pulsar mass was recently measured [5D], which rules out some 
exotic hyperon and boson condensate equations of state. 

Noting that velocity v = f2 x R, and the angular speed J7 — 2tt/P, then 
the light-cylinder radius at which the speed v — c is 

Rlc = Pc/2n . (20) 

Even the rotation of a simple misaligned dipole field leads to unusual geometric 
features, for example, the footprint on the neutron star surface of the field lines 
that open into the light cylinder [5T] . 

An expression for the surface polar magnetic field can be obtained by 
equating the rotational spindown energy loss rate with the magnetic dipole 
radiation power. The former is 

dE... d n_j^,\i_n^^p^ (21) 



dt dt \2 J P3 

with the moment of inertia / depending on the mass and equation of state of 
neutron star matter. This implies a characteristic age 

T = P/2P , (22) 

for P much longer than the initial spin period. The magnetic dipole power 
can be estimated from the Poynting flux of a dipole radiation at the light 
cylinder, namely. 



- (4^i?Lcc) -2^Ns[^^)coc^. (23) 



Thus 



Bns ^VPP. (24) 
The magnetic field i?LC at Rlc is therefore 

/ , • \ 1/2 

/ 247r'i/P \ , , 
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The Shapiro delay is a regular general-relativistic change in the light travel time 
from the pulsar as the radio photon travels through the gravitational field of the 
binary system; see |49j . 
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We can imagine a time-independent configuration for a spinning, mag- 
netized neutron star if we wait long enough until all current flows have 
asymptotically relaxed to their steady-state values. This defines the condi- 
tion of a force-free magnetosphere, where the vanishing of the Lorentz force 
F = Q[{v/c) X B + E)] = implies the existence of a magnetosphere filled 
with plasma with Goldreich- Julian 52] density 

Pulsar models build on these elementary concepts. For generic polar cap 
models, charge depletion in the polar field lines that open to the light cylinder, 
compared with the Goldreich-Julian density, generates strong electric fields 
that induce vacuum breakdown and current flow. In outer gap models, strong 
fields are generated at the gaps created by the surface deflned by the condition 
fl ■ B = 0. The primary 7-ray production and attenuation processes are 
curvature radiation, Compton radiation, and magnetic pair production. Slot 
gap models study a more sophisticated realization of the pair-starved region 
and how a slot gap forms along field lines with different acceleration potentials. 
Besides discriminating between magnetospheric models, pulsar studies help 
disentangle the geometry of a pulsar. For more on pulsar models, see |53l 1511 



Properties of Fermi pulsars 

Two different approaches are taken to discover pulsars, either by period folding 
with data from a pulsar established at other frequencies, or to perform a blind 
search for the pulsation. In the first method, the timing parameters of pulsars 
already known at radio or X-ray frequencies are used to search for evidence 
of 7-ray pulsations in the 7-ray data. For the blind-search technique, spectral 
power is calculated from time-tagged 7-ray photon event data by searching 
through likely values of P and P values. This is extremely computationally 
intensity, and not feasible for MSPs found in binary systems. More compu- 
tationally efficient techniques compute the arrival time differences of 7 rays, 
which cluster at values corresponding to multiples of the pulse period |55j . 
For 7-ray astronomy, with limited photon statistics, only the time-differencing 
technique for blind searches is effective at finding pulsars. 

The first blind-search detection with Fermi was a pulsar with period P — 
0.317 s in the shell of the SNR CTAl, coincident with unidentified X-ray and 
3EG JOOlO-l-7309 7-ray sources [M]- The spin-down power of this pulsar is ~ 
4.5 X 10^^ erg s~^ , and the pulsar's inferred age is 14,000 yr, which is consistent 
with the SNR age. Of the 46 Fermi pulsars reported in the First Fermi Pulsar 
Catalog, 16 of them were found using blind-search techniques in sources that 
are positionally coincident with unidentified EGRET sources and supernova 
remnants |57] . These 16 pulsars are all young and highly magnetized, with 
inferred magnetic fields between « 1 - 10 TG. 
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Fig. 4. Vela light curves at optical, X-ray, and 7-ray energies [58], binned to 0.01 of 
the pulsar phase. The main peaks PI, P2 and P3 are labeled in the top right panel. 
The bottom left panel shows the 8-16 keV RXTE light curve [Spl along with the 
radio pulse profile (dashed lines). At lower right, the 4.1 - 6.5 eV FST/STIS NUV 
light curve [60] is shown. 



The rotational energy-loss rates of pulsars in the First Fermi Pulsar Cat- 
alog range from ^ 3 x 10'^^ erg s^^ to 5 x 10'^® erg s~^, with the young, 
highly magnetized pulsars typically having rotational energy-loss rates exceed- 
ing ^ 10^^ erg s~^. Comparing with the phase-averaged apparent isotropic 
7-ray luminosity implies efficiencies for the conversion of spin-down power to 
7-ray energy in the range from ~ 0.1% to « 100%. About 75% of the 7-ray 
pulses have two peaks, but a third emission structure, P3, shows up in the 
Vela pulse profile and moves to later phases with increasing photon energy 
[55] . As Fig. |4]shows, the main peak, PI, which is dominant from at least UV 
to GeV energies, becomes less intense compared to P2 at multi-GeV energies. 

The Fermi LAT phase-averaged spectrum of 7-ray pulsars can be well fit 
by a generalized exponentially cutoff power law, given by 
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Fig. 5. Phase-averaged spectrum for 0.1 < i?(GeV) < 60 [HI]- 



N{E) cx eiip[-{E/E,f] , (27) 

with hard photon number indices generally in the range 1 ^ F-y ^ 2, and 
cutoff energies Ec between « 1 and 5 GeV. This form has been found to apply 
to all types of pulsars, whether radio or 7-ray selected, or normal or millisecond 
pulsars. For Vela itself, as seen in Fig. [5] = 1.38 ± 0.03, E^ = 1.36 ± 0.15 
GeV, and b = 0.69 ± 0.02 (quoting statistical errors only) [ST]. The sub- 
exponential (6 < 1) rather than super-exponential (6 > 1) cutoff that would 
be produced if the 7 rays were made deep within the polar cap of the neutron 
star magnetosphere, and the detection of emission at tens of GeV, is evidence 
against a polar-cap model. A phase-averaged spectrum fitted with equation 



(27) is the composite phase-resolved spectra that individually can be fit by 
simple exponential behaviours |611 146| . Significant variations of cutoff energy 
with phase are needed to reproduce the phase-averaged spectra. 



Millisecond pulsars and globular clusters 

Pulsed 7-ray emission was detected with the Fermi LAT from J0030-I-0451, 
making it the first firm detection of a MSP in 7 rays although there 
was a marginal EGRET detection of PSR J0218+4232 ^ that has been 
confirmed with the LAT (55] . Nine months into science operations, the Fermi- 
LAT Collaboration reported 8 7-ray MSPs [SI] , establishing a new population 
of 7-ray emitters. As noted above, it is not computationally feasible to per- 
form blind searches for binary MSPs, which introduces too many possibilities 
into the timing solutions. This makes it harder to know what fraction of the 
unidentified sources are MSPs, or if there are any radio-quiet MSPs. 

Globular clusters are quasi-spherical stellar systems consisting of hundreds 
of thousands of old (ages of « 10"'^'' yr) metal-poor stars. The high-stellar den- 
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sity at the centers of these systems makes for unusual multi-body encounters. 
According to the most common scenario, binary neutron-star systems in low- 
mass X-ray binaries can be spin up by accretion over long time spans to 
periods of a few ms, provided that the polar surface magnetic field is ^ 10^ 
G, which is ~ 3 orders of magnitude smaller than the typical TG fields of 
young 7-ray pulsars. Even though they have very different magnetic field and 
rotation rates, the 7-ray spectra of MSPs and young pulsars are similar, and 
well-represented by a hard power-law with a modified exponential cutoff given 
by equation (27 1. 

The Fcrmi-LAT collaboration has reported the detection of the famous 
globular cluster 47 Tucanae at a significance of 17cr [H^. The 7-ray spectral 
shape is consistent with a large populations of MSPs. With a typical lumi- 
nosity of ;S 3 X 10^'^ erg s~^ in 7-rays, ~ 20 - 40 MSPs are required to give 
the measured 7-ray luminosity. Using 587 days of data, the number of Fermi 
sources associated with globular clusters has grown to 8 and now fl in 
the 2FGL gD]. 

Both 7-ray detected "garden-variety" pulsars and MSPs have the highest 
values of magnetic field at the light cylinder, B^c^ suggesting that similar 
emission mechanisms operate. The 7-ray luminosity grows with spin-down 
energy oc E at E ^ 10'^^ erg s^^ and Lj oc \fi] at i? ^ lO'^^ erg s^^, with 



tXL ±^ ^ig) "J C^lUa ±Jy 

large scatter in this relation due, at least, to distance uncertainties (Fig. 6 in 
[35]). 



Pulsar wind nebulae 

Given the accuracy of timing parameters for a pulsar, the rotational energy- 
loss rate is known to the uncertainty of the moment of inertia of the neutron 
star. The spin energy not coming out as 7 rays, which usually accounts for no 
more than ~ 10% of the spin-down energy, must come out in a different form, 
e.g., as field energy in the form of a relativistic Poynting wind with such small 
baryon loading. The wind Lorentz factors attains values of ~ 10*. The inter- 
action of the outflowing MHD wind with the ISM makes a termination shock 
where additional particle acceleration can take place. The termination shock 
separates the cold upstream MHD wind with a turbulent flow downstream 
into the ISM. The boundary between the wind and ISM is highly structured 
because of various fluid instabilities. Cold particles can be directly injected 
into the ISM with the wind Lorentz factor, or accelerated at the termination 
shock. 

Identification of pulsar wind nebulae (PWNe) with the Fermi LAT depends 
on finding steady extended GeV radiation. Three PWNe were reported with 
Fermi LAT using 16 months of survey data, namely the Crab nebula, Vela-X, 
and the PWN inside MSH 15-52 [67j. Searching in the off-pulse emission of 
the pulsar yields a candidate PWN related to PSR J1023-5746 and coinci- 
dent with HESS J1023-575. The sources with GeV-detected PWNe have the 
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Fig. 6. COMPTEL 0.7 - 30 MeV data (gray filled squares) and long-term average 
Fermi-LAT fOO MeV - 300 GeV 7-ray data (open circles), and February 2009 (open 
squares) and September 2010 (filled squares) flares of the Crab nebula [71] . 



highest spin-down flux and tend to have the highest spin-down power, with 
Erot ^ 10^^ erg s'l. 

Crab nebula and flares 

The Crab pulsar and nebula is associated with a supernova that took place 
in 1054 CE at a distance of ~ 2 kpc. It is the best example of a center- 
filled plerionic SNR, where the interior is illuminated by synchrotron emission 
powered by the wind of a young, 33 ms neutron star. The spin-down energy 
of the Crab pulsar is 4.6 x 10'^^ erg s~^, of which sa 30% is converted into 
the radiant power of the nebula. The mean spectral index of the nebula in 
the LAT range is -4 between « 100 MeV and 400 MeV, hardening to 
an index of = —1.64 at higher energies [68 (see Fig. [6|. The harder, high 
energy component is the synchrotron self-Compton component (predicted by 
Gould [69]), and well-explained as SSC if the mean magnetic field is « 200 /iG 

[701 Eg. 

Year-scale variations by a several percent are found in Crab nebular fluxes 
measured with RXTE, Swift BAT, and INTEGRAL. The fluxes also declined 
during the first two years of the Fermi mission, consistent with measurements 
of an « 7% decrease between 12 and 500 keV as measured with the GBM 
[7^ . In the period between 2008 August and 2009 April, no variations in the 
100 MeV - 30 GeV flux were found with the LAT [55] • 

Although some evidence for flux variations had been noted during the 
EGRET era [73], Fermi-LAT detected > 100 MeV flaring episodes from the 
Crab lasting for 16 d in February 2009 and 4 d in September 2010 [71] ■ The 
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September 2010 flare was first announced by the AGILE team [74], and they 
also reported an earlier October 2007 flare [7S] . Typical 7-ray flare luminosities 
are « 4x 10'^^ erg/s, and reach an integral photon number flux above 100 MeV 
that is w lOx brighter than the mean flux. 

The PSF of front events of the Fermi-LAT is w 0.12° - 0.15° at > 10 
GeV, so imaging can be achieved to « 0.12°/8 ~ 1' (at 2 kpc, 1' and 1" 
correspond to w 0.6 pc and « 3 x 10^^ cm, respectively). The extent of the 
radio nebula is « 10', and the extent of the Chandra X-ray torus from which 
perpendicular jets emerge is ~ 1.3'. A bright inner X-ray ring in the Chandra 
image is ~ 0.5' in extent. Zooming into the central region with both Chandra 
and HSlj^ reveals a dynamic wind, wisps, and ^ 1" knots that brighten and 
dim. For comparison, one light day corresponds to « 3 x 10^^ cm, far smaller 
than the smallest resolvable knots, yet the 7-ray flares radiate « 1% of the 
spin-down energy. 

The crucial feature of the Crab spectrum may well be that the quiescent 
spectrum displays a strong softening at sa 20 MeV. This value is a factor f« 10 
smaller than the maximum synchrotron frequency at w 200 MeV obtained 
by balancing the timescale for electron synchrotron losses with the Larmor 
timescale on which particles gain energy through Fermi particle acceleration 
mechanisms (see Section 5.5), for a nonrelativistic flow. A mildly relativistic 
flow can enhance the emission by large factors and could be compatible with 
changes in the flow proflle at the termination shock when the ultra-relativistic 
wind slows to mildly relativistic speeds. Considerable theoretical interest has 
focused on perturbations of the flow induced by statistical fluctuations [75] , by 
the generation of large amplitude MHD waves [77] , and by outflowing knots 
formed at an oblique termination shock [TS]. No evidence for the flare energy 
generation is found in the pulsar spin-down behavior [TV . 

Pulsar physics with Fermi 

Three years after launch, over 100 7-ray pulsars are known, and 7-ray emission 
is associated with 8 globular clusters counting 47 Tuc [SB]. In NGC 6624, a 
single MSP dominates the 7-ray energy output [75]. The 7-ray pulsars are 
divided into 31 radio-selected MSPs, 37 blind search 7-ray pulsars, and 38 
young radio-selected 7-ray pulsars. Of the 37 blind search pulsars, extensive 
follow-up observations reveal pulsed radio emission in only 3 cases. Radio 
followup of Fermi pulsars is described, e.g., in [80] . 

Some of the interesting open questions relate to the fleld geometry and ra- 
diation mechanisms for pulsar 7 rays. The significant fraction of Geminga-like 
blind-search pulsars is generally interpreted in terms of a larger 7-ray than 
radio cone. The detection of pulsed emission in the Crab to w 125 GeV [81] 
means that this emission has to be made high in the magnetosphere. Emission 
to these energies requires extreme parameters for a curvature radiation origin. 
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and may imply Compton-scattering in the Crab's magnetosphere. (The cur- 
vature radiation mechanism remains a feasible mechanism for pulsars older 
than the Crab.) The cx i?" behaviour is not well-explained from the lowest 
MSP powers to the highest pulsar powers. Pulsar 7-ray emission can make a 
significant, tens of %, fraction of the Galactic diffuse emission. This fraction 
could be very different in early-type galaxies with only MSPs now active. 
For reviews of Fermi pulsars, see [82] and [83] . 



3 Fermi AGN Catalogs 

The discovery of the 7-ray blazar class is one of EGRET's lasting legacies. Of 
the 66 high-confidence sources associated with AGNs in the 3EG catalog [17] , 
all except one — the nearby radio galaxy Centaurus A — were associated with 
blazars. 

Blazars are sources that exhibit violent optical variability (e.g., chang- 
ing in flux by ~ 50% in one day), high optical polarization (exceeding sev- 
eral per cent), and flat radio spectra with radio spectral index Ur < 0.5 at 
GHz frequencies {Fn oc y~°'^). Superluminal motion at radio frequencies and 
highly luminous and variable 7-ray emissions are also typical blazar properties. 
Blazars themselves are interpreted to be relativistic jet sources powered by a 
supermassive black hole, like radio galaxies, though with the observer looking 
nearly along the jet axis. The variety of multiwavelength SEDs displayed by 
blazars and their misaligned populations can often, though not exclusively, 
be attributed to orientation effects amplifled by the Doppler boosting of the 
jetted radiation (Section 4). 

3.1 LAT Bright AGN Sample (LB AS) and First LAT AGN 
Catalog (ILAC) 

Three lists of AGNs detected with Fermi have now been published by the 
Fermi Collaboration. These are the LAT Bright AGN Sample (LB AS) [38], 
and the First and Second LAT AGN Catalogs, ILAC [3f and 2LAC [51], 
respectively. 

The LBAS is based on 3 months of science observations taking place be- 
tween 2008 August 4 and 2008 October 30, and consists of 106 high Galactic 
latitude (|6| > 10°) sources associated with AGNs. These sources have a test 
statistic TS > 100, corresponding to ^ lOa signiflcance, and are a subset of 
the 205 sources listed in the BSL [36) . 

By comparison, the 3EG \l7\ and the EGR [53] list 31 sources with signif- 
icance > lOcr, of which 10 are at high latitude. Remarkably, 5 of the > 10a 
EGRET sources are not found in the BSL. These are the flaring blazars NRAO 
190, NRAO 530, 1611-1-343, 1406-076 and 1622-297, the latter of which is the 
most luminous blazar detected with EGRET |85) . 
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The ILAC [37] is a subset of the 1451 sources in the IFGL [34] derived 
from analysis of the first 11 months of LAT science data. There are 1043 
IFGL sources at high latitudes, of which 671 are associated with 709 AGNs, 
with the larger number of AGNs than ILAC sources due to multiple asso- 
ciations. Associations are made by comparing the localization contours with 
counterparts in various source catalogs, for example, the fiat-spectrum 8.4 
GHz CRATES (Combined Radio All-Sky Targeted Eight GHz Survey; [86]) 
and the Roma BZCAT blazar catalog [57]. The probability of association is 
calculated by comparing the likelihood of chance associations with catalog 
sources if randomly distributed. 



Table 4. Classes of 7-ray emitting AGNs and galaxies in the ILAC and 2LAC 
"clean" samples 



Class 


Number in 


Characteristics 


Prominent Members 




ILAC (2LAC) 






All 


599 (885) 






BL Lac objects 


275 (395) 


weak emission lines 


AO 0235+164 


. . .LSP 


64 (61) 


< 10" Hz 


BL Lacertac 


...ISP 


44 (81) 


10" Hz < < 10^^ Hz 


30 66A, W Oomao 


...HSP 


114 (160) 


i^Z" >''lO''^ Hz 

pk 


PKS 2155-304, Mrk 501 


FSRQs 


248 (310) 


strong emission lines 


30 279, 30 354.3 


...LSP 


171 (221) 




PKS 1510-089 


. . .ISP 


1 (3) 






...HSP 


1 (0) 






New Classes^ 


26 (24) 






. . . Starburst 


3 (2) 


active star formation 


M82, NOG 253 


. . . MAGN 


7 (8) 


steep radio spectrum AGNs 


M87, Gen A, NGG 6251 


. . .RL-NLSls 


4 (4) 


strong Fell, narrow permitted lines 


PMN J0948+0022 


. . . NLRGs 


4 i-f 


narrow line radio galaxy 


4G + 15.05 


. . . other sources^ 


9 (11) 






Unknown 


50 (156) 







^ Total adds to 27, because the RL-NLSl source PMN J0948+0022 is also classified as FSRQ in 
the ILAG 

^Includes PKS 0336-177, BZU J0645+6024, B3 0920+416, CRATES J1203+6031, GRATES 
J1640+1144, GGRaBS ,11647+4950, B2 1722+40, 3C 407, and 4G +04.77 in ILAG Clean Sample 
^Glass designation deprecated in 2LAG 



Of the 671 associations, 663 are considered "high-confidence" associations 
due to more secure positional coincidences. The "clean" sample is a subset 
of the high- confidence associations consisting of 599 AGNs with no multiple 
associations or other analysis fiags, for example, evidence for extended emis- 
sion. As listed in Table [4| these subdivide into 275 BL Lac objects, 248 flat 
spectrum radio quasars, 26 other AGNs, and 50 AGNs of unknown types. 
The "New Classes" category contains non-blazar AGNs, including starburst 
galaxies and various types of radio galaxies, e.g., narrow line and broad line. 
An AGN is classified as an "unknown" type either because it lacks an optical 
spectrum, or the optical spectrum has insufficient statistics to determine if it 
is a BL Lac objects or a fiat spectrum radio quasar (FSRQ). In comparison 
with the 671 AGNs in the ILAC, EGRET found 66 high-confidence (> 5a) 



Sources of GeV Photons and the Fermi Resuhs 29 

and another 27 lower-confidence detections with significance between 4ct and 
5a, as noted earher. Thus the ILAC aheady represents an order-of-magnitude 
increase in the number of AGNs over EGRET. There are « 300 unassociated 
and therefore unidentified high-latitude Fermi sources in the ILAC. 

3.2 Classification of radio-emitting AGNs and unification 

Different classes of extragalactic AGNs are defined according to observing 
frequency. We have already noted the association of Fermi sources with BL 
Lac objects and FSRQs, which are based on an optical classification. The 
precise definition used by the Fermi team is that an AGN is a BL Lac object 
if the equivalent width of the strongest optical emission line is < 5 A, and the 
optical spectrum shows a Ca II H/K break ratio < 0.4 in order to ensure that 
the radiation is predominantly nonthermal (the Ca II break arises from old 
stars in elliptical galaxies). The wavelength coverage of the spectrum must 
satisfy {Xmax — Xmin)/Xmax > 1-7 in Order that at least one strong emission 
line would have been detected if present. This helps guard against biasing the 
classification for AGNs at different redshifts where the emission lines could 
be redshifted out of the relevant wavelength range. For sources exhibiting BL 
Lac or FSRQ characteristics at different times, the criterion adopted is that 
if the optical spectrum conforms to BL Lac properties at any time, then it is 
classified as a BL object. 

The criterion for classification of radio galaxies according to their radio 
properties stems from the remarkable correlation between radio morphology 
and radio luminosity |88| . The twin-jet morphology of radio galaxies is seen in 
low-power radio galaxies, whereas the lobe and edge-brightened morphology is 
found in high-power radio galaxies, with a dividing line at « 2 x 10^^ W/Hz-sr 
at 178 MHz, or at a radio luminosity of « 2 x lO''^ erg s^^ for the current 
cosmology. Besides a radio-morphology/radio-power classification, radio spec- 
tral hardness can also be used to characterize sources as flat-spectrum and 
steep-spectrum sources. Furthermore, radio galaxies can be subdivided ac- 
cording to the widths of the optical emission lines into broad- and narrow-line 
radio galaxies. Correlations between radio-core dominance and 7-ray lumi- 
nosity supports a scenario where the jet 7-ray flux is greatest along the jet 
direction [g^ . 

Blazars and radio galaxies can also be classifled according to their broad- 
band SED when there is sufficient multiwavelength coverage to reconstruct a 
spectrum from the radio through the optical and X-ray bands and identify a 
uFi, peak frequency i'^^" of the lower energy, nonthermal synchrotron com- 
ponent of the spectrum (see Section 4) . When the peak frequency z^^^" of the 
synchrotron component of the spectrum is < 10^^ Hz, then a source is called 
low synchrotron-peaked (LSP), whereas if the SED has z/^^" > 10"'^^ Hz, then 
it is referred to as high synchrotron-peaked (HSP). Intermediate synchrotron- 
peaked (ISP) objects have 10^^ Hz < i/^^." < lO^^ Hz. SEDs of the bright 
Fermi LBAS sources are constructed in 90J. Essentially all FSRQs are LSP 
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blazars, whereas BL Lac objects have large numbers in all LSP, ISP, and HSP 
subclasses. 

According to the standard unification scenario for radio-loud AGNs ^T\ . 
radio galaxies are misaligned blazars, and FRl and FR2 radio galaxies are the 
parent populations of BL Lac objects and FSRQs, respectively. To establish 
this relationship requires a census of the various classes of sources that takes 
into account the different beaming properties for the Doppler-boosted radia- 
tion of blazars. Even if analysis of data of radio galaxies and blazars supports 
the unification hypothesis, this paradigm still does not explain the reasons for 
the differences between radio-quiet and radio-loud AGNs, or between BL Lac 
objects and FSRQs. 

New classes of extragalactic Fermi sources found in the ILAC include star- 
burst galaxies (Section 6.1), narrow line radio galaxies (NLRGs), radio-loud 
narrow-line Seyfert Is (RL-NLSls), and radio-quiet AGNs. Five NLRGs are 
reported in the ILAC. These objects have narrow emission lines in their opti- 
cal spectrum, suggesting that they are observed at large angles with respect 
to the jet direction, with the surrounding dust torus obscuring the broad line 
region (BLR). 

RL-NLSls have also been recently established as a 7-ray source class [5^ . 
These objects show narrow H/3 lines with FWHM line widths ^ 1500 km s""'^, 
weak forbidden lines {[OIII]/H/3 < 3) and an Fe II bump [SS]- By comparison 
with the ~ IO^Mq black holes in blazars, the host galaxies of RL-NLSls are 
spirals that have nuclear black holes with relatively small (~ 10® - 10® Mq) 
mass that accrete at a high Eddington ratio. The detection of these objects 
challenges scenarios (e.g., [Mj) where radio-loud AGNs are hosted by elliptical 
galaxies that form as a consequence of galaxy mergers. 

The ILAC includes 10 associations with radio-quiet AGNs. In 8 of these 
cases, at least one blazar, radio galaxy, or CRATES source is also found close 
to the 7-ray source. In the remaining two cases, the association probabilities 
are weak. Thus none appear in the ILAC "clean" sample. In some of these 
candidate radio-quiet 7-ray sources, such as the Sy 2 galaxies NGC 4945 or 
NGC 1068 [SS], which are also starburst galaxies, the 7 rays could be made by 
cosmic-ray processes rather than from a radio-quiet Sy nucleus. More extensive 
searches for 7 rays from Swift-BAT AGNs has not established, however, that 
radio-quiet AGNS are GeV 7-ray emitters [SS] . 

3.3 Properties of Fermi AGNs 

Various correlations are found by comparing 7-ray properties of Fermi AGNs 
according to their radio, optical, or SED classification. Probably the most 
pronounced correlation is between the > 100 MeV 7-ray spectral index 
and optical AGN type. FSRQs have significantly softer spectra than BL Lac 
objects, with (r^) =S 2.40±0.17 for FSRQs and (r^) ^ 1.99±0.22 for BL Lac 
objects in the LB AS [SS]. The SED classification shows that the mean 7-ray 
spectral index (F^) = 2.48,2.28,2.13, and 1.96 when the class varies from 
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Fig. 7. Gamma-ray spectral slope F-y of BL Lac objects (open blue circles) , FSRQs 
(open black squares), FRI radio galaxies (red circles), FR2 radio sources (green 
squares), and star- forming galaxies (magenta diamonds), are plotted as a function 
of their 100 MeV - 5 GeV 7-ray luminosity L^. 



FSRQs to LSP-BL Lacs, ISP-BL Lacs, and HSP-BL Lacs, respectively. The 
progressive hardening from FSRQs to BL Lac objects can be seen in Figure [7j 
which also compares with values for radio galaxies and star-forming galaxies 



Fermi data reveal complex 7-ray blazar spectra. FSRQs and LSP-BL Lac 
objects, and most ISP blazars with sufficiently good statistics, show breaks 
in the « 1 - 10 GeV range [98]. This was already apparent from the first 
observations of the bright blazar 3C 454.3 [22] i to be discussed in more detail 
below. The HSP blazars, though, are generally well-described by a flat or 
rising vF^ SED in the GeV range, with vF^, peak frequencies between « 100 
GeV - TeV energies implied by VHE data. 

Only 121 out of 291 BL Lac objects had measured redshifts at the time of 
publication of the ILAC. For sources with measured redshift, BL Lac objects 
are mostly found at low redshift, z ^ 0.4, with only a few HSP BL Lac objects 
at higher redshifts. By contrast, the FSRQs span a wide range from z « 0.2 
to the highest redshift ILAC blazar with z = 3.10. 
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This significant redshift incompleteness hampers interpretation of AGN 
properties, in particular, Fig. [7j which can only display sources with known 
redshifts. For these sources, the hard-spectrum BL Lac objects typically have 
much lower than the FSRQs. This divide has been interpreted as a change 
in the accretion regime at « 1% of the Eddington luminosity [55] ■ In addition, 
the nearby radio galaxies with z ^ Q.l inhabit a separate portion of the 
vs. plane, and are characterized by lower 7-ray luminosities than their 
parent populations. The two more distant steep spectrum radio sources, 3C 
207 (z = 0.681) and 3C 380 (z = 0.692), and the FR2 radio galaxy PKS 0943- 
76 {z — 0.27) fall, however, within the range of 7-ray luminosities measured 
from FSRQs. Indeed, steep spectrum radio sources are thought to be slightly 
misaligned FSRQs. 

3C 454.3 and FSRQs 

The FSRQ 3C 454.3, at z = 0.859, underwent giant flares and became the 
brightest 7-ray source in the sky for a week in 2009 December and 2010 April 
|100j . and again in 2010 November [3D]. The outburst in 2010 April triggered 
a pointed-mode observation by Fermi. During the December outburst, its 
daily flux reached -F_8 — 2200(±100), corresponding to an apparent isotropic 
luminosity of L^^o « 3 x 10^^ erg s^^, making it the most luminous blazar yet 
detected with Fermi. In its 2010 November outburst, it reached F„8 = 6000 
and Lisa ~ 10^" erg s~^, becoming k, 5 times brighter than the Vela pulsar. 

Figure [s] shows the light curve of 3C 454.3 [30], which can also be found 
at the public website of Fermi monitored sources]^ The fluxes are plotted in 
durations of one day over the course of the Fermi mission. The flux rises to a 
plateau level preceding 7-ray outbursts, with the 2008 July outburst showing 
strong resemblance to those in 2009 August and December, and 2010 April and 
November. Intense flaring occurs during periods of enhanced activity, which is 
to say that the flares are not isolated events, but seem to occur during periods 
of enhanced accretion activity. 

As noted already in the initial Fermi report [55]) the spectrum of 3C 454.3 
breaks strongly by 1.2(±0.3) units at E\jr « 2 GeV. Such a break is incon- 
sistent with simple radiative cooling scenarios, which predict a break by 0.5 
units. The more recent analyses of data from the major outbursts of 3C 454.3 
[1001 [30] confirm the strong spectral break and finds that E^r is very weakly 
dependent on the fiux state, even when the fiux changes by more than an or- 
der of magnitude (see Fig. |9]) . No consistent pattern expected in acceleration 
and cooling scenarios 1102) is found in the spectral index/flux plane. 

The origin of the spectral break in 3C 454.3 bears on several important 
issues in FSRQs: the location of the 7-ray emission site; the source of soft 
target photons in Compton-scattering models; and the relation of FSRQs and 
BL Lac objects in view of the disappearance of such breaks in ISP and HSP 
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Fig. 8. Fermi LAT daily light curve of 3C 454.3 |30j, showing the giant flares 
in December 2009, April 2010, and November 2010 (MJD 55200 corresponds to 4 
January 2010). Inset shows the light curve on a logarithmic scale, with black points 
from EGRET [ITj, and blue points from AGILE [ToT] . 



blazars. Such a break would be readily understood if the target field was suf- 
ficiently strong to attenuate the blazar radiation by 77 absorption processes, 
but the intense Ly a radiation field at 10.2 eV observed with GALEX [103] 
would make a break at Ebr ^ 30 GeV |104j . Models employing photon at- 
tenuation deep within the BLR by He II recombination and Ly a photons 
with E > 54.4 eV [105] have been proposed to explain this break, as have 
nonthermal electron scattering models with accretion-disk and BLR photon 
targets |106j . The spectral break could also be due to Klein-Nishina effects in 
scattering, as has been proposed to explain the SED of PKS 1510-089 |107j . 
The KN break due to upscattered Ly a radiation occurs at a few GeV, and 
the observed break energy is insensitive to the Doppler factor. This scattering 
problem is treated in Section 4. 

PKS 2155-304, Mrk 501, and BL Lac objects 

PKS 2155-304, an X-ray selected BL Lac object at ^ = 0.116 and an EGRET 
source |108j . is one of the most prominent representatives of the HSP blazar 
population. During a period of extraordinary flaring on 2006 July 28, PKS 
2155-304 exhibited a succession of 7-ray flares varying on time scales as short 
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Fig. 9. Fermi LAT vF,^ SEDs for four different time periods corresponding to pre- 
flare (magenta; lower spectrum), plateau (green; third from top), flare (blue; top 
spectrum), and post-flare (red; second from top) periods, along with the fitted func- 
tional forms, including a broken power law (solid), log-parabola (dashed), and power- 
law with exponential cutoff (dashed-triple-dotted) functions [3D] . The inset displays 
inbreak (red) and Scutoff energies (blue) as a function of flux for the different periods 
(filled symbols) and for historical flares (open symbols). 



as « 5 min with apparent luminosities > 10 erg/s — larger still when EBL 
absorption is included [109) . This is an extreme, hypervariable source in the 
sense of equation ( 14 1 . 



An 11 day campaign between 2008 August 25 and 2008 September 6 was 
organized early in the Fermi mission to measure the SED at optical (ATOM; 
Automatic Telescope for Optical Monitoring), X-ray (RXTE and Swift), and 
7-ray (Fermi and HESS) frequencies. It turned out to be in a low state, well 
fit by a standard one-zone synchrotron/SSC model, with Doppler factor Sd — 
32, magnetic field B' = 0.018 G, and comoving size R' = 1.5 x 10^^ cm 
(corresponding to a variability time of 2 d; |110[ )p^ When flaring, however, 
the one-zone synchrotron SSC model for PKS 2155-304 completely fails, or at 
least requires 6u ~ 100 [113[ . 

One-zone synchrotron/SSC models with Su ~ 10 give good fits to the long- 
term average spectra of other HSP BL Lac objects such as Mrk 421 and Mrk 



501. The spectacular multiwavelength SED of Mrk 421 shown in Figure 10 rep- 



The reader is assumed to be familiar with the synchrotron/SSC model; see, e.g. 

[iiTi[iia[Tia[3]. 
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Fig. 10. SED for Mrk 501 averaged over aU observations taken during the mul- 
tifrequency campaign performed between 2009 March 15 (MJD 54905) and 2009 
August 1 (MJD 55044) |114| . The legend reports the correspondence between the 
instruments and the measured fluxes. 



resents one of the most detailed multiwavclcngth blazar SEDs yet assembled 
|114j . The optical and X-ray data have been corrected for Galactic extinction, 
but the host galaxy, which is clearly visible at IR/optical frequencies, has 
not been subtracted. The MAGIC and VERITAS data have been corrected 
for EBL absorption (Section 7.3). A single-zone synchrotron/SSC model with 
12<Sb^ 22, B' ^ 15 - 30 mG, and i?' w 8 - 50 It-day, gives a good fit to 
the nonflaring SED shown in Figure \T0\ By comparison, radio galaxies are fit 
by synchrotron/SSC models with Doppler factors of order unity. We return 
to this point in Section 7. 



3.4 Second LAT AGN Catalog (2LAC) 

The 2LAC ^84i is based on the first two years of scientific data taken with the 
Fermi Gamma ray Space Telescope. It contains 1016 |6| > 10° sources that are 
associated at high confidence with AGNs (1319 of the 1873 IFGL sources are 
at |6| > 10°). The 2LAC clean sample comprises 885 sources, consisting of 395 
BL Lac objects, 310 FSRQs, 156 blazars of unknown type, 8 misaligned AGNs, 
4 RL-NLSyl galaxies, 10 AGNs of other types, and 2 starburst galaxies (NGC 
4945 has fallen out of the list); see Table |4j The photon index distribution 
of the blazars of unknown type suggests that they comprise roughly equal 
numbers of BL Lacs and FSRQs. Of the 395 BL Lac objects, 220 (- 55%) 
lack measured redshifts, and this fraction is roughly the same for LSP, ISP, 
and HSP BL Lac objects. 

As shown in the 2LAC, threshold sensitivity in terms of photon flux is 
strongly dependent on source spectral index, whereas energy flux is not, and 
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detection of FSRQs and BL Lac objects is complete to an energy flux of 
« 5 X 10~^^ erg/cm^-s. Because more and deeper surveys have taken place in 
the northern hemisphere, the smaller fraction of southern-hemisphere BL Lac 
objects indicates that some 60 unassociated sources at negative declination 
are expected to be BL Lac objects. 

Besides validating the results of the ILAC with additional data and im- 
proved analysis techniques, the 2LAC shows that the mean fractional vari- 
ability of FSRQs is greater than that of BL Lac objects. The duty cycle, as 
defined by the fraction of time when the average flux exceeds 1.5 standard 
deviations above the mean flux, is « 5 - 10% for bright blazars. The num- 
ber of BL Lacs and FSRQ have increased by 44% and 35% from the ILAC 
to 2LAC clean samples, respectively, with the future increases in the num- 
ber of FSRQs expected to be even more modest. This is due to cosmological 
effects, as reflected in the flattening in the log A^-log S distribution, and K- 
corrections that pushes the GeV peak of the vF^, spectra to lower energies 
where the Fernii-LAT effective area rapidly declines and sensitivity to FSRQs 
gets worse. 

Blazars make up 97% of the extragalactic 7-ray sources detected with 
the Fermi-LAT. The number of misaligned sources has stalled, with the same 
number of radio galaxies — eleven — in the 2LAC as in the misaligned AGN 
paper. But these are not the same objects. The radio galaxies 3C 78, 3C 111, 
and 3C 120 are not found in the 2LAC, evidently because a variable jetted 
emission component that contributed to the 7-ray emission in the past has 
gone quiet. Now found in the 2LAC are the FRI radio galaxies Fornax A 
and Gen B, and the head-tail radio galaxies IG 310 |115j . the latter of which 
is also a MAGIG source [116) . The radio/7-ray connection [117] and GeV- 
TeV synergy |118j resulting from the Fermi-LAT for AGN studies, let alone 
Galactic sources, cannot be adequately reviewed here; see [84] and [40] . 

4 Relativistic Jet Physics 

Here we develop and apply relativistic jet radiation physics 3J to some puzzles 
arising from the Fermi observations. 

4.1 GeV spectral break in LSP blazars 

As already described in Section 3.3, an important new result in blazar physics 
not anticipated before the launch of Fermi is the prevalence of a spectral soft- 
ening in the 7-ray spectra of low synchrotron-peaked (LSP) blazars, including 
both FSRQs and LSP BL Lac objects. The spectral softening occurs generally 
between 1 and 10 GeV [Ml 133 • For 3G 454.3 [mdnD], the break is > 1 unit, 
and the energy of the break is rather insensitive to flux (Fig. |9]) . A softer spec- 
tral break, possibly consistent with radiative cooling, is found at multi-GeV 
energies in the Fermi/MAGIG spectrum of FSRQ 4C +21.35 [120]. If a GeV 
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spectral softening and significant GeV flux variability accompanied by modest 
spectral variability is typical of blazars, then some physical process should be 
able to explain this behaviour. 

7 rays from external Compton scattering 

We treat the problem that blazar 7 radiation arises from the scattering of a 
quasi-isotropic target photon field [121 in the Klein-Nishina regime. Compton 
scattering takes place in the Thomson limit when A-fe' ^ 1 (eq. 6.123 in .3J), 
where 7 and e' are the electron Lorentz factor and photon energy in the 
comoving frame. 

Consider an external isotropic, monochromatic photon field with pho- 
ton energy e^, (in mgC^ units). The average energy in the comoving frame 
is (e') = Fci,. The upscattered photon energy in the Thomson regime is 
e'g, ^ (4/3)72e'. Hence e'e'^ = {4:/3)-/'^e''^ = (47e')^/12. Thus scattering is in 
the Thomson regime provided 12e'eJ^ ^ 1, or e* £ {Sn/ r)/[12ec{l + z)], im- 
plying Ec{GeV) ;S 12/ [^^^(eV)] for (5d = F. Away from the endpoints, there- 
fore, the scattered photon energy marking the transition to the KN regime is 
very weakly dependent on Doppler factor in highly beamed relativistic jets. 

The spectrum of an isotropic nonthermal power-law electron distribution, 
when transformed to the stationary frame, remains a power-law with the same 
index as in the comoving frame, but with a strong angular dependence. The 
endpoints of the distribution are the low-and high-energy electron Lorentz 
cutoffs, 7^ and 72 respectively, boosted by 5u [122\ . Consequently, the shape 
of the scattered spectrum cannot depend on provided that 71 = (5d7i ^ 
7KN ^ (4e*)"^ < S-DI2 = 72- 

If the break energy observed in 3C 454.3 at « 2 GeV is due to the transition 
to scattering in the KN regime, then the underlying target photon energy 

« 6 eV. This is close to the energy of the Ly a photon at 10.2 eV, or H/3 
at 2.55 eV. If the origin of the spectral break is due to scattering of nearly 
mono-energetic line radiation, then this would (1) place the scattering site 
within the BLR, and (2) explain the near constancy of the spectral break, 
independent of flux state. 

We use the method of Georganopoulos et al. (2001) |122| . applied to 
isotropic radiation fields in the stationary frame of the supermasssive black 
hole and BLR. The differential Compton-scattered spectrum when isotropic, 
monoenergetic electrons Compton upscatter isotropic, monochromatic target 
photons, is given for the full Compton cross section in the head-on approxi- 
mation by 



•c 



3 caTUp /eA Nej-f, fis)d^ 



4 dl UJ 



7 



(28) 



Ll + (l/^e) 



(eq. 6.129; [3]), with 



N'ij') 



7' = 7/fc 



(29) 



38 Charles D. Dermer 



(eq. 6.124), 

= + z)e = e, , Q,^Q^ , (30) 

Fc(g,re) = 2glng + (1 + 2q){l - q) + l^i^^^(l - q) , (31) 

^"^^ , r, = 476, , (32) 



and q is restricted to the range (47^)-! < q < 1 1121 (eqs- 6.74-6.76, 
6.125-6.127; citedm09). Restriction to the Thomson regime occurs when <C 
1 or 476* <^ 1. The final term in eq. pl[ ) dominates for scattering in the KN 
regime. (A simpler form for analytic work is the isotropic Thomson cross 
section, 

FtM^) - ^ (1 - e) (33) 

(eq. 6.71 [?), with e — es/47'^e* (eq. 6.70), which assumes isotropic scattering 
in the electron rest frame.) 

Thus an isotropically distributed nonthermal electron distribution in the 
comoving frame gives a Compton-scattered spectrum 



- ' « /" ^ /" "7 ™ ^c(,. r.) , ,34) 

4 'iTTUi^ Jo ^* Jjmin 7 

(eq. 6.123; considering only upscattering), with 



7™» = |(l + yi + ^) (35) 



(eq. 6.90). Substituting eq. (29) for an isotropic power-law electron spectrum 
in the comoving frame, so that N'{j') — N^oj' ^pH{j; 71, 72) [^scattering an 
external monochromatic line spectrum u*(e*) = u*(5(e* — e*o) gives 

f^''" = ! f?^ (-)' r d7 7"^^+^) Fc{q, r,) . (36) 

This is numerically solved to give the results shown in Fig. |11[ which were 
compared in 110.0! with the rapidly falling spectrum of 3C 454.3. 



Fig. 12 shows a model where the 7-ray continuum observed from 3C 
454.3 with the Fermi LAT [lOOj are formed by power-law nonthermal jet 
electrons, with number index p = 3.17, that Compton-scatter Ly a line pho- 
tons. The model is insensitive to values of lower and upper comoving electron 
Lorentz factors ^min and jrnax provided jmin ~ 10^ and jmax ~ 10''. The 
Klein-Nishina softening from a power-law electron distribution gives a poor 
fit to the data, but if one assumes that the underlying electron distribution is 
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Fig. 11. Spectrum of Compton-scattered Hne radiation. Left: Dependence on 
Doppler factor Su for sources at different redshifts. Right: Dependence on endpoints. 
Within the endpoints, the spectral shape of the function is independent of So. 
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Fig. 12. Model for the 7-ray spectrum of 3C 454.3 when isotropic power-law elec- 
trons with index p — 3.17 in a relativistic jet Compton scatter Ly a photons [lOUj . 
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curved, as in the case of a log parabola function 125 , for example, it may be 
possible to obtain a good fit to the 3C 454.3 data for this model. 
For integrations over blackbody spectra, substitute 

u^,bb{£*) = — r| 7 — rpTs r (37) 



into eq. (34). The dimensionless temperature of the radiation field is 6* = 
k^T / rrif.c''' . The CMBR can be approximated as an isotropic, monochromatic 
radiation field u*(e*) — wo5(e* — eo), with uq = 4 x 10^^'^(1 + z)* ergs/cm"^ 
and eo 2.70Ocmbb. = 1-24 x 10^^(1 + z). The condition q < 1 implies 



Compton emissivity and electron energy-loss rate 

In the calculation of time-dependent blazar spectra with evolving electron 
distributions, or in pair cascade calculations, it is necessary to consider the 
full Compton cross section in scattering and electron-energy loss calculations. 
When 7^1, the head-on approximation, where the incident photon is as- 
sumed to be directed opposite ("head-on" to) the direction of the electron 
when transformed to the electron rest frame, can be employed. 

For general treatments, the time-dependent emissivity esjci^s, ^s'^fTt) is 
calculated. For external isotropic photon spectra, the angle-dependent emis- 
sivity when jetted nonthermal, isostropic electrons with electron spectrum 
ne(7, fie) Compton-scatter photons of an isotropic monochromatic photon 
field is given, in the head-on approximation, by the expression 

esjcies,n^) = ^,carUoC-Y H "^^^^^Fdq, re)Hiq; -^A) (38) 
4 e* Ji 7 

(eq. 6.74). Here the electron spectrum is written in the stationary frame where 
the radiation field is isotropic. The vLi, Compton luminosity for a single elec- 
tron is given by 

^.ic(..) = 5^(^)Vc(,,r.) Hie,, -1^) . (39) 

4 7 1 + e*/7 1 + ^ 

The electron energy loss rate by Compton scattering is given by 
-7c = de, Lcies) = ^^^^^^ j^^ de. Fciq^r^) 



The Heaviside functions H{x~a) and H{x; a, h) are defined such that H{x~a) — 1 
if a:: > a and H{x — a) = otherwise, and H{x; a,b) — 1 if a < x < b and 
H{x; a,b) = otherwise. 
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' n^r ^3 Fciq,r,)- (40) 



12ccrTWo7^ 



Using the isotropic Thomson kernel, equation ( |33[ ) becomes in the Thomson 
regime, /"e < 1, 

caTUf) ,es.2r 3 a;^ 1 r x -, 7 



(41) 

Here a; = £5/7 so that q = x/[re{l — x)]. To lowest order, 

recovering the familiar result in the limit 7 ^ 1, e*/7 ^ 1. Expanding the 
energy-loss rate expression, again using the isotropic Thomson kernel, gives 
the expansion 

- 7C.=1 ^ IcaTUoi' (1 - ^r, + + 0[r!]) (43) 

correct to 0{r^). 



4.2 Leptonic jet models 

A photon flux of F^g — 2200 from 3C 454.3 aX z — 0.859 implies an apparent 
> 100 MeV isotropic 7-ray luminosity of = 3.3x10^^ ergs"^ [TUn]. The flux 
reached a factor k, 6 larger in the November 2010 outburst [30], corresponding 
to an apparent luminosity of ~ 2 x 10^° erg/s, which is the record-holder for 
all blazars, including PKS 1622-297 from the EGRET era [85) . 

Using the measured flux and a one-day variability timescale at the time 
that the most energetic photon with E 2Q GeV was detected implies a 
minimum Doppler factor of <5D.min ~ 13, as we show. Assuming that the 
outflow Lorentz factor F w 20, consistent with the inferred value of (5D,min 
and with radio observations at a different epoch |126j , then simple arguments 
suggests a location R £ cr^tvar/(l -I- z) « 0.2(r/20)^(tvar/day) pc, which is 
at the outer boundary of the BLR. Flux variability on time scales as short as 
3 hr was measured at another bright flux state [30], which suggests that the 
7-ray emission site would be even deeper in the BLR. This stands in contrast 
to inferences based on coherent optical polarization changes over timescales 
of weeks prior to a 7-ray flare in 3G 279 |127) , which places the emission site 
at much larger distances. The MAGIG detection of VHE emission from PKS 
1222-1-21 (4G -1-21.35) |128| is even stronger evidence for 7-ray production at 
the pc scale or farther from the central engine. 
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Jet Doppler factor 

Emission from bulk magnetized plasma in relativistic motion provides the best 
explanation for the large apparent isotropic luminosities Liso, energy releases 
Siso, and short variability times tvar from cosmological 7-ray sources. The 
comoving emission-region size scale i?' is related to ivar through the expression 
R' ^ cFtvui/il + z). Combined Fermi LAT and GBM observations give the 
best GRB data from which to determine the minimum bulk outflow Lorentz 
factor r^iri through 77 opacity arguments. For blazars, multiwavelength X- 
ray measurements combined with Fermi-LAT data are important for inferring 
the minimum Doppler factor. 

It is simple to derive -Tmm in a blast-wave formulation, noting that the 
internal photon energy density 



var 



using R fa /^^civar/(l + z). The optical depth for 77 — > e+e processes is 
r^^(ei) = i?'crT(ei/2)u;(2/ei)/(mec2), where i?' = R/T and e' = 2/e[ from 
the threshold condition. The condition Tjy{e'i) < 1 with the relation re'i/{l + 
z) = ei implies 



aTrfi(l + ^)V.-ei 



1/6 



^ (45) 



(l + ^)2ei 



|129j . Here is the vFi, flux at photon energy rrieC^e, which is evaluated at 
e = e due to the peaking of the 77 cross section near threshold. 

The minimum Doppler factor 6min {So ~ [^(1 ~ P cos6)]~^) defined by 
the condition Tjj(ei) = 1, can be estimated to w 10% accuracy compared to 
results of more detailed numerical calculations through the expression 



aTdl{l + z)y,e^^'^' 



(46) 



where is the vFi, spectrum measured at frequency v — m^c^e/h, ty is the 
variability time, and Ei = meC^ei is the maximum photon energy during the 



period in which is measured. The vF^ flux in eq. (46) is evaluated at 
e = e = 2(5q/(1 + z)'^ei from the pair-production threshold condition. 



For 3C 454.3, writing = 10 "^°/-io erg cm ^ s ^ in eq. (46 ) gives 5min ~ 
16.5 [/_io£^i(10 GeV)/t„(6 hr)]^/^ and£;( keV) ^ 3.4(,5^,„/15) Vi;i(10 GeV). 
Contemporaneous Swift XRT observations give the i/F^, flux between 0.2 and 
10 keV. From Swift pubhc data [TOS], vF^iA keV) w 6 x 10"" erg/cm^-s, so 
that Smin ~ 15. Because 5d > Smin, we write S — 20(52o and S20 ~ 1- This value 
compares favorably with — 24.6 ± 4.5, bulk Lorentz factor F = 15.6±2.2, 
and observing angle 9{°) — 0.8 ±0.2 obtained from superluminal observations 
[126j . derived from measurements made at a different time. 
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Variability time scale 

In the case of 3C 454.3, the emission region for one of its major flares is con- 
strained in a colliding shell geometry to be at distance R ^ 2cr^tyl{\ -\- z) ~ 
O.ir'|ot„(6 hr) pc, assuming Tao = r/20 = 1. This is within the BLR [T^. 
Ways to avoid this conclusion are to assume that F ^ 20, which would be 
associated with large jet powers, or for only a small portion of the emis- 
sion region to be active, which would lower the jet's radiative efficiency given 
that only a small fraction of the emitting surface is radiating. This might 
be implausible in view of the already large apparent isotropic 7-ray lumi- 
nosity of 3C 454.3. RecoUimation shocks (e.g., |131| ) at the pc scale might 
reduce the characteristic size of the emission region, though the contrast 
(c(5t„/[(l -I- z)_Rj,c(pc)] '-^ 2 X IQ^^ty-pc/ Rpc) between the size of the emis- 
sion region and location R = Rpc pc seems unexpectedly small even for a 
narrow jet. Alternately, flaring episodes with short variability times might 
take place within the BLR, whereas the more slowly varying emissions could 
be radiated by jet plasma at larger distances, but then we would expect large 
spectral variations due to the different target photon sources. 

The uses of variability studies to infer properties of blazars, GRBs, or 
other sources is an ongoing concern, because flow properties are inferred from 
the variability timescale. The auto-correlation function can be used to infer a 
characteristic variability timescale, as can breaks in the power density spec- 
trum of sources [132] . Other techniques for variability analysis of the extensive 
compilation of blazar light curves are currently in development. 



Equipartition field and jet pow^er 



Here we perform an equipartition power analysis for the giant outbursts from 
3C 454.3 [1001 130j . In the blob framework, where a spherical emitting region 
with radius R' in the comoving frame moves with Lorentz factor F at an angle 
9 to the line of sight, the absolute jet power is [1331 1134j (Section 5.5, 3J), 
including the photon power. 



27rr'2/3c 



F' 



Stt 



3^3/2 



(47) 



Here — 10^^^_9 erg cm^^ s^^ is the measured bolomctric photon energy 
flux (47rd|# = 3.6 x lO'*^ erg s"^ for = 1; = 10 and 50 during the 
December 2009 and November 2010 outbursts, respectively), and the factor 
F^ /26\) in the last term of this expression recovers the increased energy re- 
quirements due to debeaming when 6 ^ 1/F oi 5d ^ r, in the limit that the 
jet opening angle is ^ 1. The factor x is the deviation from equipartition, so 
that S^B — xiSoBeq) and 



B, 



c2(l + Cpe)ln(e2/ei) 



4^(1 + z)epkUlcaTV^ 



2/7 



(48) 
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(eq. 7.80, 3^). Here /c^^ is the i/Fi, peak of the synchrotron spectrum, Ucr — 
B'^^/Stt, Bcr — 4.414 X 10^'^ G, and is the relative energy in protons to 
nonthermal electrons which emit synchrotron radiation with a spectrum a — 
0.5 between observed photon energies ei and €2- 

Letting {1 + (pe) \n{€2 / ei) — 10^ A2, the peak synchrotron frequency i^pk = 
10"'^^j^i3 Hz, and /c 4. = 10~^"/-io erg cm~^ s""'^, substitution of parameters 
for 3C 454.3 gives " 

(/-10^2)2/' 



^20 ^13 (d) 

i.e., an equipartition field of a few G. The corresponding jet power is 



"d i^i^ "Ja ' 

+ 5 X 10-^5^] . 

The Eddington luminosity for a lO^Mg Solar mass black hole is 1.26 xlO'''' 
erg s~^. Estimates for the black-hole mass in 3C 454.3 range from Mg w 0.5 
[103] to Mg « 4 [135]. Thus we see that the jet power from 3C 454.3 would 
be super-Eddington if x departs from its equipartition value by more than a 
factor sa 4. Moreover, the system cannot be severely debeamed, though this 
is already unlikely from other arguments, e.g., core dominance parameter and 
superluminal motion observations |126j . 

By taking the ratio of the synchrotron and Compton vFi, peaks, using the 
relations epk,syn = {S/2)j^^^{B/B^r)S/{l + z) and epk,c = {'i/5)jl^rej/{l + 
z), we have = {i^pk,c / repk,syn){B / B„) or 

i..(eV)^4.6^^--^(^""^-^)^(^^^ (51) 

1^13-^20 

This suggests that the soft photon energy scattered to make the GeV emission 
is a few eV, which would correspond to either atomic line radiation or scattered 
Shakura-Sunyaev accretion-disk radiation. 

The energy density of the external radiation field can also be estimated 
from the 3C 454.3 observations. Let Ac represent the ratio of the energy 
fluxes in the Compton and synchrotron components. If u^ denotes the energy 
density of the target radiation field in the stationary frame, and ub denotes 
the comoving magnetic-field energy density, then for 7^ ^ 1, 4m*/^^/3 « 
AcU'b, or u*(crg cm^^) ^ 0.007(v4c/10)B2(3G)/r|o This can be compared 
to energy densities in the BLR using relations between the BLR radius and 
line luminosity [103] . The implied energy densities are an order of magnitude 
larger, suggesting that either i? > 3 G or < 20. 
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4.3 Hadronic jet models 

Measurements of air showers induced by cosmic rays impacting the atmo- 
sphere give the best evidence for the existence of powerful accelerators of 
UHECRs with E ^ 10^* eV, with Larmor radii so large that they almost 
certainly originate from beyond the Galaxy. 



Adiabatic losses, and photopair and photopion losses on the 
CMBR 



The radial scale factor 

B^j^ (52) 

at redshift z. In an adiabatic expansion process, dQ = = dU +pdV and the 
energy content U — uV , where u is the energy density and V is the volume. 
Thus udV + pdV = -Vdu, or {u+p)dV = -Vdu, so 

dV dx dy dz dR du 

— = 1 1 = 3— = . (53) 

V X y z K u+ p 

For a relativistic gas with p = u/'i, 

3dR du „ „ 3c?lnu . 
^ = -3dlni?= , (54) 



R if) 

which implies u cx Because the energy density U = uV cx R^^ , 7 oc R~^ , 
so -7 cx -{l/R^)dR/dU (see Section 9.4 in [3]). 
With R given by equation ( 52 ) , 



-j^{l + z)-^\ — \J^Ho^/n^lT7)^Ti2^, (55) 

for a flat cold dark matter cosmology with cosmological constant A. Thus 

din 7 1 dz 



dt (1 + z) ' dU 



= Ho , (56) 



The mean-free path for energy losses due to photopair production is given 
in Section 9.3.2 of |3| by the expression 

cE c 



\dE/dt\ \d\n^ldt\ 
+ Gpc., , npc, (57) 



{l + zfF^,[E^8{l + z)] " [l + zfT^M 

where 
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= 0.74+ 1.78 ln(y/2) + ^l^lP , (58) 

taking the asymptote j/ ^ 1 or Ei^ ^ (1 + z)^^. If the ankle in the spectrum 
is a consequence of photopair losses of ultra-high energy protons, as advocated 
by Berezinskii and colleagues |136j , then the energy of the ankle would reflect 
the redshift epoch of greatest source production of the UHECRs, noting that 
r^e oc (l + z)~^ at the maximum energy loss of protons due to photopair losses 
on the CMB. 

Approximately, then 

r rflnJ^ n (l + z)3[0.74 + 1.781n(y/2)] 

L dt J-^^ yy(Gpc/c) ' ^ ' 

when y = Eifi{l + z) > 1. The conversion (Gpc/c) = 1.03 x lO^'^s = 3.2 x 10^ 
yr. 

The mean-free-path for a proton to lose energy through photopion losses 
on the CMBR is given by Ref. [3], Section 9.2.4, by 

N c 13.7exp[4.0/£;2o(l + z)l 

\d\\iE/dt\^-„ (1 + zy^yl^- A.\)/E2o[l + z)\ 

which, though derived in the limit i?2o ^ 4/(H-z), gives a good approximation 
even at higher energies. 



Photopion production efficiency 

From the relations between time elements in the stationary frame (starred), 
the comoving frame (primed), and the observer frame (unscripted), R — 
f3ct^ = pert', so f = RIpPc for relativistic flows, and t = {I + z)R/pr^c, 
since t = {1 + z)t' /F. The comoving dynamical time scale is therefore given 

by trf,„ = R'/c = R/Ec. 

The energy flux 

^ _ d£ _ L^. _ d£^ _ 

" dldt ~ 4^ ' * " " ' ^ ' 

noting F£' = and dt^ = dt' /E. Therefore ei(e) = 47rd|/e = E'^e'L{e'), and 

= , D2 2 ' so n' e' = , ^ , . 62 

The causality constraint for rapidly variable emissions from a relativistic blast 
wave is 

, , Ar ^ At' At,. EAt , , 

At' = —<—^-^ = , 63 

I C 1 C (1 + Z)C 

and t^ar = (1 + z)R/E^c. 



Sources of GeV Photons and the Fermi Resuhs 



47 



We write the target vF^, photon spectrum as 

h = A,. [{—TH{e,^ -e) + { — tH{e - e,,)] (64) 

The threshold for photopion {(j>TT) processes is 7^6' ^ mTr/nie = 400, using the 
approximation of |137| that the product of the photopion inelasticity and cross 
section is Kp^a^^^ = 70 above comoving photon energy e' = 400. Because 
7' = Ep/ruipC^ and e' = (1 + z)e/r, the threshold escaping proton energy is 



j^thr ^ iOOmpC^r^ j^pk ^ mmpc^r^ ^ 3.7 x 

^ il + z)e P il + z)epk {l + z)epk 

where 



1000 

The timescale for energy loss due to photohadronic processes is 

de'n'(e'). (66) 



Thus 



tpy KK) = ..^^9 / ■ (67) 



Using eq. ( 64 ) gives 



^ P '~ R^m,c^repkil-b)il + z) 



The photopion production efficiency 



At Ep^Ef, 



Thus 



„pfc - n fRPfc^ - iKp'yCrp^)dlfe^^ 



i^V-' , Ep < Ef 



(68) 



Wi^.)^^. (69) 



VpyiEp) - <! f (1^)1- , a < 1 . , (71) 

, Ef<Ep 

1 , a > 1 
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Table 5. Band function fits to GRB 090926 

Time interval f-e'^ a b ^pk 

(a) : 0.0 - 3.3 s 3.5 1.58 ± 0.03 -0.64^0 gl 0.66 

(b) : 3.3 - 9.7 s 6 1.45 -0.46 ' 0.56 

(c) : 9.7 - 10.5 s 5 1.41 -1.69 0.41 

(d) : 10.5 - 21.6 s 1 0.30 -0.80 0.36 



"^Intervals (c) and (d) arc best fit with additional cut-off powcr-law/powcr-law component 
^Units of lO"** erg cm"^ s'^ 



To illustrate these results, we apply them to Fermi observations of GRBs. 
For GRB090510 at z = 0.903 and = 1.8 x lO^^ cm, Ef'^2x 
10i^(/|/epfc) eV and rjP^ ^ 0.03/„5//3'i_2epfc(l - b), and is thus at the ~ 1 
- 10% level. For GRB 080916C at z = 4.35 and ^ 1.25 x lO^^ cm [HS], 
Ef^7x lO^^ri/epk) eV and t??^ ^ 0.0015/_6/r|t„ar(s)epfe(l - b). 

For GRB 09G926A at z = 2.1062 and d-L = 16.54 Gpc, the estimate above 
gives Ef ^ 10i7(r|/epfc) eV and t?^^ ^ 0M{f.e/6)/rii,arepk{l - b). The 
spectral parameters to derive the photopion efficiency are given for the Band 
function in Table [s] [139' . Here the variability time scale t^ar is scaled to 
0.15t„ar s, which is the characteristic FWHM time of the well-resolved pulse. 
One percent efficiency can easily become 100% if /s ss 0.3, and more for 
higher energy protons, making escape difficult. If long GRBs accelerate UHE- 
CRs, then the blast-wave Lorentz factor must be close to that given in the 
simple 77 opacity limit, otherwise the GRB would be highly radiative by 
hadronic processes. With the strong dependence on I^, photohadronic inter- 
actions become 100% efficient when F k, 300, so if GRBs accelerate UHE- 
CRs, an accompanying photohadronic 7-ray spectral component is predicted 
[14(11 [Mil [T421IT43] . 



Proton models 



The comoving synchrotron cooling timescale of an ion with atomic mass A 
and charge Z is given by (eqs. 7.49 - 7.50, 

_A^ fmpV GTrnieC 

Balancing the synchrotron energy-loss timescale with the shortest energy- 
gain timescale allowed in standard first- or second-order Fermi acceleration 
processes gives the relation 

A^ ( ra,\^ dTTiTLpC 



B'r^ -^A i^] (73) 



crT<; 



where (/> 1 is the number of radians over which a particle gains ^ e of its 
original energy. The observed peak synchrotron photon energy is 
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r^'syn ^ 3r Z fm^ \ B' 



(1 + z) 2(1 + z) A \m.pj B. 
r 27 m„ lO^ra A 



1' 



(l + z)08a/TOe (1 + Z)(0/1O) Z2' 



(74) 



Overlooking other limitations on particle acceleration, equation ( 74 1 shows 
that proton synchrotron radiation from GRB jets and blazars can reach « 50 
TeV {r w 10^) and « 5 TeV {F w 100), respectively. Proton synchrotron 
models have been developed in Refs. [1441 1145] for 7-ray blazars, and in [146] 
for GRBs. We consider application of this type of model to these two source 
classes. 

GRBs 

Eliminating 7' from these equations using the observer time tsyn = (1 + 
z)t'gy^/r for the emission to be radiated at measured energy nieC^esyn = 
rnieC^e'syn/i^ + ^) — lOOE'ioo MeV imphes a comoving magnetic field |147| 

[US] 



11 \ /Q \ 2/3 

1 + Z \ I OTTmeC ^ 



1.3 X 10^ /yl5/3\ / I 



rBf^jpCgyji J \^T^ syn 

3\ / 1 ^ , \ 1/3 



- t%l{s) [zv^J [rsE,oo) ' ^^^^ 

and an isotropic jet power, dominated by magnetic-field energy, of 

1 1 „3n6^2 n/2 o2 4 10/3 „3 /™ \ 10/3 

* *'''~2 ~2 (1 + Z)2 213/3 Z14/3 (1 + ^)4/3 

pl6/34.2 / fiTrm r\^/'^ 41^/3 7^16/3^2 (\ 

X ^-273% ( T^] ^ 2.5 X 10^^ ^473 ' 2/3 erg/s, 

(76) 

letting the blast-wave radius R = r'^ct^arjiX + A- For Fe (A = 56, Z = 26), 
the pow er r equirements are reduced by a factor k, 0.17. 

Eq. (76) shows that Lq oc 1^16/3. The absolute energy release is 

At ^ ^ , 2.5 X 1059 ^10/3 rfihlg,[s) ^ , , ^ 
^^'""^ = TT^ ^"^^ - (TT^ zw;^ k^M^ '''' ^ ^ 

where /f, is the beaming factor. Using a 100 kG magnetic field as a fiducial, 
so S = lO^iJs G, then the absolute GRB energy release is £*^ahs — 1-3 x 
10^^ B^At{s)tlgj.{s) fbF^ / {I + z)^ erg. Such large energy requirements pose a 
problem for highly magnetized GRB models. 

In the model of Razzaque et al. (2010) |146j . taking a lower limit /a = 0.5 
implied by naive 77 opacity arguments, gives the absolute energy requirements 
for GRB 080916C ^EM of 
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- 2 X lQ^^E-;^^\r:,/Q.bY^/^{e,/l Aeg)H%l{s) erg, 

after multiplying by a two-sided jet beaming factor = 1.5 x \Q~^(Oj/\ deg)^. 
The precise values oltyar, tsyn and At depend on model interpretation. For the 
external shock model [HH applied to GRB 080916C, the rapid variability may 
be made by irregularities in the complex surrounding medium, so t^ar = 0.1 s, 
the proton synchrotron cooling timescale to make the delayed onset is tsyn — 4 
s, and the GRB duration At = 10 - 50 s. 



From equation ( 77 ), strong proton/ion synchrotron radiations can be emit- 
ted in the Fermi range near 1 GeV or even 100 MeV due to strong cooling 
in an intense magnetic field, but the absolute energies are large without a 
very small beaming factor. The jet break time with apparent isotropic energy 
release « 2 x 10^'' ergs is t^r = 0.3(0^/1 deg)-^^/^no d. For such a narrow 
jet, the jet break would have taken place before Swift slewed to GRB 080916C 
at « To 17.0 hr [T5Dj . 

More complicated geometries might, however, relax the bulk Lorentz factor 
requirement further |151) . If the inner engine makes the prompt MeV radiation 
and residual shell collisions at larger radii make LAT 7-ray photons, then F 
could be as low as ^ 300 ^152]. In this case, the absolute energy release from 
GRB 080916C could easily be < 10^^ gj.g 

For a proton synchrotron model of GRBs to be viable [146] . a narrow 
jet opening angle of order 1° along with a value oi F ^ 0.5 gives £^,abs — 
few X 10^^ ergs, within ranges implied by interpretation of radio and 7-ray 
observations and beaming |153j . The detection of distinct components in GRB 
spectra suggests that a cascading interpretation be more carefully considered 
|154j . The local rate density of Type Ib/c supernova progenitors of long GRBs 
like GRB 080916C or low-luminosity GRBs like GRB 980425 or GRB 060218 
cannot exceed « 300 Gpc~^ yr~^ |155j . and the local beaming-corrected rate 
density of long duration GRBs is « 10 - 50 Gpc'^ yr"! [1561 HIT]. The local 
Type Ib/c rate is 91:1 X Gpc-3 yr-i [155 . This agrees with the star- 
forming galaxy SN Ib/c rate of ~ 0.28 century"^ per 300 Mpc^ per galaxy 
« lO'' Gpc-3 yr-i. These imply a beaming rate of « (10 - 50)/9,000 w lO^^ 
- 0.006. With the sa 1° opening angle required to explain GRB 080916C, then 
all Type Ib/c SNe would have to make GRBs, presenting another challenge 
for strong magnetic-field jet models of GRBs. 

Blazars 

Proton synchrotron and photohadronic models for blazars have been de- 
veloped, as noted, in |144l I145j . The protons and ions make a 7-ray com- 
ponent that originates from pion-induced cascades, and the synchrotron ra- 
diations of pions and muons make additional 7-ray emissions. Primary elec- 
tron synchrotron radiation is still usually required to make the nonthermal 
radio/optical/X-ray synchrotron blazar emission. 

Anita Reimer gives |158j a detailed fit to the famous HBL Mrk 421 at 
z = 0.031 at luminosity distance (11 — 134 Mpc. Mrk 421 is one of the first 
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Whipple TeV blazars and the first one with sub-hour [159} variabiUty detected. 
For the composite SED of Mrk 421 averaged between 2009 January 9 (MJD 
54850) to 2009 June 1, Reimer makes a proton fit with the parameters shown 
in Table El 



Table 6. Proton blazar model for Mrk 421 [l58] 



Parameter Value 



B' 50 G 

Su^r 12 " = r/12, B',o = B'/50 G, S12 = <5d/12 

4 X 10^* cm 
Ema. 7.2 X lO^^ZAaS^o eV" 



The Hillas criterion [160] for a single-zone blob model requires the comov- 
ing Larmor radius to be smaller than the comoving blob radius, that is, 

(78) 



With eq. ( 16 ) giving the maximum escaping particle energy Emax ~ rE'^^^ 
ZeB'r'i^, we see that this model blazar can accelerate super-GZK (i? ^ 5 x 10^® 
eV) particles if charged ions, but this is only possible for protons with slightly 
larger magnetic fields. These protons have maximum comoving Lorentz factors 

lp..ma. « 6.4 X . 

For particles well-trapped by the magnetic field of the plasma, Jp^max — 10^79 
with 7g ^ 1. 

For the 2-sided jet power as defined in eq. (47), taking (5d ~ gives 

Lj = 27rrf /Scr^ XI ' « = e, B , (79) 

i 

and the total Poynting power dominated by magnetic field energy is 

Lj^B = ^rf /3cr2B'2 ^ 4.3 x lO^^r^^B'^l erg/s , (80) 

which is not excessive, considering that the Eddington luminsoity is L^dd ~ 
10^ ''Mg erg/s, and A/g « 0.2 - 0.9 for the supermassive black hole power Mrk 
421. 

We still need to calculate the particle power and the efficiency for hadronic 
production of 7-ray emission. In the case of proton synchrotron radiation 
amounting to total energy fiux <P (erg/cm^-s), the particle power is 

Lj^p - znr^ pel [ — )-- - , [SL) 
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for a spherical comoving emission region. 
According to our standard approach [3], 

The proton-synchrotron energy loss rate is 

where U'^, = B'^/Sn. 

For a large range of proton spectra, the 7-ray power from photohadronic 
processes can be approximated by the number Np^Q of protons with Lorentz 
factor 7' w 7p , from which follows that 



max 



and 



where the bolometric 7-ray energy flux supplied by the proton synchrotron 
process is 10~^°^_io erg/cm^-s. The particle power is reasonable, even if 
7g — >■ 0.1. The proton synchrotron frequency, from eq. (74), 

^P.yn = ^^tI^ ^ |^7;^™a. = l-l X 10S^^i?^0^12, (85) 

so 

Ep^syn = 5.679' B50'^i2 GeV . (86) 



The equipartition field is, from eq. (49), or eq. (7.81) [3] 

Beq{G) ^ 130 «/7^13/7 1/7 ■ (^'') 

In the case of Mrk 421, j/ig = lO^'j/ij, = IQylio, id = (1 + 2)rft/c(5D = 
1200/(5i2 s, implying that 

.2/7 ,2/7 

Be,(G)-4.5^^^i4f , (88) 

'5l2J^17 

SO that this model is magnetically dominated with a magnetic field « 10 x Beq. 

The minimum bulk Lorentz factor from 77 constraints, using eq. (45), is 
r^^.^ = 20.3/_io[i;i/(10 TeV)]i/^ 13.8/_lo[^;l/(TeV)]l/^ for a maximum 
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measured photon energy Ei and using a variabiHty timescale of« lO'^ s. Note 
that /_io represents the i/F^, flux for the target photons and, though both 
have the same units, is a differential quantity compared to the bolometric 
energy flux ^. The condition for 7-ray transparency is somewhat inconsistent 
with the detection of multi-TeV photons from Mrk 421. 

The cooling regime for ultra-high energy protons through synchrotron cool- 
ing is determined by the ratio p = fgy^/t'^y^^, with strong synchrotron cooling 



when p <C 1 and weak cooling when p 3> 1 . From eq. ( 72 ) we have 



3(mp/me)3 140 
4ctt C/s-T^^b 79 

Thus this proton synchrotron model for Mrk 421 is in the weak cooling regime. 
Alternately, we can use the expression .146) for the saturation Lorentz factor 

3 X 10^° 

for which particles with larger Lorentz factors are in the strong cooling regime 
(here aj = /he w 1/137 is the fine structure constant). 




We can use eq. ( 71 1 to determine the photopion production efficiency for 
this model. For e^fc = 10'^''vn/1.24 x 10^° = 8.1 x IO'^vit, a ^ and 6 = -1, 
we obtain 

<^- ^^";r^-" ^ 4.3x10-^^, (91) 

taking = (1 + z)r^/c(5D — O.G13/(5i2 d. The comoving Lorentz factor of 
protons interacting with target photons with energy e'pf. is 7^^. = 6 x 10^. 
Protons with larger Lorentz factors interact primarily with the a = portion 
of the vFi, spectrum with frequencies less than Vpk- The efficiency of proton 
energy loss with 7' ^ 7^^, approaches 

VM) « vf, (^) - 7.4 X 10-3:^ . (92) 

For photopion losses (as for proton synchrotron losses) , the energy loss is most 
efficient for the highest energy particles, reaching ^ 1% for 7g ~ 1. Somewhat 
larger efficiencies are allowed — provided these allowances are consistent with 
the variability data — if the protons are trapped on long times compared with 
the crossing time, and the emission region is slow to expand. 

In conclusion, hadronic models for blazars and GRBs face no insurmount- 
able objections based on power or energetics. 



4.4 Cascade halos and the intergalactic magnetic field (IGMF) 

The magnetic field i?iGM in the IGM is bounded by a primordial field gen- 
erated by quantum fluctuations in the early universe or decoupling transi- 
tions of the fundamental forces |161) . Dynamo processes amplify the seed 
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fields. Gamma ray astronomy provides a method for measuring tlie inter- 
galactic magnetic field (IGMF) by exploiting the effects of the 77 —>-e~^ e~ 
pair-production process on 7 rays from extragalactic TeV sources interacting 
with soft photons of the EBL. 

When VHE 7 rays interact with ambient radiation fields and the EBL, 
then blazars, and therefore radio galaxies (i.e., misaligned blazars) are sur- 
rounded by anisotropic jets of relativistic leptons made when the 7 rays mate- 
rialize into energetic e+e" pairs. These leptons spiral in the ambient magnetic 
field to make extended synchrotron radiation halos |162| . and GeV 7 rays by 
Compton-scattering photons of the CMBR [163] . Arrival time information 
in pair cascades generated by impulsive sources of high-energy, multi-TeV 7 
radiation can also be used to infer the strength of the IGMF |164( I165j . 

Table 7. Derived Limits on Bigm for the source lES 0229-1-200 assuming persistent 
emission 



lES 0229+200 


di (rad) 


Bigm(G) 


Ref. llMj 


TT 


> 3 X lO-^'' 


Ref. [167] 


0.1 


> 5 X 10"^^ 


Ref. [I68] 


0.03 


> 2 X 10"^^ 


Ref. [Ml 


0.1 


> 5 X 10"^^ 



Because the emission of 1 - 10 TeV photons from a source at redshift 
z ^ 1 is attenuated by, primarily, the IR EBL, then the cascade spectrum 
can be calculated for a given spectral model of the EBL and properties of of 
-BiGM- For sufficiently weak magnetic fields, the pairs travel rectilinearly while 
Comptonizing CMB photons to 7 ray energies. The absence of this cascade 
signature in joint Fermi-HESS observations of candidate TeV blazars implies 
a lower limit on -Bigm under the assumption that the blazar persistently 
emits high-energy radiation for arbitrarily long time. Neronov & Vovk [166) 
and Tavecchio et al. [167| argued that nondetection of the TeV blazars lES 
1101-232 (z = 0.186), lES 0229+200 (z ^ 0.14), lES 0347-121 (z ^ 0.188), 
and H 2356-309 (z = 0.165) by Fermi implies a lower limit Bigm ^ 3 x 10"^^ 
G (Table [7]). 

The high-energy electrons and positrons also undergo deflections in the 
ambient magnetic field, so the emissions arriving latest generally come from 
leptons that have been most severely deflected. This will cause steady sources 
to be surrounded by an extended angular halo formed by leptons deflected 
back into the line of sight. Ando & Kusenko [170 ^ claim that ^ 30' halos are 
found in stacked Fermi data for 170 AGNs that are bright at 10 - 100 GeV. 
By fitting the angular distribution of the halo, they deduce that -Bigm ~ 
10^^^(Acoh/kpc)~^/^, for magnetic-field correlation length Acoh ~ 10 - 100 
kpc. This claim is disputed in [171) . in part because of the use of P6_v3 
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response functions which are known to have inaccurate PSF from on-orbit 
calibration data [?7] . 

More reahstic limits on Biqmf based on evidence that the blazar was 
operating only during the time frame over which it was observed, was proposed 
in |172j . and independently, in |169j . 

Cascade radiation from EBL attenuation of TeV photons 

Consider a source and observer separated by distance d, as shown in Figure 



13 Photons with dimensionless energy ei = hui/m^c emitted at angle 9i to 
the line of sight between the source and observer, travel an average distance 
A^-y = A-y^(ei,z) before materializing into an electron-positron pair via 77 
absorption with photons of the EBL. After production, the pairs cool by 
scattering CMB radiation, which is detected at an angle to the line of sight 
to the source when the secondary electrons and positrons (hereafter referred to 
as electrons) are deflected by an angle 0dfl- Thus 9i — 9dfi — 0. For the purposes 
of the argument, we neglect redshift effects for the TeV blazar sources under 
consideration, which limits the treatment to sources at z ^ 0.2 (see |173j for 
redshift corrections). 

The time delay At between the reception of photons directed towards the 
observer and those which undergo the process described above is given through 
the expression 

. , J \ , c'sin(6'dfl -6) 

cAt = \j + X — d = \j H — d = 

sm Udfi 

Xjj{l — cosOdti) — d{l — cos6) , (93) 
noting that x = d sin 9i/ sin ^dfl and A-y^ = d sin 6 / sin ^dfl • In the limit of small 



observing and deflection angles, equation ( 93 1 implies 



At^^elAi-^)- (94) 




Observer 



Fig. 13. Sketch of the geometry of the event. A photon with energy -EtoV TeV 
emitted at angle 61 < 9j to the line of sight, where 9j is the jet half-opening angle, 
interacts with an EBL photon to create an electron-positron pair. The electron is 
deflected by angle Odd and scatters a CMB photon to energy -EgoV GeV and angle 
6, which is detected as a source photon by the Fermi LAT when falls within the 
angular point spread function 6'ps/(GeV) from the source direction. 
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Moreover, this time delay is observed by photons received at 

0^9^^^ (95) 
a 

to the line of sight. 

The deflection angle depends on the Lorentz factor 7 = 10^76 of the pro- 
duced electrons, and can be written in terms of the received photon energy 
E — EqcV GeV. The average photon energy of the CMB at low redshift is 
eo ~ 1.24 X 10~^ in mgC^ units, so that the mean Thomson-scattered photon 
energy is ex ~ (4/3)eo7^- Thus, an electron with Lorentz factor 7 scatters 
CMB radiation to photon energy E when 70 ^ l.l\/i?GcV- The characteristic 
length scale for energy losses due to Thomson scattering is 

SmgC^ ,0.75, ^ , 
At = = ( ) Mpc , (96) 

4crTUCMB7 76 

where ucmb — 4 x 10""'^^ erg cm~^ is the CMB energy density at low redshifts. 
While losing energy, the electron is deflected by an angle 6b = At/?'l in a 
uniform magnetic field of strength Bigm — G oriented perpen- 

dicular to the direction of motion of the electron, where the Larmor radius 
Tl = nieC^^/eB = 0.55(76/B-i5) Mpc. Thus the deflection angle for an elec- 
tron losing energy by scattering CMB photons to energy i? in a uniform field is 
9b = At/j'l — l-lS-is/i^GcV- Introducing a coherence length Xcoh that char- 
acterizes the typical distance over which the magnetic field direction changes 
by w 7r/2, then the deflection angle 




(97) 



with w = H{Xcoh - At) +_vA^WXri?(AT ~ Xcoh)- 

The EBL model of [T71] for sources at z = 0.14 gives X-y-y{E) = 200 Mpc, 
125 Mpc, and 70 Mpc at E = 1, 3, and 10 TeV, respectively. A low EBL based 
on galaxy counts |175j gives values of Xyj{E) = 280 Mpc, 150 Mpc, and 85 
Mpc at E = 1,3, and 10 TeV, respectively. Thus we write A^^^ = lOOAioo 
Mpc, so that Aioo ~ 1 gives a minimum value of At and 9 for the reprocessed 
TeV radiation. 

For a source at distance d — dcpc Gpc, with dcpc ~ 1 corresponding to 
z ~ 0.2, the time delay for emission observed at angle 



from the line of sight is given from equation ( 94 ) by 
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Equation ( 99 1 shows that small time delays are implied when A-y^ is small 
or A^^ ~ d. When A-^-^ ^ At, an additional delay ~ •^T^dfi/"^ arises during the 
time that the electrons are losing energy and being deflected by the IGMF 
|1761 1140j . Such small values of A^-y ~ 1 Mpc are only relevant at low redshifts 
to ^ 100 TeV photons pair producing within « 1 Mpc of their source which, 
however, may be in the galaxy cluster environments where the magnetic field 
is not representative of the dominant volume of the voids. Thus we can dismiss 
such an origin of a short time delay without assuming special properties of 
the TeV sources. The case A^^ ^ d formally gives short time delays, but this 
corresponds to the case when the mean free path for 77 pair production is 
about equal to the source distance, which occurs for 7 rays with energies of 
several hundred GeV when d ~ 1 Gpc. In this case, the secondary electrons 
which take half the energy of the high-energy photon scatter CMB photons to 
<C 100 MeV. Even so, the attenuation of the high-energy photons takes place 
over the entire distance d, so that the upscattered photons with short delay 
time comprise only a very small fraction of the incident flux attenuated close 
to the observer. 

The remaining alternative to avoid assuming that TeV blazars are steady 
on timescales of millions of years is to suppose that either Bigm ^ 10~^^ 
G or that Xcoh ^ 1 Mpc. If Xcoh ^ 1 Mpc, this contradicts the claim that 
the IGMF has been measured to be Bigm ~ 10~^^ G. If Xcoh < 1 Mpc, then 
the field must be even larger in order that the electrons are deflected away 
from the direction of the photon source. Thus TeV blazars must be assumed 
to be steady emitters on long timescales. Here we relax this assumption, and 
suppose most cautiously that the blazar has been operative over the last few 
years that the they have been observed. This reduces the implied lower limit 
by several orders of magnitude, as will be seen. Here we treat the blazar as a 
point source; see [172] when the Fermi-LAT PSF constraint is included (Fig. 



13) 



Derivation of the cascade spectrum 

Consider a source that emits TeV-scale photons from within a photon beam 
of half-angle Oj, and an apparent isotropic spectral luminosity e*L*(e*) within 
the beam ("flat-topped" jet). If the TeV photons can escape from the nuclear 
environments without 77 absorption, then they may still be attenuated by 
interactions with EBL photons. After 77 attenuation, e^-e" pairs lose energy, 
primarily by Compton scattering the CMBR. This upscattered radiation will 
be detected by an observer until the cooling pairs are deflected out of the 
beam. 

We use the notation of Ref . [3 . From eq. Q , denote the vF^, spectrum by 
/e = I'K = , e = hv/nieC^ , = (1 + z)e. 



The correction due to EBL attenuation is 
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f^^ 4^^2 cxp[-T^^(£;z)j = ^— p exp[-T^^(e,z)J , (100) 

At low redshifts, z ^ 1, e* « e, and ei(e) = meC^ N {e) , where 
A'^(e) is the photon injection function. The number of photons surviving to 
the observer is 7V(e) exp[— T^^(e; z)], so the number of absorbed photons is 
Nabs{f) = N{e)[l — exp[— T^-y(e; z)], and this also represents the electron in- 
jection function Ninj{'ji), with normalization 



d-/^ N.njili) = 2 de Nabsi^) . (101) 
Jo 

because each photon makes two leptons with 7^ = e/2. So Ninj{'ji) — 
4iVa6s(e) = 4A^(e){exp[T^^(e, z) — 1)]}. Using equation ( 100 1 gives the injection 
function ^ 

N^nJh^) = ^^^^^ {expK^(e, z)] - 1} , e = 27, . (102) 

This injection source of e"'"-e~ loses energy by Compton-scattering photons of 
the CMBR to GeV energies, and this cascade GeV component is considerably 
dimmed when the leptons are deflected out of the beam. 

Photons with energy w 1 TeV make leptons with 7 10^ which scatter 
the CMBR to 7^60 ~ 10^, or to photon energies ^ 500 MeV. Scattering is 
in the Thomson regime for 4769 <C 1 or 7 <C 2 x 10^, that is, electrons with 
energies < 100 TeV. 

The Thomson energy-loss rate 

^7, 4 Mo 2 _ 2 /1nQ^ 

with uo/rrieC^ = 4.9 x 10^^ cm"'^. The solution to the steady-state electron 
continuity equation is 

1 r°° 

7V(7) = ^ / d7'iV(7')- (104) 

l^T7 J-y 

Limiting 7 > 7^// , where 7^/; is the deflection Lorentz factor where the lepton 
is deflected out of the jet opening angle by the IGMF. Taking the luminosity 



spectrum from Compton scattering [3], eqs. (102) and (28) give 



Und^caTUo ^e^y ^ FHqJ}) r /, {exp[r^^(e, z)] - 1} 



VTrrieC^ ^eo' J 7"* 



11 \ 

"7i 

7 



(105) 



At high energies, most of the energy is taken by one of the leptons. See, e.g. 

[ml. 
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with e = 27i. The interior integral is the injection function from high-energy 7 
rays absorbed by photons of the EBL, and depends on a model of the optical- 
depth function Tjj^e, z) for a photon with detected energy e that was emitted 



by a source at z. In general, one uses the expression eq. (31 ) for the Compton 



kernel Fc{q, F^), though the isotropic Thomson kernel, eq. (331, is sufficiently 
accurate our purposes here. 

The lower limit ^low in eq. ( 105 1 is the maximum of various constraints 



given by kinematic factors, engine duration, and 7d/;. The kinematic Lorentz 
factor 7fe„m given by eq. (35). In the Thomson regime, ^knm (l/2)^/es/eo. 



The deflection Lorentz factor is determined by the condition that the 
energy-loss timescale is equal to the timescale for e"'"-e~ pairs to be deflected 
out of the beam. The Thomson energy loss timescale — 7t 
tT = I/utI- The deflection timescale 



vt^'^ implies 



tdfi = = t'ji — , when At < Koh 

c eB 



(106) 



and tdfi = Oj [meC^/eB) ^ At l^coh when At > \coh- Solving gives 




B \2/3 



(107) 



This constraint implies 7d/; = 3.7xlO'^yB_i5/6'_i for 7 > 7.5x 10^/Aco;i(Mpc), 

and 7d/i = 6.2 x lO6A^/|^(Mpc)(B_i5/0-i)2/3 for 7 > 7.5 x lOVAcoh(Mpc). 

To avoid solving a time-dependent electron continuity equation, we intro- 
duce electron Lorentz factor limits 7(^i) to define the time the engine was 



operating. Following eq. (34 1, noting that /e 



,7d/i,7(zlt) 



7 



esL{es)/4:T:d^, gives 



1 



47260' 



X / d7i 
'7 



/, {exp[T^^(e, z)] - 1} 



(108) 



with e = 27i. 

The minimum Lorentz factor 7 related to period At of activity of central 
engine is determined by equating the time delay with the extra pathlength 
followed by photons. Thus 



At 



2~c ^''f' ' 



(109) 



and Xtot = A^-y(ei) -I- At (7) = lOOAioo Mpc, Odfi = w9b, Ob = Xt/tl- From 
this we derive 

j{At) = 
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eB 



Ai^^l/4 ^ 9.9x10«a;/>-15 7.5 xlO--' , ^ 

m^WT \2cAt) ^ [Zit(s)]i/4 ' Aeoh(Mpc) ^ ' 

/- A,„,A.„,. Nl/3 ^ 2.3xlO^^A^/gBy/,A;/f,(Mpc) 7.5x10^ 
V 2i'TAt ) — [zit(s)]i/3 ; 7 A,<,fe(Mpc) 

(110) 

Note that there is no separate constraint that the electrons have cooled 
long enough to scatter significant emission in the given band. In the station- 
ary frame, one can define the cooling electron Lorentz factor jcooi giving the 
characteristic Lorentz factor of electrons that have cooled in time At, ex- 
pressed in terms of the GeV photon energy EceV to which CMBR photons 
are Compton-scattered CMB photons. At low redshifts, 

. 3 17^ ^7.7x IQi^s _ 2.2 Myr 



AcaTUojcooiiAt) 2c(TtUo V e 7 VEceV ' 

(111) 

we see that it takes nearly a million years for the generated electrons and 
positrons to cool and begin to make their strongest emissions in the GeV - 
10 GeV band of Fermi. But this applies for pair halo emissions made at large 
angles to the jet axis. The cooling timescale compared to the engine starting 
time, as measured by an observer within the jetted emission cone, is instead 



given by equation ( 109 1 , which we write as 



Atis) ^ 60 Aioo T^w^ . (112) 



The lack of distinct GeV echoes from impulsive or flaring high-energy sources 
can limit parts of the BiGM-Acoh parameter space. Better yet, the discovery 
of such echoes in GRB light curves, as originally proposed by Plaga (1995) 



|164| . would finally reveal the primordial magnetic field. From equation ( 111 ), 
the Thomson cooling length of relativistic electrons scattering the CMBR is, 
at redshift z, given by 

2.3 X lO^Ocm 0.7 Mpc 
7(1 + ^)^ "71^' ^'''^ 
where the last expression applies at low redshifts. 



Fig. 14 shows calculations with this semi-analytic model PJ2j- Here we 
use a photon-energy dependent expression for X^j and perform an integration 
over the CMBR spectrum. If the jet is persistent on long time scales, then 
the jet opening angle is limited to be ;$ 0.4, as can be seen in the middle 
panels of this figure. Restricting TeV activity of lES 0229-H200 to « 3 - 4 
years, during which the source has been observed, leads to a more robust lower 
limit of BiGM ~ 10~^^ G, which can be larger by an order of magnitude if the 
intrinsic source flux above « 5 - 10 TeV from lES 0229+200 is strong. 

If there were no intergalactic magnetic field at all, then the pairs made 
from a source at distance d, minus the distance A.^.^ over which they are 
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Fig. 14. Model of cascade radiation spectrum, equation (100 1, applied to HESS, 
VERITAS, and Fermi observations of lES 0229+200, using model spectra (soHd 
curves) and EBL model of |174) to give attenuated source spectra, (a) Cascade 
spectra for lES 0229+200 assuming persistent TeV emission at the level observed 
with HESS and VERITAS, for different values of Bigm and \coK = 1 Mpc (solid) 
or AcoJi = 100 kpc (dot-dashed) for a jet opening angle Qj =0.1. The intrinsic TeV 
source spectrum is given by power-law with vF^ index = 4/5 with super-exponential 
cutoff oc exp [-(S/5 TeV)^]. The PSF constraint for the Acoh = 1 Mpc case is shown 
by the dot-dashed curves, (d) Same as for (a), except that Qj = 0.3. (e) Same as 
(a), except that = 1.0. (b), (c), (d): Here the TeV engine operates for 3 yr, 3 yr, 
and 100 yr, respectively, but the intrinsic TeV source spectrum differs. In (b), it is 
the same as (a). In (c), it has an exponential cutoff oc exp(— (i?/10 TeV). In (f), the 
intrinsic source spectrum use parameters of Dolag et al. (2011) |169j . In (b), (c), and 
(d), \aoh = 1 Mpc and cascade spectra are calculated for different values of Bigm, 
as labeled. 
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made, will have cooled to electron Lorentz factors 7 given by d — A^-y — 
At = 750 Mpc/(7/10^). It is interesting to think that the pair injection pro- 
cess will have seeded the voids of intergalactic space with an ultra-relativistic 
nonthermal electron-positron component that has already cooled to low ener- 
gies. Given how active TeV sources are in our present declining phase of the 
universe (in terms of star formation activity), TeV sources must have been 
vigorous operating in the early universe. 

Without making severe assumptions about the 7-ray duty cycle and ra- 
diative behavior of blazars like lES 1101-232 and lES 0229-^200 on long 
timescales, the best limit to the strength of the IGMF is ^ 10~^® G for 
Xcoh ^ 1 Mpc [nil [1691 • Constraints on the value of Bigm ;S 3 x IQ-^^ G can 
be obtained from a search for pair echos in the analysis of GRB data [17811179] . 
so if the larger field is correct, then no GeV echo radiation is predicted from 
GRBs. 



5 7 Rays from Cosmic Rays in the Galaxy 

Cosmic rays are the most energetic particles in the universe, and sources of 

• the light elements Li, Be, B; 

• the Galactic radio emission; 

• the Galactic 7-ray emission; 

• Galactic pressure; 

• collisional excitation of atoms and molecules; 

• terrestrial -'^^C and ^°Be, with half lifes of sa 5700 yr and « 1.5 Myr, 
respectively; and 

• astrobiological effects. 

Cosmic rays are composed mainly of protons and ions, but also include en- 
ergetic electrons, positrons, and antiprotonsj^ They make up an important 
particle background in the space radiation environment and contribute to the 
space weather. Cosmic-ray electrons emit radio and X-ray synchrotron ra- 
diation. X-ray and 7-ray bremsstrahlung radiation, and X-rays and 7 rays 
from ("inverse") Compton scattering [180j . It is believed with good reason 
that Galactic GeV-PeV cosmic rays are accelerated by supernova remnants 
[1811 1182| . The origin of the UHECRs, and its relation to the origin of the 
cosmic rays, is an open question; the hypothesis that their origin involves 
rotating black holes is developed in _3]. 

The difficulty to solve the problem of cosmic-ray origin for the '-^GeV 
cosmic-ray protons and ions that carry the bulk of the cosmic-ray energy con- 
tent is that, being charged particles, they do not point back to their sources 
as a consequence of intervening magnetic fields that deflect them in transit. 



High-energy neutrinos could also, depending on definition, be included. 
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Sites of high-energy particle interactions are identified by 7 rays, but ascer- 
taining whether the emission is made by hadronic cosmic rays is comphcated 
by the possible leptonic origin of most of the 7 rays. GeV-PeV neutrinos, by 
comparison with charged cosmic rays, unambiguously point to the sources of 
the cosmic-ray hadronic interactions, but are faint and difficult to detect. 

Cosmic-ray models use strong nuclear interaction cross sections and ac- 
celeration, loss and transport physics to derive the cosmic ray intensity and 
spectrum, whether in supernova remnants, diffuse and extended clouds of gas, 
galaxy clusters, the Galactic Center region, or wherever there is a significant 
column of gas with an illuminating cosmic-ray flux. The Sun and Solar flares 
are especially instructive for study of particle acceleration, transport, and 
radiation physics from 7-ray observations. 

5.1 7 rays from Solar system objects 

Besides direct observations, the effects of cosmic rays in the Solar system are 
traced via the cosmic-ray induced 7-ray flux of the Sun, moon, and Earth. 
In the Solar cavity, the cosmic-ray intensity is modulated by the outflowing 
Solar wind plasma, making an anti-correlated decrease and lag in the ^ 10 
GeV/nucleon cosmic-ray flux reaching in the Solar cavity with a period of the 
11-year sunspot cycle (one-half the 22 year Solar cycle). Interaction of these 
cosmic rays with Solar system objects make GeV 7-ray flux that varies on 
these timescales. With 3 years of data taken with Fermi, this only amounts 
to ;S 30% of the sunspot cycle, the first 2.5 yrs of which were taken with the 
Sun in a deep and extended Solar minimum. 

Due to its proximity, the Earth is the strongest 7-ray source for the Fermi- 
LAT, which is why zenith-angle cuts are made on source spectral recon- 
struction to eliminate the interfering effects of the cosmic-ray induced 7-ray 
"albedo" emission (recall footnote [5]). The 7-ray spectrum of the Earth albedo 
depends on angle to the nadir. At Fermi-LAT's h = 565 km orbit, it views the 
horizon at 9ndr ~ arcsin(l-|-/i/rE)~^ = 66.5°; the Earth's radius is — 6378.1 
km. Fermi is therefore exposed to (1/2)(1 — cos(7r — 9ndr)) ~ 70% of the full 
sky, with the Earth occulting « 30% of the sky. 

Observations of Earth albedo flux reveals a number of interesting things. 
One is the exposure bias toward the North Ecliptic Pole due to favored rock- 
ing to the North. Recall that Fermi rocks about zenith, early in the mission 
by 39°, which was increased to 50° later. Because of the intensity of the 
Earth albedo, a standard analysis cut is to accept photons only within 105° 
of zenith. For studies of the albedo, then, an acceptance of Ondr ^ 80° gathers 
mostly cosmic-ray induced terrestrial 7-ray emission, i.e., albedo. Knowing 
the albedo spectrum gives, empirically, a better characterization of the 7-ray 
background, and is useful for Fermi-LAT calibration. Deconvolving the 7-ray 
spectrum with an Earth atmosphere model and knowledge of the interaction 
cross sections would give the primary cosmic-ray spectrum and information 
about the deflection of primary cosmic rays by the geomagnetic field. 
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Fig. 15. Exposure maps (top) and intensity maps (bottom) of terrestrial 7-ray 
albedo emission 



In [183] . two data sets are gathered: one during the first 90 days of the 
Fermi mission, and a second for a two-orbit pointing at the Earth's limb. The 
analysis covers a total of 6.4 x 10^ events giving the albedo spectrum in the 
100 MeV - TeV range, with 218 ph(> 100 GeV), and 16 ph(> 500 GeV). 

From two-dimensional intensity maps with increasing energy range, a 
beautiful high-energy ring forms above 3-10 GeV due not to Fermi's energy- 
dependent PSF but to beaming of the emission during formation (see Fig. 
15). The bright limb at the Earth's horizon are 7-ray light showers made by 
grazing incidence cosmic rays coming directly towards the LAT. 

The next interesting feature of albedo emission is the angle-dependent 
spectrum. At the nadir, the spectrum is soft, and generated by 7-rays that are 
backscattered by large angles, which is considerably less likely due to threshold 
and beaming effects at 1 - 10 GeV/nuc, where the particle Lorentz factor 
7 ^ a few. The deflection of cosmic rays by the geomagnetic field makes an 
east-west effect that fades out at high energies. The spectral intensity rises 
and hardens until the Earth's limb is reached. See |183j for details. 

For the moon, without an atmosphere, a steep secondary nuclear produc- 
tion spectrum with flux F_8 = 110 ± 20 |184j is made as cosmic-ray GeV 
protons and ions impact the surface of the moon, confirming the EGRET de- 
tection |185j . A search for 7-ray emission from asteroid populations and other 
Solar system rocks and dust |186[ I187j is currently in progress. 

The 7-ray emission from the Sun consists of two components: cosmic rays 
impacting the surface of the Sun to make an albedo-type emission, and cosmic- 
ray electrons Compton-scattering solar photons to 7-ray energies [1881 1189j . 
Both components of emission are sensitive to the phase of the Solar cycle. 
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and both components have been detected. In analysis of 18 months of data 
|190| . the solar disk emission is found at the level of w 46. The measured 
integral flux of the extended non-disk emission from a region of 20° radius 
centered on the Sun is F-s ~ 70. So the Sun and moon are really bright 7-ray 
sources, and one has to be alert to data contamination and spurious variations 
when the Sun or moon drift past. 

The Fermi-LAT measurement of the Solar emission confirms Seckel's model 
|191| . but at a much higher (by a factor w 7) flux. The disk emission is 
practically flat in a i^F^, spectrum up to « 8 GeV, where it begins to fall 
off. The observed spectrum and angular profile is in good agreement with 
theoretical predictions for the quiet Sun emission |188j . 

The first Fermi-LAT as well as nuclear 7-ray line flare in Solar Cycle 24 
is the 2010 June 12 M3 flare [HI]. The LAT emission lasted for only 50 s, 
compared to past long-duration Solar 7-ray flares. Combined Fermi GBM and 
LAT analysis reveals a rich 7-ray line spectrum superimposed on a continuum 
with a hard-energy tail that is consistent with either a separate nonthermal 
bremsstrahlung component, or pion production in energetic nuclear events. 
The hard bremsstrahlung X-rays and the ~ 300 MeV 7 rays flare up within 3 
s of each other, placing strong requirements on acceleration theory and target 
properties. As we enter the active portion of the Solar cycle, a wealth of new 
Fermi data on particle acceleration physics of large Solar flares is anticipated. 

Also interesting and timely is the use of Fermi for the study of terres- 
trial 7-ray flashes [193j that accompany thunderstorms found mainly in the 
sub-tropics and tropics. Experiments in nadir-pointing modes to increase sen- 
sitivity to TGFs are currently underway. Interesting analysis effects having to 
do with the shortness of the pulse duration and the strong soft X-ray emission 
affecting the ACD, which is also important for Solar flares, reveals Fermi's 
capability for this science. 



5.2 GeV photons from cosmic rays 

The cosmic-ray induced 7-ray glow of the Milky Way is the most pronounced 



and distinctive feature of the Fermi sky (Fig. 16). A dominant fraction of the 
Galactic 7-ray emission is believed to be truly diffuse, and made by cosmic- 
ray bombardment of gas and dust in the interstellar medium (ISM). The most 
important hadronic process for 7-ray production is secondary nuclear produc- 
tion from the collisions of cosmic-ray protons and ions with ISM particles. The 
most abundant secondaries are pious (others are kaons and heavier baryons 
and baryonic resonances). The pious decay according to the scheme 

tt" ^- 27 , 7r+ ^ /i+ + ly^, , TT^ ^ fj.^ + Pf, , (114) 

/i+ ^ e+ + + i^e , fi^ ^ + 1^^ + i?e . (115) 

For the purposes of Galactic cosmic-ray physics, the decays are essentially 
instantaneous, and result in 7-ray, electron, and neutrino injection emissivities 
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Fig. 16. Fermi-LAT 2FGL 100 MeV - 10 GeV all-sky map, using 2 years of sky 
survey data (2008 Aug 4 - 2010 Aug 1) HOj. Scale in units of 10^'^ erg/cm^-s-sr. 
The 2FGL consists of 1873 sources, with 1170 associations and 127 identifications, 
of which 1319 are at |6| > 10°. 



(as defined below, eq. 118) proportional to the cosmic-ray intensity and gas 
density at that location. The secondary electrons and positrons, as well as 
those accelerated at cosmic-ray sources, lose energy during propagation mainly 
by synchrotron and Compton losses at high energies, bremsstrahlung losses 
at intermediate energies, and ionization and Coulomb losses at low energies, 
and can also in principle be reaccelerated. 

The diffuse Galactic 7-ray glow is the superposition of all the radiations 
made by 7r°-decay 7 rays, 7 rays from cosmic-ray electrons that Compton- 
scattered the available radiation fields, and electron bremsstrahlung 7 rays. 
Electron synchrotron radiation would only contribute to the Galactic 7-ray 
emission from localized sources, such as pulsar wind nebula. 

A simple expression for the demodulated cosmic-ray proton intensity in the 
local interstellar space inferred from measurements of the near-Earth cosmic- 
ray intensity is 

Jp{Ep, Hp) = 2.2E-^ '^^ CR p/cm^-s-GeV-sr (116) 

|194j . so that the GeV cosmic-ray flux represents tens of cosmic-ray protons 
per cm^ per s. The cosmic-ray kinetic-energy density is dominated by the 
kinetic energy Tp(= Ep — nipC^) of ~GeV protons, and is given by 

47r r°° , 
ucRKE = — dTp Tp Jp{Ep, Hp) = 0.7 eV/cm^ , (117) 
c Jo 

to which ions contribute another factor ^30 - 50%. The dominant elementary 
hadronic process is p+p — >■ tt^'^ , so astrophysical studies of secondary nuclear 
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production focus on this process. For more detail on radiative processes and 
cosmic-ray propagation, see |1951 11961 11971 1198j . 

Background modeling 

Reliability of source detection is improved if the intensity of the diffuse back- 
ground radiation is known, because a clumpy gaseous structure illuminated 
by cosmic-ray induced 7-ray fluorescence could be mistaken for a point source. 
Moreover, knowledge of the background is required to recognize dim sources, 
because the background dominates for all but the brightest sources. The back- 
ground model for Fermi analysis takes into account diffuse Galactic 7 rays 
from interactions of cosmic rays with material found in the various phases of 
the ISM, including the neutral hydrogen, HI, molecular hydrogen, II2, ionized 
hydrogen, HII or 11+ , and the dark gas phase |199j . The distribution of neutral 
atomic hydrogen (HI) is traced by 21 cm line surveys. The molecular hydrogen 
(H2) distribution is derived indirectly, most commonly by using 2.6 mm line 
observations of carbon monoxide (CO). (The ratio of H2 to CO — the so-called 
'X'-factor, X = A^(i?2)/M^co[^is derived from 7-ray observations.) The total 
atomic and molecular gas column density can also be traced indirectly from 
extinction and reddening by dust, depending on the relative fraction of dust 
and gas. Cosmic rays also interact with ionized hydrogen. The low-density 
ionized gas can be inferred from dispersion measures of pulsar signals in the 
radio band. 

Cosmic rays move in large-scale galactic magnetic field and diffuse by scat- 
tering off magnetic turbulence. In the thick disk of the Galaxy, where the bulk 
of the diffuse 7-ray emission is made at |6| ^ 5° - 10°, the cosmic-ray intensity 
changes with galactocentric distance since the source distribution is peaked 
at star-forming arms at 4 6 kpc. The cosmic-ray intensity also changes with 
the distance from the Galactic plane because of escape, and because the diffu- 
sion coefficient is energy dependent, with high-energy particles diffusing faster 
through the Galaxy. In addition to the CMBR, the different distributions of 
background optical and IR fields, and location-dependent Galactic magnetic 
fields, means that cosmic-ray electrons suffer position-dependent energy losses. 

For the electronic component, gradients in the particle distribution can 
be severe and contribute a Galactic background contribution to the 7- 
ray emission. The GALPROP (GALactic cosmic ray PROPagation) model 
[2011 ISni HOI Uni, started in 1996 and developed independently of Fermi, 
determines cosmic-ray diffusion coefficients from fits to cosmic-ray data. The 
spatial and momentum diffusion equations for cosmic-ray transport are solved, 
taking into account source injection, energy and fragmentation losses, and en- 
ergy changes for cosmic-ray protons, ions, and electrons. The stellar optical 
field, assorted IR and PAH lines in the 10/i valley, and a FIR dust peak at 



N{H2) is the column density of molecular hydrogen, and Wco is the brightness 
temperature of CO integrated over velocity |200j . 
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~ 100^ |202j provide target photons to be Compton scattered by relativistic 
electrons. The GALPROP model is constrained by the energy dependence of 
the B/C and ^Be/^°Be ratio, from which predictions for the e+, e~, p, and 
7-ray spectra and intensity can be made. 

Template modeling of 7-ray emission from Gould belt clouds in Cassiopeia 
and Cepheus |205j shows a weak Galactocentric gradient from the Gould belt 
to the Perseus arm, and an increase in the X-factor from = 0.87 x 10^" cm^^(K 
km s~i)~i in the Gould belt clouds to = 1.9 x 10^° cm'~^{K km s^^y^ in the 
Perseus arm. The dark gas represents « 50% of the mass traced by CO. 



Diffuse Galactic 7 rays from cosmic rays 

There can be little doubt that cosmic-ray interactions make a large fraction of 
the 7 rays observed with Fermi. This is established most clearly by Fermi-LAT 
observations [2(76j towards a region in the third quadrant between Galactic 
longitude 200° ~ 260° and latitude 22° - 60° that contains no known molecular 
clouds. After subtracting point sources and Compton emission, the residual 
100 MeV - 10 GeV 7-ray intensity exhibits a very strong linear correlation 
with atomic gas column density. 

According to the model of Cordes & Lazio (2002) [207], the N(HII) column 
density is only (1 - 2) x 10^° cm~^ and fairly smooth in the third quadrant 
regions analyzed in Ref. [206] . The contribution from ionized gas has a small 
effect on the measured emissivity, as the N(HI) column density ranges from (1 
- 13) X 10^° cm~^. The measured integral 7-ray emissivity measured with the 
Fermi LAT is 1.63±0.05 x 10~26 photons/s-sr-H-atom and 0.66±0.02 x lO^^e 
photons/s-sr-H-atom above 100 MeV and above 300 MeV, respectively, with 
an additional systematic error of ~ 10%. 

These numbers are explained in the first approximation if cosmic rays 



pervade the gaseous disk of the Milky Way with the same intensity, eq. ( 116 1, 
as observed locally. The pion production rate per unit volume, or differential 
photon emissivity, is 

npH^TT^yF^) = 47rnH / dTp Jp{Tp,f2p) — , (118) 

Jo '^-'tt 

and the bolometric tt" 7-ray production rate per H atom is, therefore 



97 = dtdVdQ ^^^" - 2C / dTp Jp{Tp, Qp)(Tpp^^o(Tp) , (119) 



including a factor for two photons per tt*^ and a metallicity correction C, ~ 1.5. 
Between Ep k, 1.3 GeV, just above threshold, and Ep « 10 GeV, the inclusive 
t:^ production cross section is approximately linear, and can be written as 

crpp^^o(mb) 4 mb{Ep/l.i GeV) . 
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from which, with eq. (116), we find that 



9^ 2.2 X 10^26 [s-sr-H-atom]-! . (120) 

roughly agreeing with the Fermi-LAT measurements |206| quoted above. 
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Fig. 17. Galactic diffuse emission intensity as measured in the latitude range 
10° < \b\ < 20° |39) . (left) LAT Galactic diffuse 7-ray intensity data are given 
by red dots and red cross-hatched error regions; EGRET data are given by the blue 
crosses and blue shaded regions, (right) LAT data compared to a model consisting 
of tt" decay 7 rays (red); bremsstrahlung from primary and secondary electrons and 
positrons (magenta); Compton-scattered soft photons to 7-ray energies by cosmic- 
ray electrons (green); unidentified background consisting of point sources, isotropic 
diffuse Galactic, and the isotropic extragalactic 7-ray background. 



The intermediate latitude, 10° < |6| < 20°, Galactic diffuse 7-ray emission 
intensity spectra as measured with the Fermi-LAT and EGRET are shown 



in the left panel of Fig. 17 |39]. The prominence of the 7r° decay feature 
is apparent, as also are the large systematic differences at both low, ;S 300 

the LAT 



17 



MeV, and high, ^ 300 MeV energies. In the right panel of Fig. 
spectrum shown in the left panel is compared with spectra from an a priori 
model (based only on local cosmic-ray data) for the diffuse Galactic 7-ray 
emission, updated from GALPROP. As the Fermi LAT accumulates data, 
and analysis becomes even more accurate below 100 MeV for the diffuse class, 
we can foresee using Fermi data to determine the best interstellar cosmic-ray 
proton and He ion spectra rather than the other way around. 
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Other Fermi LAT cosmic ray results 

CR e spectrum: The Fermi LAT can be used as a cosmic-ray electron detector 
by using its anticoincidence dome to identify incoming charged particles, and 
distinguishing between leptons and ions from the tracks in the tracker and the 
calorimeter. The Fermi telescope has an acceptance of > 2 m^-sr for combined 
cosmic ray e+ and e~ (CR e), since it cannot distinguish charge. 

The extraordinary statistics of the Fermi enabled a precise measurement of 
the CR e spectrum, showing a featureless spectrum consistent with a power 
law of number index = -3.04 between w 25 and 900 GeV [105]. The CR 
e spectrum is harder than a GALPROP model prediction where the diffu- 
sion coefficients tuned to the low-energy electron spectra are extrapolated to 
higher energies. The measurement disagrees with the ATIC report [209] . The 
PAMELA measurement of CR e+ |210| and the cosmic-ray electron spectrum 
inferred at TeV energies with HESS data [211] . in addition to the Fermi mea- 
surement of the cosmic-ray electron spectrum, place constraints on cosmic-ray 
electron models that require either (or a combination of) modifications to the 
propagation characteristics; local sources, most likely pulsars [2121; or contri- 
butions from dark matter (see lectures by Prof. Bergstrom). 

The 2010 Fermi analysis gives the CR e spectrum from w 7 GeV to 1 TeV, 
with a slight hardening above 100 GeV and a softening above 500 GeV |213j . 

EGRET excess: This term refers to EGRET measurements of the diffuse 
Galactic 7-ray intensity that were found, irrespective of direction, to be in 
excess of that predicted using the local demodulated cosmic ray spectrum 
and measured target gas mass "214'. Possible explanations included an 

1. unusual location and local cosmic-ray spectrum measured here at Earth; 

2. nuclear physics wrong; 

3. addition of 7 ray signal from annihilating dark matter; or 

4. EGRET miscalibration. 

With the launch of the Fermi Gamma ray Space Telescope, measurements [39] 
of the diffuse Galactic emission have been found to favor the latter hypothesis, 
namely, that EGRET was poorly calibrated above « 5 GeV. The differences 
in the EGRET and Fermi LAT diffuse Galactic 7-ray emission are shown in 



the right panel of Fig. 17 



The GALPROP model decomposes the Galactic plane emission into pionic, 
electron bremsstrahlung and Compton fluxes, point sources, and an isotropic 
diffuse background. The tt*^ — )■ 27 signature is clearly seen in the spectrum 
of the diffuse Galactic 7 radiation [31]. Discrepancies between EGRET and 
Fermi spectra still remain at ^ 100 MeV, where systematics effects become 
severe for both the EGRET spark chamber and the LAT tracker; see Fig. 
[TT] The latest analysis of diffuse Galactic 7 radiation in the Third Galactic 
Quadrant finds weak evidence at best for a Galactocentric gradient in the 
cosmic- ray intensity |215| . 

Fermi has also measured the positron flux and fraction using the Earth as 
a magnetic field |216j . 
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5.3 Fermi bubbles 

It is not possible to cover, or even mention, all Fermi-LAT diffuse studies in 
this short review, e.g., the Galactic center region, the Cygnus Loop |217j . and 
7-ray emission from massive star clusters and OB associations like the cosmic- 
ray filled cocoon region in Cygnus |218j , but we would be remiss not to mention 
the Fermi bubbles that D. Finkbeiner and colleagues find |2191 1220] , In their 
analysis, backgrounds determined by template fitting are subtracted from the 
7-ray data, leaving (apparently) well-defined structures — the Fermi bubbles — 
at ss 1 - 10 GeV, and which probably extend to much higher energies. They 
are symmetrically arranged north and south of the Galactic Center with a 
width of « 40°, and extend some ~ 50° in height. They are correlated with 
the WMAP haze, which is excess 20 - 40 GHz emission found after correcting 
for the dipole anisotropy and subtracting out emission traced by Ha and soft 
synchrotron radiation extrapolated from the 408 MHz Haslam survey |221j . 
The 1-100 GeV luminosity of both bubbles, if emanating from the Galactic 
Center region, is « 4 x 10'^'' erg/s. With a minimum lifetime of w lOkpc/c, its 
minimum energy content is 4 x 10^^ erg. The bubble luminosity represents 
- 5% of the 100 MeV - 100 GeV luminosity of the Galaxy [2^ . 

The haze itself, if not an artifact of template fitting or mis-extrapolation 
of the radio emission [223j . could be formed by dust, spinning dust, or dark 
matter, but a synchrotron origin seems favored, especially given the Fermi 
bubbles. Indeed, a leptonic model for the bolometric SED of the bubbles, with 
synchrotron radiation for the WMAP haze and Compton-scattered CMB 7 
rays for the Fermi-LAT emission, can be arranged by tailoring the electron 
distribution ^220j . A detailed model including inverse Compton scattering off 
the CMB, FIR, and optical/UV radiation fields fits the average LAT 7-ray 
spectrum |224| . which is practically flat between 1 and 100 GeV. There is 
furthermore no energetics problem in a leptonic model, but it requires reac- 
celeration and a low-energy cutoff in the electron spectrum. The sharp edges 
reported for the Fermi bubbles, to be confirmed in ongoing Fermi-LAT analy- 
ses, stand in contrast to the the WMAP haze. A specific hadronic model [225j 
suffers from weak radiative efficiency in the dilute wind from the Galactic 
center, and lateral diffusion would smear the edges of the bubbles. 

More fundamental, perhaps, is what the existence of the Fermi bubble 
says about our Galaxy. Was it a giant Galactic explosion, the residuum of a 
more dynamic period in the history of our Galaxy, a superwind from a past 
starburst episode, or something else? 

5.4 7-ray supernova remnants 

The cosmic-ray power Pcr can be estimated as 
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~ 10''" erg/cm^ , (121) 

using a residence time of ~ 20 Myr from analysis of ^"Be abundance in cosmic 
rays [THS], and the Milky Way volume Vmw ~ 7r(200 pc)(15 kpc)^ - 4 x 10^6 
cm^. Thus the sources of the Galactic cosmic rays need to supply ^ 10^° erg/s 
in the form of nonthermal particle power throughout the Galaxy. One Galactic 
supernova every thirty years supplying « 10^^ erg in outflowing kinetic energy 
will inject with 10% efficiency w 10*^ erg/s of cosmic-ray power, so significant 
inefficiency can be accommodated if SNRs accelerate the cosmic rays. For the 
cosmic-ray power implied by Fermi-LAT 7-ray observations of the Milky Way 
for leptonic and hadronic cosmic-ray models, see |222j . 

The substantial time-averaged kinetic powers of SNRs make them the fa- 
vored candidate source population for the hadronic cosmic rays in our Galaxy. 
This power is greater than the time-averaged power available from, for exam- 
ple, pulsars, galactic X-ray binaries, stellar winds, or novae. In the electromag- 
netic window, the confirming signature of enhanced cosmic-ray activity is the 
7r° 7-ray bump [18111182] peaking at to^o /2 « 70 MeV in a number spectrum, 
and at several hundred MeV (depending on the spectrum) in a vFy represen- 
tation. The confirming signature of tt" production at ^ 200 MeV competes 
with systematics and inaccuracies in the background model which rise due to 
the large LAT PSF in the lowest decade of its energy range (Table [1]). 

Nearly 300 SNRs are known, mostly through radio detections. The Chan- 
dra X-ray catalog contains nearly 100 SNRsj^ and ROSAT Galactic SNRs 
must number several dozen. The poor EGRET PSF made it impossible 
to identify specific SNRs with sites of enhanced > 100 MeV emission, 
though several associations could be made, including IC 443, W28, and W44 

In the IFGL, 41 SNRs are associated with LAT 7-ray sources, and 3 
are, through morphological features, identified with SNRs (Table [2]), namely 
W51C, W44, and IC 443. By the time of 2FGL, 62 associations and 9 (mor- 
phological) identifications of SNRs and PWN with Fermi LAT and other 
wavelengths have been made (Table [2]). The identifications represent 6 SNRs, 
adding also W28, W30, and the Cygnus Loop |217j . and three PWNe, namely 
J0835.3-4510 (Vela X), J1509.6-5850 (MSH 15-52; possibly the remnant of a 
supernova in 185 CE), and J1826.1-1256 (HESS J1825-137). A fourth PWN, 
MSH 15-52, was found in a search of off-pulse emission of PSR B1509-58 [^ . 

Some important questions to be answered are 

• Are cosmic-ray protons and electrons accelerated at SNR shocks, and with 
what relative efficiencies? 

• With what efficiency is shock kinetic energy converted to cosmic-ray en- 
ergy? 

• What is the spectrum of accelerated particles escaping from the shock 
acceleration site? 
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• What is the maximum particle energy, and do SNRs accelerate CRs up to 
or beyond the knee of the cosmic-ray spectrum? 

• Does magnetic field amplification take place in the vicinity of SNRs? 



Table 8. Historical Supernovas |230) 



SN 


Type" 


Distance 
(kpc) 


Size*" 


7-ray detection'^ 
GeV/VHE 


Pulsar'* 


1987A 


II 


55 


0.8" 






1680; Gas A"" 


II 


3.4 


5.9' X 5.5' 


Y/y 




1604; Kepler 


la 


5 


4.5' X 3.7' 






1572; Tycho 


la 


2.4 


8.7' X 8.6' 


Y/y 




1181; 3C58 (?) 


II 


2.6 


6' X 3' 




J0205-H6449 


1054; Crab 


II 


2 


2.3' X 2.1' 


y/Y 


Crab 


1006 


la 


1.6 


32' X 26' 


?/Y 




393; RX J1713.7-3946 


II (?) 


1 


^ 40' X 40' 


y/Y 




386? 












369? 












185; MSH 15-52 (?) 


II 






y/Y 


1509-58 



''Chandra X-ray sizes^^ except for SN 1987A (HST), RX J1713. 7-3946 [233], and 

3C 58 (radio; [^) 

'^Flux measured at GeV or VHE energies; Y: larger i/Fi, flux, y: smaller i/Fi, flux 

"^See Tabled 
''Detection by Flamsteed disfavored |230| 



Fermi has now detected historical/young ( ^ 3000 yr; see Table |8| SNRs, 
namely Cas A (1680 CE), Tycho (1572 CE), and RX J1713.7-3946 (393 CE), 
intermediate age (w 10^ yr) SNRs, e.g., IC 443 and W51C, and middle-aged 
( > 2x lO-* yr) SNRs, including W44, W28, the Cygnus Loop, and G349.7+0.2. 
All SNRs detected with Fermi, other than Cas A, show evidence for molecular 
cloud interactions. 

Detecting the 7r° decay feature in Fcrmi-LAT SNR spectra depends on 
good background modeling to expose the low-energy ( ^ 200 MeV) flux. Shock 
compression of pre-existing cosmic rays could be argued to account for appar- 
ent cosmic-ray production at SNR shocks, but this seems less plausible now 
that we have measured GeV 7-ray luminosities ranging from w 10'^^ erg/s for 
the Cygnus loop, seen as an « 1° radial loop and the largest resolved 7-ray 
SNR, to w 10^4 - 10^5 erg/s for young SNRs, and « lO^^ - few x lO^^^ erg/s for 
SNRs showing strong molecular cloud interactions [235| . This implies cosmic- 
ray kinetic energy on the order of £'cit(erg) ~ 10"^^ Lae^pp = 10^^i36/n(cm~'^), 
in accord with the SNR hypothesis for cosmic-ray origin. This estimate is not, 
however, so strong, given that distance uncertainties can easily lead to factor 
of 4 or more uncertainty in L. 
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Fig. 18. Combined Fermi-LAT and TeV spectra of selected SNRs [236] 



For Cas A, one of the youngest SNRs detected with Fermi '237' , combined 
VHE (MAGIC, VERITAS, and HEGRA) data show a vF^' peak at a few 
GeV, with ambiguous evidence for a low-energy cutoff, but clear evidence for 
a TeV steepening (see Fig. 18). Cas A, with an extent of « 6', appears almost 
as a point source for the LAT. The somewhat older SNR RX J1713. 7-3946 
is famous for being the first SNR for which a VHE (HESS) map was made 
|238| . The TeV emission is extremely well-correlated with X-ray maps, e.g., 
HESS and Suzaku [239] . The Fermi-LAT map '240^ displays an asymmetry 
towards the north, in the direction of a molecular cloud where the X-ray and 
TeV emission is enhanced. The GeV spectrum itself is remarkably hard, with 
spectral index = 1.5 ± 0.1. The vFi, spectrum of this remnant peaks between 
« 200 GeV and 1 or 2 TeV, with a rapid cutoff at higher energies. Nevertheless, 
the spectrum of RX J1713. 7-3946 extends to ~ 100 TeV, making it the source 
of the most energetic 7 rays yet detected. Its hard spectrum cannot be made 
by conventional hadronic models. One-zone leptonic models can fit the X- 
ray/7-ray spectra of RX J1713. 7-3946, as can two-zone models. 

The prototypical Type la Tycho SNR has also been recently reported as a 
GeV source with a GeV- TeV 7-ray spectrum implying, if hadronic, a proton 
spectrum with a —2.3 number index |241j . Claims have been made that the 
combined LAT/ VERITAS spectrum requires hadrons |242j . but a two-zone 
leptonic model with both bremsstrahlung and Compton processes can fit the 
spectrum as well |243j . Because of the weakness of the GeV/TeV 7-ray flux, the 
vFi, peak of Tycho can be constrained to be ^ few GeV, and the spectrum is 
not well-resolved in the crucial energy range ^ 300 MeV where the tt" bump 
might be evident. For this purpose, a better SNR might be IC 443 [244] . a 3 - 
4 kyr old core collapse SN (like Cas A and RX J1713) surrounded by clouds 
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of gas. With an 0.75° radio radius, it is resolvable with Fermi. The SED of IC 
443 is a flat vF^, spectrum between w 300 MeV and 5 GeV with a peak at 1 
or 2 GeV. The spectrum extends to VERITAS/MAGIC VHE energies hke a 
power-law with i/Fiy slope equal to w —0.6. 

The middle-aged shell SNR W44 [21S] has an age of w 20 kyr. The 7-ray 
shape closely following the 4.5 /i Spitzer IR image which traces shocked H2. 
Its spectrum peaks at 1 GeV or less, and it falls off sharply at higher energies, 
so is not detected with the VHE telescopes. Its physical extent precludes the 
7-ray emission from originating from a pulsar, though a pulsar, B1853-I-01 
with age « 20 kyr, is found in this SNR and could be associated with W44. If 
so, W44 is probably the result of a core-collapse SNR, and therefore Type II. 

As a final example of Fermi-LAT observations of SNRs, consider W51C, 
the first resolved GeV SNR at GeV energies. It is an w 10 kyr old remnant 
with extended GeV emission compatible with the location and shape of the 
ROSAT emission [246] . The spectrum of W51C is similar to that of W44, 
peaking below « 1 or 2 GeV, and falling off steeply at higher energies. HESS 
weakly detects emission at about 1 TeV from W51C, (whereas W44 is not 
detected at VHE energies), implying a vF^ spectral index « —0.5 between 10 
GeV and 1 TeV. 

In terms of cosmic-ray origin theory, a few remarks can be made. First, 
almost all Fermi-LAT SNRs, other than Cas A, show evidence for interactions 
with molecular clouds, in particular, those with OH(1720 MHz) maser emis- 
sion from the OH hydroxyl molecule tracing dense, shocked gas [247| I248j . 
Out of 24 known maser SNRs reported in j247j . 10 have GeV and/or TeV 
associations, and 6 have both. Detection of illuminated molecular cloud com- 
plexes in front of the SNR shock could reveal the existence of runaway cosmic 
rays that more likely would be protons than electrons. 

Second, nonthermal cosmic-ray protons lose energy on the secondary nu- 
clear production timescale tpp = {nappc)~^ ~ 35 Myr/n(cm~'^), whereas 
nonthermal electrons lose energy through bremsstrahlung on the free-free 
timescale tff = (naffc)^^ « 6 x 35 Myr/n(cm^'^), where ct// ~ a/UT — 
(Tpp/6. In either case, the process is radiatively inefficient unless the target 
density n 3> 1 cm~'^, but the nonthermal electron bremsstrahlung is less ef- 
ficient by no more than an order of magnitude than the nuclear production 
cross section. The nonthermal bremsstrahlung or secondary nuclear produc- 
tion model face energetics problems if they requires ^ 10^" erg in cosmic 
rays. 

Third, ionization/ Coulomb losses can harden a nonthermal electron spec- 
trum, making 7 rays from a leptonic bremsstrahlung model masquerade as 
a tt" feature. A full spectral model for multiwavelength production from 
SNRs requires, most generally, particle acceleration at the forward and reverse 
shocks, zones of different magnetic field strength at which particle acceleration 
can occur, leptonic bremsstrahlung and Compton scattering, and secondary 
nuclear production. 
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Perhaps there is a trend in the SEDs of SNRs that can reveal the hkely 
7-ray production mechanism (Fig. 18). For leptonic models, both (i) nonther- 
mal bremsstrahlung Compton scattering, are potentially feasible. For hadronic 
models, only (iii) secondary nuclear/particle process is effective in the SNR 
environment. Confrontation of these three nonthermal processes with multi- 
wavelength SNR data imply power and spectral and morphological constraints 
that can in principle identify the dominant radiation process, though in prac- 
tice, this has not been so simple. But establishing a trend from hard to soft 
7-ray spectra with age, even scaled to the Sedov age, is oversimplified, as is 
evident from the rather soft spectrum of Tycho. 

For useful studies, a plot of the VHE (100 GeV - 10 TeV) flux divided 
by the GeV (300 MeV - 30 GeV) flux vs. SNR age (to avoid uncertainties 
in distance measurements) for Type la and Type II SNe might provide some 
insight on the evolution of particle acceleration with remnant age. At present, 
it still seems premature to claim that the problem of the origin of the Galactic 
cosmic rays is solved. 



5.5 Nonrelativistic shock acceleration of electrons 

First-order Fermi acceleration is highly developed as a mechanism to acceler- 
ate cosmic rays [2491 1250] . Here we illustrate this process by calculating the 
nonthermal spectrum of test particles accelerated at a shock-discontinuity in 
density and velocity. We use a continuity equation approach}^ which com- 
plements the approach to shock acceleration using the convection-diffusion 
equation and probability arguments 3 . 

The geometry we consider is particle acceleration at a discontinuity in ve- 
locity and density, as illustrated in Fig. |19| The upstream unshocked material 
has speed u in a frame where the supersonic flow intercepts gas at rest. For a 
strong nonrelativistic shock with compression ratio x = 4, the upstream (— ) 
flow approaches with speed w_ — /3_c = 4/3c/3 and the downstream (-f ) flow 
recedes with speed u+ — Pj^c ~ jic/Z in the comoving primed frame station- 
ary with respect to the shock. Consequently u = — = 13c. Here we treat 
the acceleration of relativistic nonthermal particles with Lorentz factor 7 S> 1 
and speed PparC « c. 

In the general case, the compression ratio 

X = - = — ■ (122) 
For a lossless adiabatic shock wave, 



^" 7-1 + 2/X? ^ = 1 + 3/X2 (123) 



The treatment of |251| assumes energy-independent escape timescale representing 
some second-order Fermi acceleration scenarios. 
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Fig. 19. Shock geometry in rest frame of shock. In this figure, S is the shock width. 

(eq. (13.8), [3]), where 7 is the adiabatic index, and 7 ~ 5/3 for a nonrela- 
tivistic monatomic gas. When AI1 S> 1 and 7 = 5/3, x ^ 4. It follows from 
eq. ( 122 ) and the condition u — u_ — u+ that 



x-i 



u+ = 



X-1 ' 



(124) 



so that in the limit Al 3> 1, u_ — > 4u/3 and u+ — > u/3 in the frame comoving 
with the shock. 

The change of energy of a relativistic particle with dimensionless momen- 
tum p = /3"f = \/"/^ — I — > 1, when 7 ^ 1, in a complete cycle of Fermi 
acceleration is, from eq. (12.10) 3J, given by 



I ^ ^Fl 3P ' 

where u — f5c = — m+ . For second-order Fermi acceleration. 



(125) 



(126) 



and both expresssions receive an added boost cx -T^, at least in the first cy- 
cle of acceleration, for shock acceleration by a relativistic flow with speed 

The rate of energy gain, or acceleration rate, in nonrelativistic first-order 
Fermi acceleration is given for particles 



7fi 



(127) 



where the cycle time tcyc is given by 
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tcyc = - ( 1 ) = (K- + 



=77^^ 0C7. (128) 



3u X V B+J QB^ 3u x 

|252j . where u ^ c is the particle speed, x = U-/u+ — p+/p^ is the compres- 
sion ratio, and p = mn is the mass density. The spatial diffusion coefficient 

K± = -X±v = -V±rL±v = -V±ri^±PV , 

here writing the diffusion coefficients in terms of the parameters r]± that give 
the particle mean-free-paths scaled to the values in the local magnetic field. In 
the Bohm diffusion approximation, the diffusion mean-free path is set equal 
to the Larmor radius, so 77^ = 1. Thus the Bohm approximation is 

1 



where the Larmor radius 



r-L± = r£±p = — - p . (129) 

The acceleration rate of nonthermal relativistic particles by first-order shock 
acceleration can therefore be expressed as 

4/? 4 « 
7acc = ^7 — 7 = o — 7 • (130) 



The acceleration rate is independent of p or 7 because tcyc c< 7 from eq. ( 128 1 . 
In nonrelativistic shock acceleration, the escape probability per cycle is 

Vcsc = ^13+ = = . (131) 

c c(x-l) 



|253| 1254) . The escape time is therefore 

, _ tcyc _ tcyc (^^^0\ 

~ Vcsc ~ 4/3+ ■ ^"'^ 

The steady-state particle continuity equation with gains, losses, and escape 
takes the form (eq. (CIO), [3]) 

Eq. ( |133 ) has solution, for an energy-gain process 7 > 0, given by (eq. (Cll), 
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Eqs. ( 130 ) and ( 132 ) imply 



tesc{rDiacc{i) = ^7 = ^^y-7 



so 



^3/(x-l) 
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(135) 
(136) 



noting that this represents the differential number density spectrum of particle 
accelerated at the shock discontinuity. 

Taking a (5-function source injection n(7) = rio(5(7' — 70), the injection 
spectrum downstream where the particles are no longer subject to acceleration 
at the shock is 
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where 
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is the test-particle number index well known in studies of nonrclativistic shock 
acceleration |250| . 

Now consider shock acceleration with the addition of Thomson/synchrotron 
radiative losses, so 
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Following the approach above gives 
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This is a pileup spectrum for x > 5/2, that is, when Atp is harder than 3. The 
maximum electron energy is given by 1 — k^"^ — 0, or 
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7™a. = ^ = . (143) 

Even though the spectrum piles up, the number of particles and the energy 
of the accelerated particles is convergent. 

The synchrotron emission made by this electron distribution has a maxi- 
mum value at dimensionless photon energy 

3 B 2 _ 97^e / X \ a2 {B-B+/B 



(^syn,max - ^ 7 - „ _ -. ) ^ / „ \ i3_ (1^4) 



(cf. [2SS]). Thus 



^syn^max 



(145) 



The apparent divergence due to the term % — 1 in the denominator of eq. ( 145 ) 
can be seen not to arise, considering that X"^/(x ~ 1) = U-U. The leading 
term, esyn,max — 27/8q;/ = 462, or Esyn,max — 236 MeV, represents a bound 
for nonrelativistic shock acceleration, remarkably close in value to the cutoff 
energy of the Crab pulsar wind nebula (SS]. This indicates that the pulsar 
wind nebula of the Crab is formed by a wind termination shock moving out 
at mildly relativistic velocities. Small changes in the bulk speed from a knot 
in the Crab pulsar wind, as imaged with Chandra, could produce the Fermi 
and AGILE flares from the Crab [7T]. Better imaging at w 100 MeV would 
localize the emission source, but is hardly possible with Fermi. 

Addition of a diffusion term will produce smoothed, realistic pile-up elec- 
tron injection distributions formed in first-order shock acceleration. The ad- 
dition of a diffusion term is under study in work with P. Becker. Note that 
these pileup functions differ from the Schlickeiser pile-ups |256l I257j where es- 
cape is always independent of particle energy. At relativistic energies and 
with relativistic flows, the maximum dimensionless synchrotron energy is 
esyn,max — (27/8a/)i^/(l + z), which is relevant when interpreting the max- 
imum photon energies in GRBs as a consequence of synchrotron emission 
formed by particle acceleration at an external shock. 



6 7 Rays from Star-Forming Galaxies and Clusters of 
Galaxies, and the Diffuse Extragalactic 7-Ray 
Background 

6.1 7 rays from star- forming galaxies 

Galaxies with ongoing star formation, most notably the Milky Way in which 
we live, are illuminated at 7-ray eneregies by secondary products of cosmic-ray 
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interactions with gas and dust. The characteristic 100 MeV 7-ray luminosi- 
ties of normal star-forming galaxies are ~ 10'^* - 10''" erg s~^, some factors 
of ~ 10^ - 10^" X smaller than those of active galaxies. Yet the star-forming 
galaxies vastly outnumber the AGNs. For example, the space density of a 
typical L« spiral galaxy like the Milky Way is w 3 x 10^ - lO'' Gpc~^, by 
comparison with the space density of FR II radio galaxies, which is « 2000 
Gpc~^ |258j . (The volume of 1 Gpc"^ extends to about z = 0.15 from the 
present epoch.) 

Besides the Milky Way, the Large Magellanic Cloud, detected earlier with 
EGRET with flux F_8 = 19, is now measured with Fermi LAT at the level of 
F_8 = 26 ± 2 [259]. At a distance of « 50 kpc, the LMC has about 10% of 
the mass and a supernova rate R:!20% of the Milky Way. The LAT resolves the 
LMC, and finds that 30 Doradus, its major star-forming region, is a bright 
source of 7 rays that does not consist of significant point source contributions. 
The 7-ray spectrum is consistent with an origin in cosmic-ray production. Its 
7-ray emission correlates well with massive star forming regions and an ionized 
H"*" template, but more poorly with neutral or molecular gas distribution. The 
7-ray emission morphology is surprisingly compact, and indicates that cosmic 
rays are accelerated in star-forming regions that are not very diffusive, thus 
accounting for the bright compact emission centered around 30 Doradus. 

Our other notable dwarf companion galaxy, the Small Magellanic Cloud, 
is also for the first time detected in 7 rays. The LAT measures a flux of 
F_8 = 3.7 ± 0.7 [260] from an extended, - 3° region. Unlike the LMC, the 
7-ray emission from the SMC is not clearly correlated with the distribution 
of massive stars or supernova remnants, though the emission may trace su- 
pergiant shells. 



Table 9. Properties of Star-Forming Galaxies" 



Galaxy 


d 


RsN 


Mcas 






Index 




(kpc) 


(century ^ ^ ) 


(lO^Mo) 




(10" ph/s) 


(10™ erg/s) 


MW 




2.0± 1.0 


6.5 ± 2.0 




11.8 ± 3.4'' 


1.2±0.5 2.2±0.15 


LMC 


52 ± 2 


0.5 ± 0.2 


0.67 ± 0.08 


26.3 ± 2.0 


0.78 ± 0.08 


0.041 ± 0.007 2.26 ± 0.11 


SMC 


61 ± 3 


« 0.12= 


fB 0.45 


3.7 ± 0.7 


0.16 ± 0.04 


0.008 ± 0.003 2.23 ± 0.12 


M31 


780 ± 30 


1.1 ± 0.2 


7.7 ± 2.3 


0.9 ± 0.2 


6.6 ± 1.4 


0.43 ± 0.09 2.1 ±0.22 


M82 


3600 ± 300 


20 ± 10 


2.5 ± 0.7 


1.6 ± 0.5 


250 ± 90 


13 ±5 2.2 ±0.2 


N253 


3900 ± 400 


20 ± 10 


2.5 ± 0.6 


0.6 ± 0.4 


110 ± 70 


7.2 ±4.7 1.95 ±0.4 



"Refs. [2611 [260l [2221 

'F.^ = 10-^F_8 ph (> 100 MeV)/cm^-s 

"^Fluxes and luminosities in 100 MeV — 5 GcV range 

'^Value is strongly dependent on assumed size of Galactic halo used to model MW; I222| find 

0.6 5 £-8 ~ 1-0 
"Ref. [262] 

Though not unexpected [263] , the discovery of starburst galaxies at GeV 
|261| and TeV energies is important to test predictions of cosmic-ray origin. 
The two nearest starburst galaxies, M82 and NGC 253, each at a distance of 
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« 4 Mpc, were detected at a level of F_8 = 1.6 and F-s — 0.6, respectively 
PST] . with the LAT, soon after being reported as VERITAS [IM] and HESS 
|265j sources, respectively. Properties of the brightest star-forming galaxies 
at 7-ray energies are shown in Table |9] Although correlations with mass, SN 
rate, and their product can be extracted, the situation, as exhibited by the 
distribution of 7-ray emission in the LMC, is far more complex. As can be 
seen from Fig. [Tj starburst galaxies define a separate track in the spectral 
index vs. luminosity plane. 

Now that the LAT sensitivity for star-forming galaxies has reached beyond 
our neighboring galaxies, we can expect the number to grow. The Andromeda 
galaxy, M31, has been detected at a level of F_8 = 0.9 [266] • Other starburst 
and infrared luminous galaxy detections from sources like Arp 220 |267j or 
Mrk 273 are keenly anticipated as the LAT photon statistics accumulate. 

6.2 7 rays from clusters of galaxies 

Clusters and superclusters of galaxies are the largest manifestations of the 
ongoing structure formation. Regions collapse by gravitation, driven primar- 
ily by dark matter fluctuations. The magnitude of the density fluctuation 
determines the formation time, insofar as larger structures form by accret- 
ing smaller clumps in a hierarchical merging scenario. The result, as seen in 
numerical simulation codes, is lumpy, continuous accretion forming filaments 
and webs of enhanced density and magnetic field. 

Clusters of galaxies are the most energetic events in the universe since the 
big bang, releasing as much as GM1M2/ R ^ 10^^ - 10^'' erg, here taking Mi ~ 
IQ^^Mq ^ IQM2 and R ^ 1 Mpc. Yet they are hardly the most luminous, 
since the energy is released over a Hubble timescale H^^ - 4.3 x 10^^ s - 13.6 
Gyr, implying a power of ~ 2 x 10"*^ - 2 x 10"*^ erg/s. This energy, carried 
primarily in the gravitational potential of the dark matter, goes into heating, 
turbulence, and gas motions. In the events of cluster mergers, turbulence 
is generated in wakes, and shocks are formed in the collision. Nonthermal 
particles are accelerated by shocks and turbulence. 

Evidence for nonthermal particle acceleration in merging clusters of galax- 
ies is already known from the existence of radio halos and relics. The relativis- 
tic electrons and positrons making this radio emission are either accelerated 
directly or produced as secondary pairs by cosmic-ray protons and ions collid- 
ing with particles in the intracluster medium. The unpolarized central radio 
halos, ~ Mpc size, with a morphology similar to the X-ray bremsstrahlung, 
reflect particle acceleration at the merger shocks between two merging clus- 
ters, intermittent AGN activity, and dynamical friction of galaxies in the hot, 
keV, thermal plasma. By contrast, radio relics, also about ~ 1 Mpc in size, 
which lie on the cluster outskirts, display elongated morphologies, and are up 
to ~ 50% polarized, may reflect a different origin from the radio halos. Most 
likely, they represent ongoing accretion of pristine gas from the big bang. 
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The halo radio emission of, e.g., the Coma cluster of galaxies, is very 
soft above ~ 1 GHz, and can be modeled in a merging cluster framework 
if the system is observed soon after the cluater and subcluster have merged 
[268j . By scaling nonthermal cosmic-ray proton and ion energy to the electron 
energy required to make the observed radiation, predictions for the 7-ray 
emission can be made. The identification of a nonthermal hard X-ray feature 
as Compton-scattered CMBR would mean a larger nonthermal electron total 
energy, giving more optimistic predictions for 7-ray emission from clusters of 
galaxies. In any case, clusters of galaxies act as storage volumes for the cosmic- 
ray protons |269[ . which have radiative lifetimes exceeding the Hubble time. 
Because of the long, ~ 10 Gyr, crossing time of electrons, versus a radiative 
lifetime of ~ 0.1 Myr, in-situ acceleration of pre-existing relativistic electrons 
by turbulence [2701 1271) is now favored to explain radio halos and relics. 

The Fermi LAT has many high-priority cluster candidates, including not 
only Coma and Abell 2256, but Hydra, Centaurus, and Fornax. The Perseus 
cluster is dominated by the AGN 3C 84 in the central elliptical galaxy NGC 
1275, and is detected with the Fermi LAT [272]. But this 7-ray emission is 
from a variable blazar core (though there may be steady emission at the level 
of F_g ^ 4). Radiation mechanisms involving the energetic leptons together 
with the decay of neutral pions produced by hadronic interactions have the 
potential to produce abundant GeV photons. Using data from 2008 August 
to 2010 February, upper hmits of 33 galaxy clusters, selected according to 
their proximity, mass, X-ray flux, temperature, and non-thermal activity were 
reported by the Fermi collaboration [273] . The flux upper limits, in the 0.2 - 
100 GeV range, are typically at the level of F_8 = 0.1-5. These results limit 
the relativistic-hadron-to-thermal energy density ratio to be ^ 5% - 10% in 
several clusters |2741I275) . 

As disappointing as it has been that 7 rays have not yet been detected 
from galaxy clusters, we are still learning about the injection conditions in the 
low Mach number merger shocks, the high Mach number accretion shocks, and 
the injection conditions of turbulence models. A purely secondary production 
model where all the radio emission results from secondary electrons formed in 
secondary nuclear production is likely ruled out |273[ I276j , and dark matter 
annihilation cross sections are further constrained. 

6.3 Extragalactic 7-ray background and populations 

Fig. [20] shows the spectrum of the diffuse extragalactic 7-ray background 
(EGB) obtained in an analysis of thefirst 10 months of Fermi-LAT science 
data in the range 200 MeV - 100 GeV based on the GALPROP model for the 
Galactic emission [20], alongside the EGRET EGB ^ and the EGRET EGB 
based on a GALPROP analysis of the Galactic diffuse EGRET emission [3T] . 
The Fermi-LAT intensity extrapolated to 100 MeV based on the power-law fit 
/(> 100 MeV) = (1.03±0.17) x lO^^/cm^-s-sr is significantly lower than that 
obtained from EGRET data, namely /egret(> 100 MeV) = (1.45 ± 0.05) x 
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Fig. 20. Fermi measurement of the diffuse extragalactic gamma-ray background 
(EGB) 20 (purple filled circles). Also plotted is the EGRET result [19 (filled black 
circles) and background intensity from GALPROP analysis of EGRET data [31] 
(open red circles). 



lO-Vcm^-s-sr [TO]. The GALPROP analysis of the EGRET data [HT] agrees, 
however, with the Fermi flux extrapolated to 100 MeV. Furthermore, the 
Fermi-LAT spectrum is compatible with a featureless power law with index 
-2.41 ±0.05 [2g. This is significantly softer than the EGRET spectrum, with 
index —2.13 ±0.03 To check that the different spectra are not due to the 
instrumental point-source sensitivities, a threshold flux Fig'' = 10, comparable 
to the average EGRET sensitivity, and an isotropic 7-ray intensity like that 
quoted above is measured. Therefore, the discrepancy cannot be attributed to 
a lower threshold for resolving point sources. Most likely the relatively hard 
spectral slope for the EGRET EGB was due to underestimation of the EGRET 
effective area above « few GeV [25^ , or due to lack of a model for Compton 7 
rays, which can make a significant contribution to the EGB at high latitudes. 

The origin of the EGB is an interesting open question. Because of the 
preponderance of blazars in the Third EGRET catalog [17], a large fraction 
of this emission was expected to be contributed by blazars, though some 
pre-Fermi studies [221 1277] found that blazars made only a small fraction 
of this emission. The detection of ~1000 high Galactic latitude sources [JD], 
most of which are blazars [84], has allowed a precise characterization of the 
contribution of blazars to the EGB. Such studies have shown that unresolved 
blazars contribute < 30% of the EGB emission [1751 UTS]. 

If not blazars, then other source classes must be invoked to account for 
the EGB. Possibilities include emission from star-forming galaxies, misaligned 
blazars (i.e., radio galaxies), 7 rays made by particles accelerated by structure- 
formation shocks in clusters of galaxies, |280| I281|, I282j . and dark matter; 
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see references in [278]. Star- forming galaxies such as our own Milky Way 
are known to be 7-ray emitters due to the interaction of cosmic rays with 
interstellar gas and radiation fields [261] . The superposition of the numerous 
but individually 7-ray weak star-forming galaxies can contribute a greater 
fraction of the EGB than the rare, individually bright 7-ray blazars [2831 1284j . 
Fermi detection of several star-forming galaxies, including NGC 253, M82, and 
the LMC |261j , supports this possibility, though the exact percentage remains 
highly uncertain. 

Included in the EGB are pulsar contributions, including millisecond pul- 
sars. Because millisecond pulsars (MSPs) are also found in early-type galaxies, 
they might contribute a significant fraction of the EGB. A way to determine 
their integrated diffuse emission would be through the integration of their flux 
(or luminosity) distribution in 7-rays. The typical spectrum of MSPs is hard 
(e.g. photon index of ^1.5) and shows an exponential cut-off around a few 
GeV [ST]. Thus the diffuse emission arising from MSPs should show a bump 
around a few GeV, similar to a feature found in the Galactic diffuse emission 
that has been ascribed to annihilating dark matter |285j . The left-hand panel 
of Fig. [21] illustrates a possible decomposition of the EGB into different source 
classes, including a contribution from WIMP dark-matter annihilation. 




10^ 10' 10" 10' " 1 2 2-5 3 3.5 



Energy [MeV] ^ 

Fig. 21. Left panel: Fermi-LAT measurement of the extragalactic diffuse back- 
ground spectrum (data points; see [20]), and model contributions from different 
source classes. The curves represent the integrated contributions of: 1) FSRQs (thin 
solid line), 2) BL Lac objects (thin short dashed line), 3) star-forming galaxies, 4) 
MSPs (dot- long dashed), 5) WIMP bb annihilation (long dashed) and 6) sum of all 
the previous (hatched). Right panel: Luminosity function of 7-ray selected FSRQs 
[575] . 



A preliminary 7-ray luminosity function for FSRQs detected with the 
Fermi-LAT is shown in the right-hand panel of Fig. 21 ^275] . The luminosity 
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function shows the change in the space density of FSRQs as a function of 
redshift and for different luminosity classes. What is apparent from Fig. [21] is 
that the space density of powerful FSRQs peaked at redshifts ~ 2, and then 
declined to the space density now observed. The peak of maximum growth 
occurs at different epochs for different luminosity classes, with the more lumi- 
nous sources reaching their maximum space density earlier in the history of 
the universe, while the bulk of the population (the lower luminosity blazars) 
are more abundant at present times. This cosmic downsizing behavior cannot 
persist to very large {z ^ 1) redshifts and early times, as may be reflected by 
the sudden turnoff of high-z 7-ray blazar. The largest redshift 7-ray blazar 
in both the ILAC and 2LAC is at z = 3.1, despite many candidate higher- 
redshift blazars. 



7 Microquasars, Radio Galaxies, and the EBL 
7.1 7-ray binaries 

The literature on the subject of Galactic binaries is vast; a good X-ray review 
remains Ref. [286) . Here we highlight only a few of those aspects that make 
for an appreciation of GeV - TeV 7-ray binaries in the Galaxy, of which 6 are 
well-established, namely 

1. LS 5039 (P = 3.9 d); 

2. LSI +61° 303(P = 26 d); 

3. PSR B1259-63 (P = 3.5 yr); 

4. Cygnus X-3 (P = 4.8 hr); 

5. IFGL J1018.6-5856 (P = 16.6 d) [287]; and 

6. Cygnus X-1 (P = 5.6 d). 

These sources are exclusively identified with high-mass X-ray binary (HMXB) 
systems and are identified, not merely associated, with the named system by 
period folding or time differencing with the orbital period. With its long or- 
bital period, this is harder for PSR B1259-53, but 7-ray flaring near periastron 
strengthens the identiflcation. PSR B1259-63 is an example of a Be X-ray bi- 
nary system, where the compact object, a young pulsar with a spin period 
of 48 s, is in a highly eccentric orbit around a luminous B star that displays 
emission lines in its optical spectrum (thus the designation "Be" ) . Flares take 
place near periastron at all frequencies, when the pulsar passes through the 
equatorial stellar wind [288j , and now also at GeV energies [289j . 

High mass and low mass refer to the companion star, not the compact ob- 
ject. HMXBs are X-ray binaries composed of a compact object — white dwarf, 
neutron star, or black hole — and a young massive O or B main-sequence or 
Wolf-Rayet star. The low-mass X-ray binaries (LMXBs) typically have older, 
less massive companions ^ 2Mq. If the companion compact object is a pul- 
sar, it can vary in age from young and energetic to old and slowed down 
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(and/or up). The fuel illuminating galactic X-ray binaries is accretion onto 
the compact object, which is transferred, primarily, through stellar wind in 
HMXBs and Roche-lobe overflow in LMXBs. X-ray variability is modulated 
by orbital effects, the most important being the change in the binary orbital 
separation distance. In only the case of PSR B1259-63 are X-ray pulsations 
at the pulsar period seen, which are also found for intermediate mass X-ray 
binaries like Her X-1 and other HMXBs at X-ray energies. 

Microquasars are X-ray binaries with radio jets. HMXBs number around 
7, including SS 433 and V4641 Sgr to the list above. There are 9 LMXB 
microquasars in the Paredes list from 2005 [290| . Some of the most famous 
are GX 339-4, Circinus X-1 and Scorpius X-1. 

Like the radio structures seen in extragalactic AGNs (e.g., 3C 84), jetted 
structures have been found in unusual X-ray binary sources, e.g., the Hertz- 
Grindlay source IE 1740.7-2942 lying ^ 0.7° away from the Galactic center 
|291j . The analogy with extragalactic radio sources was strengthened when 
this source and other LMXB sources such as GRS 1915-f 107 [253] were 
shown to exhibit superluminal effects. No LMXBs are yet detected at GeV or 
TeV energies. 

Prior to Fermi, three Galactic binaries were detected at TeV energies, 
namely LSI + 61° 303 with MAGIC and VERITAS, and LS 5039 and PSR 
B1259-63 with HESS [lH]. An isolated flare from Cyg X-1 was observed with 
MAGIC. No solid EGRET identification with a Galactic binary source had 
been made, though both LS 5039 [295] and LSI + 61° 303, even in the COS-B 
days [H], were suspected to be 7-ray sources. 

A central debate regarding 7-ray binaries is whether they are scaled-down 
quasars, as suggested by the microquasar label, or colliding stellar wind/pulsar 
wind systems. Leptonic microquasar models [2961 12971 1298j invoke electrons 
or hadrons in jets accelerated up to TeV energies with mildly relativistic 
outflows, accounting for the extended radio emission. Confirming evidence 
would be VHE emission from X-ray binaries with definite black holes, e.g., 
Cyg X-1, V 4641, or GRS 1915-1-105. 

The other class of pulsar/star model invokes the rotational energy of the 
neutron star and the interactions of the pulsar and stellar winds. [299|, 1300] . 
This model surely apphes to PSR B1259-63. Insofar as LS 5039 and LSI 
-|-61° 303 have compact objects with M < AMq, the question of whether the 
compact object is a black hole or neutron star is unresolved. The confirming 
evidence for a neutron star is, of course, detection of pulsations. 

When considering the different sources, it is useful to have a theoretical 
picture in mind. Most developed and perhaps most likely is a leptonic model, 
where electrons, which are accelerated in the extreme environment near the 
compact object, Compton scatter the photons of the hot star. Those 7 rays are 
in turn attenuated by the stellar radiation. To first order, a source at superior 
conjunction should be GeV luminous (more head-on collisions for Compton 
scattering) and less TeV luminous (due to 77 attenuation). It should anti- 
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correlate with the TeV emission, which has least 77 opacity and therefore 
should be brightest at inferior conjunction. 

LS 5039 

Though sitting on the edge of the Galactic plane, LS 5039 is well-detected 
and identified by phase- folding |301| . Its Fermi-LAT light curve is remarkably 
anti-correlated with the HESS TeV emission. At superior conjunction, when 
the TeV flux is in a low power-law state, the Fermi GeV spectrum is in a high 
state. At inferior conjunction, the roles are reversed. Only quasi-simultaneity is 
possible, because Fermi will necessarily integrate over longer times to achieve 
the same sensitivity as a TeV telescope. The joint spectra are also remarkable, 
with the GeV high state fit as an exponentially cut-off power law, and taking 
place at the same phase of the low-state power-law HESS spectrum. In reverse, 
the low-state GeV power law takes place when the TeV flux is elevated. 

LSI +61° 303 

The Fermi LAT light curve of this source is a rather sharp-peaked sinusoid, 
with peaks of emission around periastron |302j . The GeV radiation is modu- 
lated with a maximum close to periastron, and minimum close to apastron. 
TeV emission peaks at different orbital phase, revealing a flux anti-correlation 
with the MAGIC TeV emission that might be expected in a leptonic model. 
The phase-averaged Fermi spectrum is well described by an exponentially cut- 
off power-law with index « 2.2 and cutoff energy « 6.3 GeV. This spectrum 
is not found to vary strongly with phase. The non-simultaneous MAGIC and 
VERITAS data corresponding to different phases look, if they correspond to 
the time of Fermi observations, like a distinct spectral component, with the 
hard component extending to TeV energies with a spectra index a; 1.9 over- 
taking the exponentially decaying GeV spectrum. 

Cygnus X-3 

Cyg X-3 is a HMXB microquasar with an orbital period of 4.8 hours — short 
for an HMXB — at a distance of ~ 7 kpc from the Earth. It is found in a 
complicated region, the Cygnus arm, but its period readily identifies it in 7 
rays. Its companion is the Wolf-Rayet star V1521 Cyg. It exhibits bright radio 
flares and goes through active states. It is not known whether it contains a 
black hole or neutron star binary. Four episodes of flares were observed with 
AGILE [303] a few days before radio/X-ray flares. The Fermi LAT light curve 
exhibits periods of scattered activity, and periods of quiescence that appear to 
be anti-correlated with the X-ray emission |304j . The radio and 7-ray activity 
are also correlated. 
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IFGL J1018.6-5856 

Period-folding analysis has resulted in the recent discovery of a new high-mass 
X-ray /7-ray binary, IFGL J1018.6-5856 [1HZIES5]. The 100 MeV - 200 GeV 
Fermi-LAT data is modulated with the 16.58 ± 0.04 d orbital period. The 
Swift XRT finds an X-ray source coincident with the position of the Fermi 
source. The digitized sky survey (DSS2) image of this region shows a star 
with spectrum very similar to that of the gamma-ray binary LS 5039, and the 
Australia Telescope Compact Array detects a radio source at 9 GHz. Chandra 
and XMM-Newton observations [306] find a coincident X-ray source from 
which the column density and reddening can be determined. With 2MASS 
optical observations, this implies a star of type of 06V at a distance of 6 - 12 
kpc, large by comparison with the distances of LS 5039 and LSL 

PSR B1259-63 

PSR B1259-63 may help us understand these systems best by virtue of the 
compact object being unambiguously a neutron star. The pulsar has a 47.76 
ms rotation period and a massive, « IOMq Be star companion, SS 2883. 
The stars orbit with a 3.5 yr orbital period and eccentricity = 0.87, guaran- 
teeing episodes of strong interactions near periastron passage. Because of its 
inclination and the geometry of the equatorial stellar wind, the pre- and post- 
periastron passages are not symmetrical. EGRET on CGRO failed to detect 
it during 7 weeks of observation around the 1994 periastron passage |307j , but 
flaring at radio, optical, X-ray, and HESS TeV j294j has been monitored in 
previous passages. 

So it was with great anticipation that observers at all wavelengths waited 
for the 2010 December passage. PSR B1259-63 did not fail to perform: the 
light curves measured by an array of instruments at different wavelengths 
and epochs are shown in Fig. [22] and reveal the complexity of the encounter 
|289j . It was only weakly detected with the Fermi LAT on the approach to 
periastron and up to 20 d after (MJD 55517-55563), during which the pre- and 
near-periastron spectrum in Fig. [22] was measured. The system flared at GeV 
energies starting at w 30 d post-periastron, with the Fermi LAT measuring 
during period MJD 55575 - 55624 shown as the post-periastron spectrum in 
Fig. [22] The reason that the GeV flare seems to lag behind the peak of the 
radio and X-ray emissions at « 20 d post-periastron is not understood. 

Though the time of 2010 periastron passage was not ideal for TeV ob- 
servations due to moonlight conditions, additional post-periastron TeV data 
for the most recent passage will help illuminate this system. We can reliably 
predict interesting events in the future. 

7.2 Misaligned blazars and radio galaxies 

The misaligned AGN (MAGN) population detected with Fermi using the first 
15 months of science survey data consists of 11 sources (Fig. [7] Hi]), includ- 
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Days from periastron 



Fig. 22. Left: Radio, X-ray, GeV, and TeV light curves of PSR B1259-63 at different 
epochs around the 2011 perastron passage on 2010 December 15 (MJD 55545), as 
labeled [289]. The 120 days shown represents ~ 13% of the orbital period. Right: 
Fermi spectra of the system in the pre- and near-periastron period (blue circles) and 
in the post-periastron flare (red squares). 



ing 7 FRI and 2 FRII radio galaxies, which are thought to be the parent 
population of BL Lac objects and FSRQs, respectively, in addition to 2 steep 
spectrum radio quasars that are believed to be slightly misaligned FSRQs. The 
MAGNs are associated with objects in the Third Cambridge (Revised) and 
Molonglo radio catalogs. The MAGN sources in the 3CRR catalog have large 
core dominance parameters compared to the general 3CRR source population, 
implying that the beamed component makes an appreciable contribution to 
the 7-ray flux. This is furthermore supported by the fact that three of these 
sources — 3C 78, 3C 111, and 3C 120 — do not appear in 2LAC, evidently due 
to variability. 

Why the ratio of measured 7-ray luminosities of FRI galaxies and BL Lac 
objects span a much larger range than the comparable ratio for FRII radio 
galaxies and FSRQs is an open question, if in fact not due to statistics. The 
simplest possibility is that there are a lot more nearby FRIs, and limitations 
on Fermi sensitivity will therefore favor detection of these nearby sources. At 
z <C 1, the luminosity distance « 4200z(l -t- 2) Mpc. As noted earlier, the 
Fermi-LAT reaches a limiting energy-flux sensitivity of w 5 x 10"""^^ erg cm~^ 
s~^ for two years of observations that is, unlike integral photon flux, only 



Sources of GeV Photons and the Fermi Resuhs 



91 



weakly dependent on source spectral index [84j- Fermi can thus only detect 
sources with luminosity ^ 10'^^z^{l + 2z) erg s~^. 

The number of nearby FRI and FRII galaxies depends on the space density 
of these objects which, in the study by Gendre, Best, and Wall .258 , is inferred 
from 1.4 GHz NVSS-FIRST radio observations (recall that the Fanaroff-Riley 
dichotomy is based on 178 MHz luminosities Pns mhz)- From their Figure 12, 
FRI and FRII radio galaxies have a local {z < 0.3) space density of 35, 000 
Gpc-3 and « 2,000 Gpc~^, respectively, at F1.4 qhz > 10^^ W Hz"! sr'^ or 
Lr ^ 10^^° erg/s, implying that the local space density of FRIs exceeds that 
of FRIIs by « 20. 

The volume of the universe within z < 0.15 (c^l < 700 Mpc) is «1 Gpc"'^, 
so there are thousands of FRIs and FRIIs within z = 0.1, yet other than 3C 
111, all LAT-detected radio galaxies within this volume are FRI sources (see 



Figure 29 |308j ). Sensitivity limitations and the abundance of nearby FRIs 
could explain the MAGN population statistics, but low apparent-luminosity 
off-axis FRIIs, which are far more radio luminous than BL Lac objects and 
FRI galaxies, are lacking. With only a few thousand randomly aligned sources 
within z = 0.1, a narrower 7-ray beaming cone in FSRQs, with a more rapid 
fall-off in off-axis flux, makes detection of these nearby sources far less likely 
than a broader 7-ray emission cone in BL Lac objects, as expected for the 
different beaming factors from SSC emission and external Compton processes 



[3091 1310j . This could also reflect differences in jet structure between FSRQs 
and BL Lac objects |31H 1312] . or extended jet or lobe emission in FRIs |313j 
that is missing in FRII galaxies. 

Both core and lobes can significantly contribute to the measured 7-ray 
luminosities. In the case of Centaurus A, the values of of the core and lobes 
are comparable, with the lobe emission primarily attributed to Compton- 
scattered CMBR [41]. The significant or dominant lobe component means 
that the core luminosity of misaligned AGNs can be less than the measured 

unless the lobe emission is resolved. 

One-zone synchrotron/SSC models with Sb ;i 10 give good fits to the long- 
term average spectra of other HSP BL Lac objects such as Mrk 421 and Mrk 
501. The SEDs of radio galaxies, in contrast, are fit with much lower Doppler 
factors. The SED of the core of Gen A, for instance, can be fit with « 1 
and bulk Lorentz factors F « few |314) . Likewise, the SEDs of the FRI radio 
galaxies M87 and NGC 1275 [172] and M87 [3H are well fit with Sb ^ 2 and 
r - 4. 

The much larger values of F for BL Lac objects than for their putative 
parent population, the misaligned FRI radio galaxies, is contrary to simple 
unification expectations. Moreover, even though the 7-ray luminosities from 
FRI radio galaxies are much smaller than that of BL Lac objects (Fig.jTj), they 
are still larger than expected by debeaming the radiation of BL Lac objects 
with F 20. Additional soft target photons that can be Compton scattered 
to high energies result in a reduction of the value of (5d compared to those 
implied by the one-zone synchrotron/SSG model. These target photons can be 
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produced in a structured jet, as in the spine and sheath model |316j . Another 
soft photon source arises if blazar flows decelerate from the inner jet to the pc 
scale |317j , seemingly in accord with the mildly relativistic flows at the sub-pc 
scale found in radio observations of Mrk 421 and Mrk 501. 

7.3 The EBL 

In the most general sense, the EBL refers to the cosmic microwave background 
radiation (CMBR), the radiations from all the past stars and black holes, the 
glow from annihilating and decaying dark matter, and any residual emissions 
from exotic particle decays. The EBL at infrared, optical and UV frequencies 
originates predominantly from stellar emissions, cither directly or after being 
absorbed and reradiated by dust. The present EBL is a consequence of the star 
formation history of the universe, light absorption and re-emission by dust, 
and different types of dust extinction in various classes of galaxies throughout 
cosmic time. The IR and optical/UV EBL attenuates TeV and 10 - 100 GeV 
7-rays, respectively. Knowledge of the absorption due to the EBL is therefore 
necessary to infer the intrinsic, escaping spectra of extragalactic 7-ray sources. 

Directly measuring the EBL at IR and optical wavelengths is difficult be- 
cause of the interfering foreground zodiacal light and Galactic synchrotron 
radiation. Many attempts have been taken to predict the EBL intensity. Em- 
pirical methods |318l 13191 1320j sum optical/IR emissions from galaxy-counts 
using, e.g., luminosity-dependent galaxy SEDs, and extrapolating to high red- 
shift where data is lacking. Semi-analytic models of the EBL [3211 13221 1323) 
model galaxy formation following mergers of dark matter halos, with effects 
of supernova feedback, dust attenuation, and metal production considered. 
Other models [324J 13251 1174) are based on integrating stellar light with dust 
absorption, employing arguments for the star formation rate and the initial 
mass function of the stars. A lower bound to the EBL intensity can be de- 
termined by galaxy counts at different frequencies and redshifts. The upper 
panel in Fig. |23| shows a comparison of observations with models, and inferred 
constraints. Curves are model C from [174] (solid black curve); best fit model 
from |325j (red dotted curve); j320l (dashed green curve); the fiducial model 
of |322) (dashed blue curve); and the fast evolution and baseline models from 
|319) (upper and lower dot-dashed violet curves, respectively). Orange curves 
and inverted triangles are upper limits derived from blazar observations by 
|326) and [327J, respectively. 

Gamma-ray astronomy offers a special technique to probe the EBL at IR 
through UV frequencies, because photons of the EBL attenuate 7 rays via 
pair production through the 77 — > e"'"e~ process. The threshold condition 
implies that 400 GeV photons can pair produce on soft EBL photons with 1 
wavelength, and 4 TeV photons with 10/i EBL photons, etc., so thatV 
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Fig. 23. Upper: EBL observations, models, and constraints |174| . See text for refer- 
ences and details. Lower: Highest energy photons observed from blazars and GRBs 
with the Fermi-LAT as a function of the source redshift. Curves showing the t^^ — 3 
optical depth for EBL attenuation are shown for different models as labeled [328) . 
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Analysis of the effects of this process have already ruled out high-intensity 
models of the EBL at IR/optical frequencies [328] . Ruling out other EBL 
models will require longer integration times to improve implementation of 
tests like the following: 

1. Flux ratio method. By assuming that intrinsic blazar spectra are indepen- 
dent of redshift, then the ratio of high to low energy fluxes will decrease 
with redshift due to enhanced EBL attenuation |329j . It is necessary to 
use separate blazar populations in this test since Fermi has established 
significant spectral differences between FSRQs and BL Lacs, but this test 
is not yet sensitive enough to identify any trend of spectral softenings due 
to EBL absorption. 

2. Deabsorbed ^-ray spectrum. By extrapolating the Fermi-LAT spectrum of 
an individual source to high energies, and assuming that the spectrum 
does not harden with energy, then deabsorbing the spectrum limits the 
range of EBL models that are compatible with this underlying assumption 
[330]. Upper limits are also placed on the allowed EBL intensity by as- 
suming that the intrinsic source spectrum is limited in hardness by some 
theoretically determined value 331, 326], which depends on details of the 
particle acceleration mechanism [332J. Specific approaches determine the 
likelihood of different EBL models to be compatible with the measured 
7-ray spectrum by assuming an intrinsic (usually power-law) source spec- 
trum. 

3. Distribution of highest energy photons. The probability of detecting high- 
energy photons depends on the 77 optical depth associated with a given 
EBL model. The detection of such photons, though rare in number, can 
place severe demands on the apparent power budget of a source. If the 
absolute power or energy if separately determined, e.g., by using optical 
breaks to infer the jet opening angle in GRBs, then the required energy 
release may be extreme and call the EBL model into question (see the 



lower panel in Fig. 23 ) 



Detection of lobes of radio galaxies. Photons of the EBL are Compton 
scattered to GeV energies by relativistic electrons in the radio lobes of 
nearby extended radio galaxies. Thus measurements of the GeV lobe flux 
can be used to determine the level of the EBL [333] . The first radio galaxy 
with radio lobes resolved by Fermi is Cen A [IT]. Due to the low signif- 
icance of the Fermi detection and the lack of data at energies ~2 GeV, 
the method cannot yet be applied. Longer exposures should reduce the 
error bars below 2 GeV and detect photons in the critical regime above 
~2 GeV where the emission from the up-scattered EBL is expected to 
dominate. Besides Gen A, Cen B and Fornax A — two other radio galax- 
ies with lobes — have recently been detected [51]. These objects, as well as 
NGC 6251 |334j . are potential targets for resolving extended lobe emission 
to measure or constrain the EBL. 
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Implicit in the above methods is the assumption that the 7 rays are made 
at the source. This need not be the case if blazars or GRBs are sources of 
UHECRs. UHECR protons that are able to escape from the blazar and struc- 
tured regions without deflection can deposit energy in transit through the IGM 
via photopair production. The cascade radiation can produce persistent TeV 
emission that would confuse the interpretation of the attenuation due to the 
EBL. This is less likely for UHECRs accelerated by GRBs, where deflection 
of the UHECRs would suppress the emission generated by the UHECRs. This 
is currently an active area of research; see, e.g., [3351 13361 13371 133S1 13391 IMP] . 



8 Fermi Observations of Gamma Ray Bursts 

The effective lifetime for GRB studies using EGRET's spark chamber on the 
Compton Gamma-ray Observatory ended ~ 4.5 yrs into mission, after 1996 
[IS] p"] The depletion of spark-chamber gas was mitigated through the intro- 
duction of a narrow-field mode suitable for pointed observations. This made 
the chance of catching a GRB, proportional to EGRET's FoV, too improbable 
without rapid automated slewing, which was not possible for CGRO. Conse- 
quently EGRET only detected a total of five spark-chamber GRBs, all early in 
the mission [533]. These are GRB 910503, GRB 910601, the superbowl GRB 
930131, the famous long-lived GRB 940217 011], and GRB 940301. In the 
wide-field mode, EGRET was sensitive to « 1/25*'' of the full sky, which is 
« 1/5*'' as large as the FoV of Fermi 521]. 

8.1 Fermi LAT observations of GRBs 

Since Fermi operations began, 13 GRBs were reported as significantly detected 
in the LAT by the Fermi Collaboration from early August 2008 through cal- 
endar year 20090 The year 2010 saw a dearth of bright LAT GRBs, but 
the rate has picked up in 2011 with the detection of a few remarkable events 
like GRB 110721A and GRB 110731A, the former of which displays classic 
"fast-rise, exponential decay"-type GBM light curves with prompt LAT emis- 
sion [3451 1346] ■ and the latter of which is the first long-duration GRB jointly 
detected with Fermi-LAT and Swift from the prompt phase into the afterglow 
|3471 1348] . In the meantime, the development of the LLE technique ([2S]; Sec- 
tion 1.2) has found LAT emission from several GRBs in the 30 MeV - 100 
MeV range. Given the addition of these GRBs, the rate of GRB detection 



^° Very bright GRBs like GRB 990123 [MB could still be detected far off the COMP- 
TEL and OSSE axes while making a signal in the Energetic Gamma Ray Experi- 
ment Telescope's Total Absorption Shower Counter. EGRET TASC and BATSE 
data were used to made the discovery of the additional hard component in GRB 
941017 [342] . 

See fermi.gsfc.nasa.gov/ssc/observations/types/grbs 
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with the LAT is about 1 per month, and the number of Fermi-LAT GRBs has 
reached nearly 30. Because analysis is in progress on 2011 GRBs, and there 
were so few GRBs in 2010, here we review GRB observations only through 
2009. 

All LAT GRBs during this period are also GBM GRBs and comprise 
the brighter GBM GRBs, as already expected from a comparison between 
EGRET and BATSE GRBs in terms of fluence [349] • The 13 Fermi LAT GRBs 
observed through 2009 include 11 long GRBs and 2 short bursts, namely GRB 
081024B and GRB 090510 {z ^ 0.903). The most studied— because they are 
brightest— GRBs are GRB 090902B (z = 1.822) [350 , which provides the first 
strong evidence for a hard spectral component in long GRBs; GRB 080916C 
{z = 4.35) [138], the first bright long GRB; and GRB 090926A (z = 2.106), a 
burst with a narrow spike from the lowest to highest energies in an SED that 
requires both a Band function and a hard power-law component to fit. GRB 
090926A also reveals an extraordinary spectral softening at ^ 1 GeV in its 
time-integrated spectrum when the hard LAT spectral component is bright. 
These long GRBs had GBM fluences J" in the 20 - 2000 keV range > 10"'' erg 
cm-2 (Fig. [24]). The bright, short GRB 090510, with J" ^ 10"^ erg cm'^, also 
shows (like GRB 090902B and GRB 090926A) a distinct high-energy power- 
law spectral component in addition to a Band component |129j . Its short 
duration, large distance, and the detection of a 31 GeV photon permit strict 
tests on quantum gravity theories that predict a dependence of the speed of 
light in vacuo that is linear with energy |351j . 

Besides these notable GRBs are the less well-known and also less fluent 
long duration GRBs 090323 (z = 3.57), 090626 and 090328 {z = 0.736), the 
widely off-LAT-axis GRB 081215A, the first LAT GRB 080825C [352], the 
unusual GRB 090217 (353j showing none of the typical properties of LAT 
GRBs, and the unremarkable LAT GRBs 091003A and GRB 091031. The 
weakest fluence GRB of the sample is the first short GRB detected at LAT 
energies, GRB 081024B t354j, with J" « 4 x lO""^ erg cm'^. The weakness of 
this GRB could be related to the high Epk « 2 - 3 MeV of its Band-function 
component, but the time-averaged Epk = 4 MeV for GRB 090510 between 
0.5 and 1 s after trigger is even higher |129j . 

For those GBM GRBs occurring within the LAT FoV, detection of GRBs 
with the LAT is almost guaranteed when T ^ 10~* erg cm~^. The detec- 
tion rate slips to less than 50% when w 3 x 10~^ erg cm~^, and becomes 
highly improbable for T <, 10~^ erg cm^-^. This behavior undoubtedly reflects 
a distribution in the ratios of > 100 MeV LAT to GBM energy fluence [355] . 

In the first 16 months of Fermi science operations, ^ 1 GRB per month 
was detected with the Fermi LAT, or « 9 GRB/year, with LAT detecting 
short GRBs at « 10 - 20% of the rate of long GRBs. GRBs are detected 
with the GBM at a rate of 250 GRB/yr, or « 500 GRB/yr (full-sky). When 
corrected for FoV, EGRET detected « 25 GRB/year (full sky), while the 
Fermi LAT detects sa 50 GRB/yr (full sky). Given the much larger effective 
area of Fermi than EGRET, by a factor w 6 (8000 - 9000 cm2)/(1200 - 
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1500 cm^)], this small rate increase is something of a surprise, compounded by 
the ongoing sparse period of Fermi LAT detections of GRBs in 2010. Part of 
this difference is the stronger detection criteria of Fermi LAT than EGRET. 
But an improvement in flux sensitivity by a factor « 6, with an accompanying 
rate increase by only a factor 2 - 3 suggests that LAT GRBs are being 
sampled in a portion of their logA^ — log distibution that is flattened by 
cosmological effects. This is consistent with the known redshifts of LAT GRBs, 
which range from « 0.7 to z = 4.35, with a very rough average redshift of 
(z) = 2 for long GRBs and (z) « 1 for short GRBs (based only on GRB 
090510). If typical, both classes of GeV-emitting GRBs would be subject to 
strong cosmological effects on the fluence and flux distributions. 
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Fig. 24. Fluence-fluence diagram showing 6 Fermi |354j GRBs (red data points) 
and 5 EGRET spark-chamber [31H] GRBs (black data points). The EGRET fluence 
is measured from 100 MeV to 5 GeV, whereas the Fermi LAT fluence is measured 
from 100 MeV to 10 GeV. Short hard GRBs are circled. 



Fluence-fluence diagram 



Figure 24 shows the fluence-fluence diagram for the 5 EGRET spark-chamber 
[349] and for 6 Fermi ^354J GRBs. Most GBM GRBs have J" < 10^^ erg cm'^, 
and are only rarely detected with the LAT. Because of the small number of 
LAT GRBs, it is not yet clear whether there is a systematic difference between 
fluence ratios of EGRET and Fermi LAT-detected GRBs. The weakest Fermi 
LAT GRBs in terms of GBM fluence are both short duration GRBs. This 
could indicate a preference for short GRBs to have a larger ratio of LAT to 
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GBM fluences than long GRBs, depending on possible triggering biases, e.g., 
increased LAT background for long GRBs. 

The apparent isotropic energies Eiso of GBM and LAT emission for LAT 
GRBs with known redshifts are in several cases ^ 10^'' erg. For GRB 080916C, 
£iso ~ 10^^ erg. The LAT GRBs tend to have the largest energies of all 
measured GRBs, and as a result are good for radio studies [153j . 

Fermi LAT GRB phenomenology 

Some distinct behaviors have been identified in Fermi GRBs, namely: 

• Extended (long-lived) LAT (100 MeV - GeV) emission compared to the 
GBM (20 keV - 2 MeV) emission, known already from EGRET observa- 
tions, especially GRB 940217 [Mi] . 

• Delayed onset of the LAT emission compared to the GBM emission, in 
both long and short classes of GRBs. 

• Power-law temporal decay profiles of the LAT extended emission, decaying 
typically cx t''^-^ [356] . 

• Appearance of separate power-law spectral components with photon num- 
ber index harder than —2. 

• Delayed onset of the lowest energy GBM emission at ss 10 keV, seen for 
example in GRB 090902B and GRB 090926A. 

• Quasi-thermal Band function components with steep Band (3 found, e.g., 
in GRB 090902B at £; > 1 MeV [?50] . 

The onsets of the > 100 MeV emission appear to be delayed by ~ O.ltgo 
compared to the 100 keV - MeV emission (with tgo measured, e.g., in the 50 - 
300 keV GBM/BATSE range). This is one of the key and unanticipated results 
on GRBs from Fermi, and it appears to operate equally for both the long- 
and short-duration LAT GRBs. There have as yet been no LAT detections 
of members of the low-luminosity/sub-energetic class of GRBs that includes 
GRB 980425 and GRB 030329, nor have any X-ray flashes or X-ray luminous 
GRBs been detected with the LAT. Because GBM's primary triggering modes 
are similar to BATSE, high Epk^ relatively low-z GRBs are more likely to be 
detected compared to Swift. 

The luminosity function can be constructed from observations of GRBs. 
This function gives information about the rate density, the opening angle, and 
the total nonthermal energy radiated by GRBs. A description of the luminos- 
ity function for GRBs, and how the luminosity density can be evaluated for 
different GRB models is now given. 

8.2 GRB luminosity function 

The luminosity function Y{L) describes the distribution of apparent isotropic 
luminosities of members of a common source class. For GRBs, one can con- 
struct the luminosity function for, at least, the long duration (classical) GRBs, 
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the short hard GRBs, and the X-ray flashes and sub-energetic (low- luminosity) 
GRBs. As defined by (357 , the luminosity function 

w-if -^'.[(^r+(£)r'H(i;i.i.). (140) 

with normalization 

JLi/Lb 

As defined by }156| . the luminosity function 

r(^)", L,/Ai<L<L, 
Ygpw{L)^co{ * , (147) 

[ i-^)^ , L,<L< A2L, 

with normalization 

coi = a-i(i-z\r") + r'(4^-i). 

Generally, the luminosity function depends on z, so that Y{L) — Y{L; z). 

Luminosity density from the luminosity function 

The luminosity density £ = d£ /dVdt at redshift z, referred to comoving vol- 
ume, is given by 

i{z) = h{z){At) d{\nL)LY{L;z) , (148) 
Jo 

where n{z) is the comoving rate density and the integration over dL or dlnL 
depends on the definition of the luminosity function Y[L). The local lumi- 
nosity density, which depends on the low-redshift luminosity function and the 
local rate density, is 

/"OO 

io = no{At) / d{\nL) L Y{L) . (149) 
Jq 

The relation of mean duration (At) to luminosity is normalized so that the 
fluence J" = L{At). 

For the model of Ref. |156) . Y{L) is dimensionless (per unit InL), and the 
expression for the local luminosity density is analytic, given by 

^GPW = noCo{At)L4{l - Zir''+"V(1 + + (^2+^ - !)/(! + /?)] > (150) 
with = a-\l - Zij^") + /3-i(l - A-^). 
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For the model of Ref. [357], Y{L) has dimensions of L ^, and the local 
luminosity density requires a simple numerical integration of the expression 



klzvd 



rfii[(f)" + (f)1 > (151) 



with 



^0 ^ 



L-^ .In L 



dx +x"^)-i 
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Elementary Cosmology: Following, e.g., [357] . we write the differential rate 
of GRBs with redshift z between z and z -\- dz and luminosity L in the range 
L to L + di by 

dN ^ n{z)_ dViz) ^^^^ ^^53) 
dtdzdL 1 + z dz ' 
The comoving volume element 

dz Ho {i + zy^n^{i + z)3 + nA ' 

To derive this, note that dV^ — cdt^dA^, so dV^/dz — c\dt^/dz\dA^. With 
dA^ = {R^xfdn^ = dliz)dn/{l + z)^ (Ref. [3], eq. (4.43); see below). The 
cosmology of the universe we inhabit is a flat AC DM universe, that is, a 
dark-matter dominated flat universe with nonzero cosmological constant, well- 
described by 

?^ = , (154) 

dz Ho{l + z)^n„,{l + z)3 + SlA 



This defines, noting that V = {l + z)^V^, equation (153). Our standard values 
are Hq = 72 km s^^ Mpc"\ Q„, = 0.73 = 1- ^a- 

The proper distance, though not directly measurable, is the distance be- 
tween two objects that would be measured at the same time t. The proper 
distance at the present epoch is just the comoving coordinate, so 

dprop — X — ct ~ c dt = c dz' 1^1(1 + z') . (155) 
Jo Jo dz 

The energy fiux for a source isotropically radiating luminosity L ~ d£ /dt at 
proper distance dprop is related to the energy fiux from a source with isotropic 
luminosity = d£^,/dt^, at luminosity distance through the relation 

^ d£ d£/dt d£^,/dt^, (1561 
~ dAdt ~ Andl^^p ~ Andl 

The fluence J" = J^^ dt ${t) measured over some time interval ti - t2 is there- 
fore related to the apparent isotropic energy release through the expression 
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For an expanding universe, d£^ = e^,dN = e(l + z)dN = (1 + z)d£, and 
dt^ — dt/{l + z). From the definition of d^, the energy flux 

d£ _ _ dA d£ dU d£ + zfdA 

~ dAdt ~ Andl ^ Airdl ^dlfj ^~dt^ dAdt ^ ' ^ ' 

Therefore dA = R'^x^dH — 47rd|/(l + z)^. At emission time t*, dA^, = 
Rlx^dQ^, and dA = R'^x^df2, so that dA^/dA = 1/(1 + z) . This implies 
from the definition of the luminosity distance and the angular diameter i?*x, 
with d£ /dt = (1 + z)~'^{d£^,/dt^,), and letting df2 = dSl^,, 

R 

dL = dUz) = (1 + z)dprop = {-^)^{R*X) = (1 + zfdA = 

it* 

-^^0 Jo f2„i{l + z')^ + Ha Ho 4 

(159) 

= 4170z(l + 0.8z) Mpc for a flat universe. 

Compare this derivation from Ref. [3] . The directional event rate, or event 
rate per sr, is 

dN ^ f ,,r ■ , ~.,dU, f°° , ,dU AR.,xfnJz) 

dn = Tj'''* = 'd^' (i + z) ' ^'''^ 

where the burst emissivity n* (z) gives the rate density of events at redshift 
z. Volume density in comoving and proper coordinates is related by the ex- 
pression n*(e*; z) — hco{^*] z){l + z)^. Assuming separability of source emis- 
sion properties and the rate density, then nco,i{z) — hQiSi{z), where Si{z) 
is the structure formation history (SFH) of sources of type i, defined so that 
Si{z <C 1) = 1, and fiQi is the local {z <C 1) rate density of bursting sources 
of type i. Thus 



dn Ho Jo (1-f 2)3^/2™ (l + z)3 + i7A ' 

This expression can be easily generalized to accommodate spectral behavior. 

Star formation rate functions: The redshift-dependent rate density by Por- 
ciani & Madau (2001; SFR2) that has roughly constant star formation at z > 2 
|359j is described by 

. exp(3.4z) 

fipuiz) = 23no , . 162 

exp(3.4z) 4- 22 



The analytic function 
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Fig. 25. Star formation rate factors [5501135110371 . 



1 + ai 

(1 + z)-'^2 + ai{l + z)''3 ' 



(163) 



|360j . with fli = 0.005, 02 = 3.3, and 03 ~ 3.0 for the lower star formation rate 
(LSFR) and oi — 0.0001, 02 — 4.0, and 03 ~ 3.0 for the upper star formation 
rate (USFR) describes extreme ranges of optical/UV measurements without 
and with dust extinction corrections, respectively. 

The SFR history (SFR3) of Hopkins and Beacom (2006) [351], which is 
intermediate between the LSFR and USFR, is given by their analytic fitting 
profile 

1 + (022/01) 

where ai = 0.015, 02 = 0.10, 03 = 3.4, and 04 = 5.5 are best fit parameters 
(Figure [25]). An update of the SFH fit of [MI] is contained in Ref. [55^ in the 
form of a continuous broken power law. 



^hb(^) 



(164) 



p*(z) = Pa 



(1 



l + z 
B 



1 + 2 

c 



C7] 



-\ 1/') 



(165) 



where a = 3A{b = —0.3, c = 
last) piece breaking at zi 



-3.5) is the logarithmic slope of the first (middle, 
1 and Z2 = 4, and the normalization is po = 
0.02Mq Mpc~^ yr^^. Using rj — —10 smooths the transitions, while 77 -> cc 
recovers the kinkiness of the original form. Here B = {1 + zi)^^"/'' = 5000, 
and C = (1 + zi)(^-'^)/'=(l + 22)^"''/^ = 9. 

Fig.[25]shows star formation rate factors used in different models for GRBs. 
For data, see |363[ I360j . An enhanced rate of GRBs compared to the global 
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star formation rate is required to fit Swift statistical data |157j . as confirmed 
in other studies |362l 1364] . 



Luminosity function in Le &i Dermer model 



We approximate the temporally-evolving SED of, for example, a GRB or 
blazar by the expression 



= lit) - H{t;0,At) 



and 



S(x) = x^Hil - x) + x''H{:x -I) , X 



(166) 



(167) 



pk 



The bolometric fluence of the model GRB for observers with 6 < 6j is 



dt 



de^^h /,^, At , 



(168) 



where Xb is a bolometric correction to the peak measured lyF^ flux. If the SED 
is described by eq. (167), then Af, — (a~^ — b~^), and is independent of e^fc. 
The beaming-corrected 7-ray energy release £^,-y for a two-sided jet is 

J" 



f*^ = 47rrf^(l - fij) - 



+ z 



(169) 



where the luminosity distance is given in equation ( 159 1. Substituting eq. ( 168 ) 

f 



for into eq. ( 169 ) gives the peak flux 



47rd|(z)(l - fij)At^ Xb 



(170) 



Finally, by substituting eq. (1701 into eq. (166), the energy flux becomes 

lit) = l^,Hi^^■,f^,,l)Hit;0,At)Six). (171) 

The observed directional event rate for bursting sources with lyF^, spectral 
flux greater than fg at observed photon energy e is given by 



dNi> h) ^ c_ 

dn Ho 



dl 



de. 



diAt^) 



dz 



dliz)hcoiz)Yi£^y,lij,epk*,At^) 



S{l-t,Six) 



(172) 



(1 + z)3 ^f2mil + z)^ + f2A 
For discrete values of 5*^, epk*, and Z\t*, the property distribution function 
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(173) 

where g{^J■j) is the jet opening angle distribution, and fg is the instrument's 
detector sensitivity. The detector threshold for Swift and pre- Swift GRB de- 
tectors is taken to be ~ 10~^ and ^ 10"'' ergs cm~^ s~^, respectively. 

In the model of Le & Dermer (2007) |157) . the form for the jet opening 
angle g{nj) is chosen to be 



lifJ-j) = 30 (1 - fJ-jY -ff(Mi;Mjmin,/^jmax) , 



(174) 



where s is the jet opening angle power-law index; for a two-sided jet, /Xjmin ^ 0- 
Normalized to unity, equation (174) gives 

go = 



1 + s 

(1 - Mjmin)l+-' - (1 



Mjmax) 



1+s 



(175) 



Integrating over ^ in eq. (172) gives the factor (1 — ^j) describing the rate 
reduction due to the finite jet opening angle. Hence, eq. (172) becomes 



dN{> ft) ^ ^ 
df2 Ho 



dft 



.9(Mj)(1 - Aij ) 



dliz) n,,{z) 6 (f: - ft^,S{x)) 
X / dz , ^ ^ , (176) 

' (i + z)3 ^Q^ii +zr + f]A 



where /e^j. is given by eq. (170). The redshift, size, and opening angle distri- 
butions can then be calculated for comparison with data. 

Fitting to redshift and opening-angle distributions of pre-Swift and Swift 
GRBs using this model gives a good fit, assuming a flat vFi^ spectrum with 
bolometric factor Af, = 5, for the following parameters: 

AU = 10s, f = 4 X 10^° erg s"\ s = -1.3, fi,mn ^ 0.765, n„iax = 0.99875 , 

(177) 

i.e., 6l™„ ^ 0.05 (2.9°), 0max = 0.7 (40°). The "true" event rate hcRB = (7.5 
- 9.6) Gpc^^ yr^^. This can be compared with beaming corrected rates of 
= 75 X 0.5 Gpc"'^ yr~^ ~ 40 Gpc~^ yr~^ |156) . The local emissivity in all 
forms of electromagnetic radiation from long GRBs is automatically implied 
in this model, namely 

ho - tigrbS*^ « (3 - 4) X 10^3 erg s'^ Mpc'^ (178) 

Luminosity function: We derive the luminosity function for the Le & Der- 
mer (200 7) mo del, which can be checked in the limit 2; <C 1 to agree with 



equation (178). The isotropic luminosity i*.iso = 4:ndj^<P = Andj^^ fe/e. 



From equations (166) and (170) 
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/epfc = L^^,so(t)S{x)/\b , 

where 

L,^,,,{t) = L,^,,o Hifi;fij,l)H{t;0,At) = ^.^ H{^i:jij,l)H{t-Q, At) 

is the apparent isotropic luminosity that would be measured within the beam- 
ing cone during time t from to At. Thus the relation between the apparent 
luminosity L = L^,^iso and — £*^/Z\t* = 4 x 10^° erg s~^ is, obviously, 
L — Z*/(l — /ij) for a two-sided jet (0 < fij < 1). Thus dL/dfij — L/{1 — /ij), 
and 

noting the extra factor (1 — fij) which accounts for the smaller opening angle 
of the more apparently luminous jets. 
From the formulation of the model, 

-J^^ = n.g(M,) = n.,9o(l - ; M™"' /^f') • (180) 

Therefore 

dN _ .dfij d,N . (1-yUj) 



s+2 



dVdtdL ' dL ' dVdtdfi 



[1- ^j) - n^go H [fij , ^ ■ , /i ■ j 



/ I \ * + 3 / 7 1 



H \ L : ^— , . (181) 

7 \ r I \ " 1 ,,rmn' i ..max \ ^ ' 

The luminosity function in this model is 

'^^^^^^ ^ S ^ t (l) ^ "'^ '''' "'^ '"'^ ■ 

(182) 



Therefore L'^<P{L) cx L-'^-" cx i^-^. From equation ( |l75[ ), go = 0.051. 

Integrating <?(i) over L gives the emissivity according to equation (1491. 
We find that integ ratio n of equation (182) recovers the luminosity density 
given by equation (178 1. The luminosity function in the form niL^^LD(-^) is 
plotted in Figure |26| The model is overconstrained, which problem can easily 
be relaxed by changing the functional dependence of the jet opening angle. 



Local luminosity density of GRBs 

The "I'F^" intensity e/c(erg/cm^-s-sr) — mgC^ e^(/)(e)/47r for unbeamed sources, 
where (/)(e) is the differential photon flux, is given by 

TTieC^ f°° dt^ elnco{e*]z) 



47r 7o 1 + 2: 
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(eq. (4.57) in Rcf. [3^). This equation also applies to beamed sources provided 
that the apparent comoving event rate density nco{z) = de^nco{e*', z) is 
increased by the inverse of the beaming factor, /(,, to get the true source rate 



density (see eq. (Ill) 



The photon luminosity density at redshift z is given by 

/>oo 

£{z) = meC^ / (ie* e*rtco(e*; ^) (184) 
Jo 

and the local luminosity density £q = £{z ^ 1). Assuming that the only 
property of long-duration GRBs that is redshift dependent is the rate density, 
then £{z) — £qX!(z). The integrated intensity I = Jg dele is related to the 
mean photon flux ^ = d£/dAdt through the relation I = 'P/4tt, so 

£o^k , (185) 
c 

where 

dz , ^^''^ = . (186) 

ii + z)^^n,n{i + z)^ + nA 



Table 10 gives values for k obtained by numerically integrating eq. (186), 
using the SFR factors shown in Fig. [25] with the standard AC DM cosmology 
with Hubble constant Hq = 72 km/s-Mpc, f2„i = 0.27 and Ha = 0.73. The 
integration was truncated at a redshift z^ax ~ 10. 



Table 10. Constant k, from eq. (186 1, for different Star Formation Rate Factors 



SFR 



Constant comoving 1.85 
LSFR ,360, 0.52 
USFR [360] 0.13 
HB [361] 0.37 
LD [TS7] 0.33 - 0.37 



The mean flux <P over the full sky can be estimated from the total fluence 
per year from GRBs, given the number of GRBs per year. For 1293 GRBs in 
4B catalog the total fluence (20-300 keV) is 63.80 x lO"'^ erg/cm^ and the total 
fluence > 20 keV is 153.63 x 10"*^ erg/cm^ (M. Gonzalez, private communica- 
tion, 2003). For the 4th BATSE Catalog, considering 666 Burst/year fuU sky 
with 1293 bursts implies 1293/(666 GRB/yr)= 1.94 yr = 61225300 s, so that 
the mean 20 - 300 keV flux 1.04 x 10"^° erg/cm^-s or 3.3 x 10"^ erg/cm^-yr. 
Likewise, for the total > 20 keV fluence was, we obtain 2.54 x lO^^*' erg/cm^-s 
or 7.9 X 10-3 erg/cm^-yr Band (2002) [5S5] obtains 550 GRBs/yr for BATSE 
exposure, which brings the flux down to 6.5 x 10^'' erg/cm^-yr. Subtracting 
10 - 20% for short GRBs gives (5 - 6) x 10"^ erg/cm^-yr. 
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Fig. 26. Long duration GRB luminosity function, plotted in the form of a differential 
luminosity density £ = hoLY{L)At, for various models. 



From eq. (185), writing <P^2 = ^/(lO ^ crg/cm -yr) wc find £o — 2.3 



lO-*-^ fc<?_2 erg/Mpc-'-yr. With <P_2 ~ 0.5 and k = 0.35, then 



in 



0.5 X 10« (-^)(^) erg/MpcV 



(187) 



in agreement with the results of Eichler et aL (2010) |366j . 

This value is at least a factor of 5 smaller than values obtained through 
statistical treatments of GRB data |157l 11561 1357) , as can be seen from Table 
11 Here we give the local luminosity densities for long GRB models derived 
either numerically or analytically from the luminosity functions derived in 
these models. 



Table 11. Local luminosity density of long GRBs for different models 



Model 


4 


Photon energy range 




(10^* erg/Mpc^-yr) 




Schmidt (2001) [363 


0.03 


50-300 keV 


Guetta et al. (2005) [156, 


0.05 - 0.08 


50-300 keV 


Le & Dermer (2007) tl57j 


0.30 - 0.38 


all photon energies 


Liang et al. (2007) [357] 


2.0 


1 keV - 10 MeV 


Wanderman & Piran (2010) '364' 


1.3" 


1 keV - 10 MeV 



"Normalizing to the Swift durations brings this value down by a factor of « 2. 



The reason for the discrepancy between the local emissivity of long du- 
ration GRBs derived from BATSE data and inferred from Swift data is sug- 
gested by Fig. 26 This figure compares different models for the local GRB 
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Fig. 27. Local luminosity density requirements on different models for UHECR pro- 
duction from long-duration GRBs with rate density following various star-formation- 
rate functions shown in Fig. [25] 



luminosity functions, plotted in the form of a differential emissivity function 
'hQAtLY{L), where fiQ is the local rate density, and the characteristic dura- 
tion of a GRB in the explosion frame is At = 10 s. Because the treatments of 
Schmidt (2001) [M] and Guetta, Piran, & Waxman (2005) [HI] are based on 
BATSE observations in the 50 - 300 keV range, a bolometric correction factor 
= 6.3 is applied to normalize them to the 1 keV - 10 MeV range (D. Wan- 
derman, private communication, 2010). In this case, the implied luminosity 
density from the treatments of the luminosity function agree within factors of 
3. 

This figure show that the integrated luminosity density depends sensitively 
on the maximum luminosity taken in the construction of the luminosity func- 
tion. When redshift information is available, very luminous GRBs are found 
with apparent luminosities reaching 10^'^ erg/s. To fit these luminous GRBs, 
the luminosity function must extend to large luminosities, and these very rare, 
very luminous GRBs make a significant contribution to the integrated lumi- 
nosity density. For such luminosity functions, the large fluence GRBs make up 
a large, even dominant part of the total GRB fluence, but such high fluence 
events are so rare that the chances to now of such a GRB with large, ^ 10~^ 
erg/cm^, fluence having been detected is low. 

The local electromagnetic emissivity of long GRBs, given in Table 11 is 
therefore found to be < lO'"' erg/yr-Mpc^ [368! • This can be compared with 
the minimum local emissivity of (uhecr ~ 4 x 10**^ erg/s- Mpc required to 
power UHECRs, as seen in Fig. [27| and discussed in more detail in the next 
section. Related estimates can be made for other source classes |357l I369j . 
including the short hard GRB class, the X-ray flashes, the sub-energetic GRBs, 
and engine-driven supernovae with relativistic outflows that lack GRB-type 
emissions. These classes, though they have much larger local rate densities 
than long-duration GRBs, are not much preferred over long-duration GRBs 
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on the basis of their local luminosity density. The so-called low-luminosity 
GRBs, for example, have a large luminosity density in ejecta kinetic energy, 
^ 10^^ erg Mpc~'^ yr~^ |370l 1371] . but smaller or comparable nonthermal 
7-ray luminosity densities than long-duration GRBs. 

8.3 Closure relations 

The relativistic blast wave theory for the afterglow has become the industry 
standard model. Familiarity with the derivation of the closure relations is 
required for basic knowledge of GRB physics. The reader may see |^ for the 
derivation, [3j for review, and (372j for the essential treatment. 

Interpreting the delayed onset of the Fermi-LAT radiation as due external 
shock emission, one approach is to suppose that the blast wave decelerates 
adiabatically in a uniform surrounding medium |373j . with closure relation 
uFi, oc t(2-3p)/4j^(2-p)/2^ where p is the electron injection index. A value of 
p « 2.5 gives a reasonable fit to the data. Another regime to consider is 
a radiative GRB blast wave |356j . where the comparable closure relation is 
vF^ cx i(2-6p)/7j,(2-p)/2^ ^j^j^ p « 2 giving a plausible fit to the data. The 
adiabatic model requires unusually low densities and magnetic fields for GRB 
080916C, and the radiative model supposes pair loading can help achieve 
strong cooling. 

Alternate leptonic models for Fermi LAT GRBs include photospheric mod- 
els with the photospheric emission passing through shocked plasma in the col- 
liding shells or external shocks of the GRB outfiow '374'. A joint Fermi-Swift 
paper examines leptonic afterglow models for GRB 090510 [358] . 



9 Fermi Acceleration, Ultra-High Energy Cosmic Rays, 
and Black Holes 

Fermi acceleration is intrinsically coupled to the presence of a magnetic field 
that governs particle motions. The energy is gained stochastically in second- 
order Fermi scenarios, for example, via resonant pitch-angle scattering, or the 
particle energy is gained systematically by those few particles that diffuse back 
and forth across a shock front multiple times before convecting downstream 
(Section 5.5). In either case, the energy gain rate is related to the Larmor 
rate v/ri^. The magnetic field and the size scale set a fundamental limit on 



the maximum particle energy through the Hillas criterion, eq. ( 78 ) 



9.1 Maximum particle and synchrotron photon energy 

In Fermi acceleration scenarios, the maximum particle energy is obtained by 
balancing acceration rate with the energy-loss rate. For relativistic electrons, 

heseweb.nrl.navy.mil/gamma/ ~dermer/notes/index.htm 



110 Charles D. Dermer 



the maximum comoving electron Lorentz factor is ^'max 

= \/67re/(^CTTS'. The 

maximum synchrotron photon energy for a relativistic jet source at redshift 



z is, from eq. ( 145 ) 



syn,max — „ , n j (l^^) 

8af 0(1 + z) 
or 

Esyn,max = 240 j^^^ GeV . (189) 

The delayed appearances of the highest energy photons from Fermi LAT 
GRBs, which typically happens at times after the GBM igo, calls into question 
a synchrotron origin for these photons, unless F remains large at the end of 
the prompt phase of GRBs. 

The maximum energy for protons and ions is obtained from the Hillas 
criterion, written as — E'//3QB' < R' . The maximum escaping particle 
energy E^ax = rE'^^^ < FE' = ZfieB'R, noting that R' = R/F from length 
contraction of the stationary frame size scale as measured in the comoving 
frame. Relating the magnetic field energy density by a factor es < 1 times 
the proper frame energy density associated with the wind luminosity L, then 
= 2eBL/{l3cR^F^), and 



^ F 

implying 



E^^x < ( ^ ) J'-^ , (190) 



L, > F' El, erg/s (191) 



'3 X 1045 

|3751 13761 1308] , noting that the apparent 7-ray luminosity < L 



Eqs. ([721 and ([74| imply 

, A F ™, / 27 \ A{F/ WOO) 



E^onsyn, max - m,C ^ ^ ^ U«/ / " '^'^'^ ^Z^l + z) ^""^ ' 

(192) 

Thus proton or ion synchrotron can make a 7-ray component in jet sources 
provided the ion power and energy is sufficiently great. 



9.2 L-F diagram 

Using the minimum Lorentz or Doppler factors implied from 7-ray opacity 
arguments, eqs. (45) and (46), along with the apparent 7-ray luminosity cor- 
responding to the time that Fmin is measured, we can construct an L-F dia- 
gram. 

Figure 28 shows such a plot for blazars and GRBs |308) . Note that GRB 
090510A is a short hard GRB, whereas GRB 080916C and GRB 090926A are 
long-duration GRBs. Except for PKS 2155-304, which uses HESS data for the 
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Fig. 28. Data shows apparent isotropic L versus Fmin for blazars, radio galaxies, 



and GRBs. Solid and dot-dashed curves plot the constraint given by eq. ( 191 1 



giant outbursts in 2006, all the data were measured with the Fermi LAT. As 
can be seen, the powerful Fermi GRBs have more than adequate luminosity 
to accelerate either protons or ions to ultra-high energies, even with Lorentz 
factors Fmin ~ 10^ as inferred from 77 opacity arguments. By comparison, the 
blazars and radio galaxies have smaller apparent luminosities and also smaller 
Fmin- Acceleration of high-Z ions like Fe to ultra- high energies is possible for 
these sources on the basis of eq. (190), but acceleration of protons appears 
unlikely, except possibly during flaring episodes. 



9.3 Luminosity density of extragalactic 7-ray jet sources 

Maintaining the intensity of UHECRs against photopion losses with the CMB 
requires a luminosity density in UHECRs of lO''^ erg/Mpc'^-yr |377| . This 
value can easily be obtained by dividing the energy density of UHECRs at ~ 
10^° eV, which is ^ 10^^^ erg/cm^, by the photopion loss timescale ~ 100 
Mpc/c. Because increases rapidly with decreasing energy, the required 
luminosity density remains at this level even when fitting to the ankle of 
the cosmic-ray spectrum at E ^ 4 x 10^* eV |378) . though precise values of 
the local luminosity density depends in detail on the assumed cosmic star- 
formation rate factor applicable to GRBs. 
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The standard argument for inferring the luminosity density of a source 
class is to assume that the nonthermal 7-ray luminosity represents, to an 
uncertain factor of order unity, the luminosity in UHECRs. This seems rea- 
sonable because some radiative losses into 7 rays are likely to occur during 
acceleration. If those losses are small, then the UHECR power could exceed 
the 7-ray luminosity, but if the particles experience severe radiative losses 
during acceleration, then they could not be the sources of the UHECRs |377j . 
This estimate is furthermore independent of beaming factor, because for every 
beamed source we detect, a proportional number of misdirected sources will 
point away from us. The argument suffers, however, from the likelihood that 
a large percentage of the 7 radiation is produced by ultra-relativistic leptons 
(especially in the case of blazars where correlated variability with lower en- 
ergy bands can be monitored). Moreover, it is not certain that MeV 7 rays 
from long duration GRBs, which comprise the bulk of the energy output, is 
entirely nonthermal, inasmuch as a thermal/photospheric interpretation can 
potentially resolve the line-of-death problem that plagues nonthermal syn- 
chrotron interpretations of GRBs |379j . Thus the nonthermal LAT emission 
might be a more appropriate luminosity with which to define the nonthermal 
luminosity density of GRBs [366] . 
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Fig. 29. 7-ray luminosity density inferred from Fermi observations of various classes 
of 7-ray emitting AGNs |308j , compared with the local luminosity density of long- 
duration GRBs (Table 111, and UHECR source requirements. 
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The luminosity density of GRBs has been evaluated in detail in Section 
8.2. The luminosity density of radio galaxies and blazars, based on the average 



100 MeV - 100 GeV flux measured in the ILAC [371, is shown in Fig. 29 



9.4 Origin of UHECRs 

The sources of the UHECRs are unknown. With the anisotropy of the UHE- 
CRs observed with Auger becomes less significant that originally reported, 
leaving the cluster of UHECRs in the vicinity of Cen A as the only out- 
standing hotspot [6], even Galactic sources have been considered as plausible 
UHECR candidates [380j . Indeed, a wide variety of possible source classes 
could contribute to the UHECRs, including magnetars and young pulsars 
|3811 1382| , particle acceleration in structure-formation shocks j383| , and new 
physics candidates. Criteria that the sources of UHECRs should meet are [384] 

• extragalactic origin; 

• mechanism to accelerate to ultra-high energies; 

• adequate luminosity density; 

• sources within the GZK radius; and 

• UHECR survival during acceleration, escape, and transport. 

On this basis, we [3j suggest that the very luminous, very energetic extragalac- 
tic blazar and GRB jet sources are the most plausible candidates, with the 
BL Lac objects and FRl radio galaxies perhaps favored [308] . Electromagnetic 
signatures of ultra-relativistic hadrons in GRBs and blazars are not clearcut. 
High-energy neutrino detections from transient or bursting 7-ray sources will 
be crucial to finally solve this puzzle. 



9.5 Black holes, jets, and the extreme universe 
Black hole physics 

The Kerr metric for an uncharged rotating black hole, written in the Boyer- 
Lindquist coordinates {i, r, 9, ip} useful for asymptotic analysis in the exterior 
region of the horizon of the black hole, takes the form [385] 

ds^ = {/3^ - a^) dt^ + 2/3^d(pdt + dr^ + jge dO^ + 7<^<^ dip'^ , (193) 

where the metric coefficients are given by 



2Mr 



_ _ -2Mrasin^6' _ 
Pip = 9tip — 2 ' "^^^ — ~A 



S'^ sin^ ( 



jge = p , and - 
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Here, a is the lapse function and /3 is the shift vector, 



2 2,2 2 , 

p — r +a cos I 



A 



^Mr + a^ , and ^ {r^ + a^f - A siii^ 



AdditionaUy 



— ^ , and 



a ^7 — sin ( 



The parameter a is the angular momentum per unit mass of the Kerr black 
hole, so aM is the angular momentum of the black hole. The horizons are 
located at r± = Af ± V^M^^^o^, as shown in Fig. 30 



0.0 

30.0 330.0 




a = 0.999 150.0 210.0 
180.0 



Fig. 30. Sketch of locations of the horizons r+ and r_ and the ergosphere TergiP) 
in the Kerr geometry. 



The set of points described by gu — defines the ergosphere 



rerg{0) = M+ A'P - a? cos^ I 
Note that gu > in the region 

r+ < r < rergiO) , 

and consequently i = dt becomes spacelike in the above region. When r > 
t is timelike and future directed. See Fig. 



30 



(194) 
(195) 

'erg. 
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In elementary electrodynamics j386l , the Poynting vector in flat space far 
away from a gravitating object is given by 

S^^iExH), (196) 



and 



where 



'^^tot = field + ^mech) , (197) 

£f^eld = ^ / d^x {E^ + B^) , and (198) 

OTT Jy 

^Sraech = f X J ■ E (199) 



V 



give the field and mechanical energies, respectively. Thus we see that J ■ E is 
an emissivity, so conservation of energy requires 

'^ + W-S^-J-E, (200) 
ot 

where the field energy density is u — {E ■ D + B ■ H)/%tt. 

Electromagnetic energy extraction through the Blandford-Znajek process 
involving Penrose processes in curved space can be expressed in the case of a 
stationary axisymmetric force-free magnetosphere as 

'^'^ 'bz = S"-^. = -H^UB'-^r (201) 



dAdt 

|3881 1387] . The fields must satisfy the Znajek regularity condition |389j 

2 



SI n' 

r+ = B' {2Mr n-a) 



(202) 



at the event horizon. We can suppose that the polodial function f2 defining 
the fields is radially independent, so 



^Ibz = 47rA/r+ d0^ ^—f2{f2H - mt , (203) 
dt Jq a 

where B!^ is the radial magnetic field at the event horizon and 



= 2M^ 

is the angular velocity of the event horizon |390| . 

The form of J7 in eq. (2031 is found by solving the constraint equation, 
originally derived by Blandford and Znajek [391] . governing fields and currents 
in a black-hole magnetosphere under the force-free condition. The constraint 
equation in the 3-1-1 formalism takes the form [387j 
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1 dH^ 

where the charge density pc and current density J is expressed in terms of 
f2 and metric coefficients. By satisfying the Znajek regularity condition, eq. 



( 202 ) , we |387j obtained solutions 



{p^ = r'^ + cos^ 9) by considering the behavior of the constraint equation 
in the far-field limit for radially-independent Q. Thus 

- 2AdTpi ' ^'''^ 

and 

^ — ^ ■ (208) 
asm tl 

The i7+ solution implies a rate of electromagnetic energy extraction of 

dS, nQW a a\ Sn.^.ii^T^ a/2M«l 
-i7 BZ = — [ arctan 7^ ) ^ -i^B^M^x{ 

(209) 

after relating Qq and Bq fS]. The f2+ solution generalizes the Blandford-Znajek 
|391j split monopole solution for all values of < a < M. Inspection of this 
solution shows that Poynting flux outflow is greatest along the equator, so the 
energy flux has a pancake geometry. The jet geometry, as found ubiquitously 
in nature, is not described by this solution. 

The exact i7_ solution implies an electromagnetic energy flux peaking to- 
wards the poles of the rotating black hole, but describes an inward energy 
flux, so is unphysical. The J7_ solution has an additional freedom in the as- 
signment of an arbitrary poloidal field A in the prescription for the poloidal 
magnetic field 

Bp = ^{-[2gdr + f2rdg) (210) 

|392| . Energy extraction is accomplished not only by Poynting fiux outflow, 
but also by the matter currents that sustain the black-hole magnetosphere. 
Analysis of the timelike geodesies consistent with the Znajek regularity con- 
dition [393], and symmetry transformation of the i7_ solution yields a dual 
class of exact solutions with positive energy extraction of Poynting and matter 
outflows. The radial magnetic field for this dual class of solutions takes the 
form 

Br = Ia^^ , (211) 
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The extra freedom in the solutions with general 0-dependent A provides formal 
solutions devoid of strict physical content, and with arbitrary a dependence. 



Performing dimensional analysis of energy extraction using eq. ( 203 1 (see 
[3871 13941 [5] for detailed treatments) yields an estimate of the BZ power given 

by 

Ibz ~ 4nMr+ ^Q^Qh - [2)3 f w TrMr+ ^/2|fBf , (212) 



dt 



a 



where is the radial component of the magnetic field threading the event 
horizon. Note that the expression peaks at = SIh/'^ [395j . Along the equato- 
rial direction, 6 = n/2, ^/y/a = U^/A. Supposing the expression is evaluated 
on a size scale somewhat larger than the event horizon, to cancel the diver- 
gence, we have 



f Ibz - nc {^f r\ « 10^^ (^)' Af| Bj erg/s , (213) 



where Mg is the mass of the black hole expressed in units of 10^ M©, and the 
magnetic field is given in units of 10* G. 

As expressed, and often found in the literature, the value of 10* G is 
artificial. The magnetic field threading the ergosphere can be scaled to the 
Eddington luminosity 



by the expression 



^Edd = 1-26 X lO^^Afg erg/s (214) 



= es«Edd -. — 2~ ' (.215) 



where the parameter .^Edd refiects not only different accretion rates but dif- 
ferent radiative efficiencies, including the reduction in radiative efficiencies in 
the low Eddington ratio, advection-dominated regime. Thus 



« es^Edd « 1 . (216) 

-t^Edd 



Jets and the extreme universe 



The time-averaged luminosities of FSRQs extend to values in excess of « 
lQi9 gj.g g-i (Figure [tJ. In the extraordinary 2010 November flare of 3C 454.3 
[30] . reached apparent isotropic luminosity of (2.1 ± 0.2) x lO^*' erg s^^ 
over a period of a few hours, making it the most luminous blazar yet observed. 
Black-hole mass estimates for 3C 454.3 are in the range 0.5 ;S Mg ^ 4 |103j . 
where IO^MqMq is the mass of the black hole powering this AGN. For this 
range of masses, the Eddington luminosity therefore ranges from « 6 x 10*^ 
erg s"^ to « 5 X 10*^ erg s~^. During this extreme outburst, the apparent 
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luminosity of 3C 454.3 was more than a factor of « 400 greater than its 
Eddington luminosity. Even its time-averaged luminosity of L-^ 10"^^ erg 
s~^ is super-Eddington by a factor of « 20. 

Assuming that the Eddington condition does limit accretion flow onto the 
black hole, which is likely to be the case for the long-term average luminosity 
if not for the flaring luminosity, then the absolute radiant luminosity is lim- 
ited to a value of Labs ~ 5 x lO^'' erg s^^. This is consistent with the large 
apparent luminosities if the emission is highly beamed. For a simple top-hat 
jet beaming factor, a jet opening angle 6j implies a beaming factor = ^|/2 
for a two-sided jet with ^ 1. A mechanism for coUimation is, however, 
required. Should this arise from the Blandford-Znajek process, then we are 
still restricted to values of the absolute Blandford-Znajek power, eq. (213). As 
we have seen, by scaling the energy density of the magnetic field to the energy 
density of accreted matter near the event horizon shows that the Blandford- 
Znajek power is likewise Eddington-limited. Making the hypothesis [396' that 
the extraction of energy through black-hole rotation collimates the jet outflow 
with COS0J sa a/M, then /& = 1 — (a/M) and 9j = a/2(1 — a/M), implying 
a/M > 1 — (LEdd/-^^7)- If the jet opening angle is a consequence of the bulk 
Lorentz factor F of the outflow, then F ^ •\/L^/2iEdd- 

For the case of 3C 454.3 in its flaring state, when L^/Lgdd ^ 10"^, this 
hypothesis then implies that a/M > 0.999 and F /i^ 23, which is consistent 
with the value Fmin ~ 14 from 77 opacity arguments |30j . In the most con- 
servative case with Mg — 4, a/M ~ 0.998 and F }t 15. This is marginally 
consistent with the limiting maximum value a/AI = 0.996 suggested by As- 
chenbach from analyses of microquasars and the Galactic Center black hole 
|397j . The smaller black- hole mass estimate, Mg = 0.5 |103j . implies a value 
of a/M that violates this limit by a large margin. Much work, both numerical 
and theoretical, has been devoted to jet formation from the Blandford-Znajek 
process, and it is unclear if jet coUimation can be described by the guess that 
cos6'j « a/M, but it is interesting to suggest this possibility, which leads to 
values of F consistent with separate inferences regarding the outflow Lorentz 
factor. 

Figure 31 summarizes some various source classes and their residence in a 
diagram of L vs. time. The extreme universe, as contrasted with the moder- 
ate, or relaxed universe, is deflned by eq. ( 14 ) , when the ratio of the apparent 
luminosity and variability timescale exceeds the Eddington luminosity divided 
by the light-crossing time for the Schwarzschild radius, L-Edd/ts = 1-26 x 10^'^ 
erg/s^. "BH*" stands for events where normal main sequence stars are tidally 



captured by a black hole, as in the 28 March 2011 event [398] . and "EMRI" 
stands for extreme mass-ratio inspiral events [399]. The luminosity limit of 
the extreme universe is given by eq. (171, namely C'/G ~ 3.63 x 10^^ erg/s. 
Through observations of blazars and GRBs, 7-ray astronomy is pushing to- 
ward this limit. A coordinated assault involving gravitational wave observa- 
tories, the Fermi 7-ray telescope, neutrino observatories and multiwavelength 
campaigns might over the next decade reach this limit of the extreme universe. 
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Fig. 31. Luminosity versus characteristic timescales for variability or duration for 
various types of astronomical sources and cosmic events. Upper limit is from eq. ( 17 1 
and the horizontal line separating the extreme and moderate universe is defined by 
eq. ([Ml. 
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